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Abstract. We consider the fractional Laplace framework and provide mod¬ 
els and theorems related to nonlocal diffusion phenomena. Some applications 
are presented, including: a simple probabilistic interpretation, water waves, 
crystal dislocations, nonlocal phase transitions, nonlocal minimal surfaces and 
Schrodinger equations. Furthermore, an example of an s-harmonic function, 
the harmonic extension and some insight on a fractional version of a classical 
conjecture formulated by De Giorgi are presented. Although this book aims 
at gathering some introductory material on the applications of the fractional 
Laplacian, some proofs and results are original. Also, the work is self con¬ 
tained, and the reader is invited to consult the rich bibliography for further 
details, whenever a subject is of interest. 
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CHAPTER 1 


Introduction 


In the recent years the fractional Laplace operator has received much attention 
both in pure and in applied mathematics. 

The purpose of these pages is to collect a set of notes that are a result of 
several talks and minicourses delivered here and there in the world (Milan, Cortona, 
Pisa, Roma, Santiago del Chile, Madrid, Bologna, Porquerolles, Catania to name a 
few). We will present here some mathematical models related to nonlocal equations, 
providing some introductory material and examples. 

Starting from the basics of the nonlocal equations, we will discuss in detail some 
recent developments in four topics of research on which we focused our attention, 
namely: 

• a problem arising in crystal dislocation (which is related to a classical 
model introduced by Peierls and Nabarro), 

• a problem arising in phase transitions (which is related to a nonlocal ver¬ 
sion of the classical Allen-Cahn equation), 

• the limit interfaces arising in the above nonlocal phase transitions (which 
turn out to be nonlocal minimal surfaces, as introduced by Caffarelli, 
Roquejoffre and Savin), and 

• a nonlocal version of the Schrodinger equation for standing waves (as in¬ 
troduced by Laskin). 

This set of notes is organized as follows. To start with, in Chapter 2, we will 
give a motivation for the fractional Laplacian (which is the typical nonlocal operator 
for our framework), that originates from probabilistic considerations. As a matter 
of fact, no advanced knowledge of probability theory is assumed from the reader, 
and the topic is dealt with at an elementary level. 

In Chapter 3, we will recall some basic properties of the fractional Laplacian, 
discuss some explicit examples in detail and point out some structural inequalities, 
that are due to a fractional comparison principle. This part continues with a quite 
surprising result, which states that every function can be locally approximated by 
functions with vanishing fractional Laplacian (in sharp contrast with the rigidity 
of the classical harmonic functions). We also give an example of a function with 
constant fractional Laplacian on the ball. 

In Chapter 4 we deal with extended problems. It is indeed a quite remarkable 
fact that in many occasions nonlocal operators can be equivalently represented as 
local (though possibly degenerate or singular) operators in one dimension more. 
Moreover, as a counterpart, several models arising in a local framework give rise to 
nonlocal equations, due to boundary effects. So, to introduce the extension problem 
and give a concrete intuition of it, we will present some models in physics that are 
naturally set on an extended space to start with, and will show their relation with 
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the fractional Laplacian on a trace space. We will also give a detailed justification 
of this extension procedure by means of the Fourier transform. 

As a special example of problems arising in physics that produce a nonlocal 
equation, we consider a problem related to crystal dislocation, present some math¬ 
ematical results that have been recently obtained on this topic, and discuss the 
relation between these results and the observable phenomena. 

Chapter 5, 6 and 7 present topics of contemporary research. We will discuss in 
particular: some phase transition equations of nonlocal type, their limit interfaces, 
which (below a critical threshold of the fractional parameter) are surfaces that 
minimize a nonlocal perimeter functional, and some nonlocal equations arising in 
quantum mechanics. 

We remark that the introductory part of these notes is not intended to be 
separated from the one which is more research oriented: namely, even the chapters 
whose main goal is to develop the basics of the theory contain some parts related 
to contemporary research trends. 

Of course, these notes and the results presented do not aim to be comprehensive 
and cannot take into account all the material that would deserve to be included. 
Even a thorough introduction to nonlocal (or even just fractional) equations goes 
way beyond the purpose of this book. 

Many fundamental topics slipped completely out of these notes: just to name 
a few, the topological methods and the fine regularity theory in the fractional cases 
are not presented here, the fully nonlinear or singular/degenerate equations are 
not taken into account, only very few applications are discussed briefly, important 
models such as the quasi-geostrophic equation and the fractional porous media 
equation are not covered in these notes, we will not consider models arising in game 
theory such as the nonlocal tug-of-war, the parabolic equations are not taken into 
account in detail, unique continuation and overdetermined problems will not be 
studied here and the link to probability theory that we consider here is not rigorous 
and only superficial (the reader interested in these important topics may look, for 
instance, at [ 15 , 100 , 101 , 99 , 32 , 96 , 37 , 14 , 54 , 19 , 18 , 16 , 127 , 10 , 71 , 72 , 
133 , 8 , 21 , 20 , 38 , 39 , 42 , 67 , 89 , 98 , 113 , 66 ]). Also, a complete discussion 
of the nonlocal equations in bounded domains is not available here (for this, we 
refer to the recent survey [ 119 ]). In terms of surveys, collections of results and 
open problems, we also mention the very nice website [2], which gets^ constantly 
updated. 

Using a metaphor with fine arts, we could say that the picture that we painted 
here is not even impressionistic, it is just nai'f. Nevertheless, we hope that these 
pages may be of some help to the young researchers of all ages who are willing to 
have a look at the exciting nonlocal scenario (and who are willing to tolerate the 
partial and incomplete point of view offered by this modest observation point). 


^It seems to be known that Plato did not like books because they cannot respond to questions. 
He might have liked websites. 



CHAPTER 2 


A probabilistic motivation 


The fractional Laplacian will be the main operator studied in this book. We 
consider a function u: M" —> K (which is supposed^ to be regular enough) and a 
fractional parameter s e (0,1). Then, the fractional Laplacian of u is given by 


( 2 . 0 . 1 ) 


(-A)®u(a:) 


C{n,s) r 2u{x) 
2 


u{x + y) — u{x — y) 

|y|n+2s 


dy, 


where C{n,s) is a dimensional^ constant. 

One sees from (2.0.1) that (—A)® is an operator of order 2s, namely, it arises 
from a differential quotient of order 2s weighted in the whole space. Different 
fractional operators have been considered in literature (see e.g. [40, 128, 111]), 
and all of them come from interesting problems in pure or/and applied mathematics. 
We will focus here on the operator in (2.0.1) and we will motivate it by probabilistic 
considerations (as a matter of fact, many other motivations are possible). 

The probabilistic model under consideration is a random process that allows 
long jumps (in further generality, it is known that the fractional Laplacian is an 
infinitesimal generator of Levy processes, see e.g. [13, 7 ] for further details). A more 
detailed mathematical introduction to the fractional Laplacian is then presented in 
the subsequent Section 3.1. 


2.1. The random walk with arbitrarily long jumps 

We will show here that the fractional heat equation (i.e. the “typical” equation 
that drives the fractional diffusion and that can be written, up to dimensional 
constants, as dtu + (—A)®u = 0) naturally arises from a probabilistic process in 
which a particle moves randomly in the space subject to a probability that allows 
long jumps with a polynomial tail. 

For this scope, we introduce a probability distribution on the natural num¬ 
bers N* := {1,2, 3, • • • } as follows. If / c N*, then the probability of I is defined to 
be 


^To write (2.0.1) it is sufficient, for simplicity, to take here u in the Schwartz space of 

smooth and rapidly decaying functions, or in C^(R^) n We refer to [131] for a refinement 

of the space of definition. 

^The explicit value of C{n, s) is usually unimportant. Nevertheless, we will compute its value 
explicitly in formulas (3.1.10) and (3.1.15). The reason for which it is convenient to divide C(n, s) 
by a factor 2 in (2.0.1) will be clear later on, in formula (3.1.5). 
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2. A PROBABILISTIC MOTIVATION 


The constant Cg is taken in order to normalize P to be a probability measure. 
Namely, we take 



so that we have P(N*) = 1. 

Now we consider a particle that moves in K" according to a probabilistic process. 
The process will be discrete both in time and space (in the end, we will formally take 
the limit when these time and space steps are small). We denote by r the discrete 
time step, and by h the discrete space step. We will take the scaling t = and 
we denote by u(x, t) the probability of finding the particle at the point x at time t. 

The particle in K" is supposed to move according to the following probabilistic 
law: at each time step r, the particle selects randomly both a direction v e dBi, 
according to the uniform distribution on dBi, and a natural number k e N*, ac¬ 
cording to the probability law P, and it moves by a discrete space step khv. Notice 
that long jumps are allowed with small probability. Then, if the particle is at time t 
at the point xg and, following the probability law, it picks up a direction v e dBi 
and a natural number k e N*, then the particle at time t + t will lie at Xg + khv. 

Now, the probability u(x, f + r) of finding the particle at x at time t + r is the 
sum of the probabilities of finding the particle somewhere else, say at x + khv, for 
some direction v G dBi and some natural number k e N*, times the probability of 
having selected such a direction and such a natural number. 


(xo + khv, t + r) 



h 


Figure 1. The random walk with jumps 


This translates into 

u{x, t + t) 


|5Pi 


sL 

km* 


u{x + khv,t) 


\k\ 


l+2s 


dn^-\v). 
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Notice that the factor Cs/\dBi\ is a normalizing probability constant, hence we 
subtract u{x, t) and we obtain 


u{x, t + t) — u{x, t) 




u{x + khv,t) 


\k\ 


l + 2s 


dH^ ^{v)—u{x,t) 


u{x + khv, t) — u{x, t) 


\k\ 


l + 2s 




As a matter of fact, by symmetry, we can change v to —v in the integral above, so 
we find that 


u{x, t + t) — u{x, t) 


Cs f U{x-khv,t) -U{x,t) ^ 


Then we can sum up these two expressions (and divide by 2) and obtain that 
u{x, t + t) — u{x, t) 

u(x + khv.t) + u(x — khv.t) — 2u(x,t) , „ i, , 

dn {v). 


2\dBi 


kem 


\k\ 


l + 2s 


Now we divide by r = we recognize a Riemann sum, we take a formal limit 
and we use polar coordinates, thus obtaining: 


dtu{x,t) 


i{x, t + t) — u{x, t) 


Cs h 

^ I JdBi 


'2\dBi 
c 


- s f 

1 kem 

p + GO r 

Jo <JdBi 

-f 

l| Jr' 


u{x + khv, t) + u{x — khv, t) — 2u(x, t) , 


\hk\ 


l+2s 


-dH'^-\v) 


u{x + rv,t) + u{x - rv, t) - 2u{x, t) 


l + 2s 

u{x + y,t) + u{x — y,t) — 2u{x, t) 


2\dBi 

= - Cn,s {-^yu{x,t) 


\y\ 


n+2s 


dy 


for a suitable Cn,s > 0. 

This shows that, at least formally, for small time and space steps, the above 
probabilistic process approaches a fractional heat equation. 


We observe that processes of this type occur in nature quite often, see in par¬ 
ticular the biological observations in [140, 91], other interesting observations in 
[118, 126, 142[ and the mathematical discussions in [94, 85, 110, 105, 108[. 

Roughly speaking, let us say that it is not unreasonable that a predator may 
decide to use a nonlocal dispersive strategy to hunt its preys more efficiently (or, 
equivalently, that the natural selection may favor some kind of nonlocal diffusion): 
small fishes will not wait to be eaten by a big fish once they have seen it, so it may 
be more convenient for the big fish just to pick up a random direction, move rapidly 
in that direction, stop quickly and eat the small fishes there (if any) and then go 
on with the hunt. And this “hit-and-run” hunting procedure seems quite related to 
that described in Figure 1. 
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2.2. A payoff model 


Another probabilistic motivation for the fractional Laplacian arises from a pay¬ 
off approach. Suppose to move in a domain Q according to a random walk with 
jumps as discussed in Section 2.1. Suppose also that exiting the domain 17 for the 
first time by jumping to an outside point y e K"'\17, means earning uo{y) sestertii. 
A relevant question is, of course, how rich we expect to become in this way. That 
is, if we start at a given point x G 17 and we denote by u{x) the amount of sestertii 
that we expect to gain, is there a way to obtain information on it? 

The answer is that (in the right scale limit of the random walk with jumps 
presented in Section 2.1) the expected payoff u is determined by the equation 


( 2 . 2 . 1 ) 


{—A)^u = 0 in 17, 

' u = uq in IR"\17. 


To better explain this, let us fix a point x G 17. The expected value of the payoff at x 
is the average of all the payoffs at the points x from which one can reach x, weighted 
by the probability of the jumps. That is, by writing x = x + khv, with v G dBi, 
k e N* and h > 0, as in the previous Chapter 2.1, we have that the probability of 


jump IS 


\dBi\ |A:|i+2' 


-. This leads to the formula 


u(x) = 




By changing v into —v we obtain 


ii(x) = 

and so, by summing up, 


\dB^ 


2 J. 

fceN* 


it(x + khv) 

it(x — khv) 

|l + 2s 


dn^-\v). 


dH^-\v) 


2u(x) = 


\dBi 


2 L 


BBi 1^1 

it(x + khv) + u{x — khv) 


\k\ 


l + 2s 


d'H'^-^v). 


Since the total probability is 1, we can subtract 2ii(x) to both sides and obtain that 


0 = 




i(x -I- khv) + u{x — khv) — 2u(x) 


\k\ 


l+2s 


dn^-\v). 


We can now divide by and recognize a Riemann sum, which, after passing to 

the limit as /i \ 0, gives 0 = —(—A)®i((x), that is (2.2.1). 












CHAPTER 3 


An introduction to the fractional Laplacian 


We introduce here some preliminary notions on the fractional Laplacian and 
on fractional Sobolev spaces. Moreover, we present an explicit example of an s- 
harmonic function on the positive half-line M+, an example of a function with 
constant Laplacian on the ball, discuss some maximum principles and a Harnack 
inequality, and present a quite surprising local density property of s-harmonic func¬ 
tions into the space of smooth functions. 


3.1. Preliminary notions 


We introduce here the fractional Laplace operator, the fractional Sobolev spaces 
and give some useful pieces of notation. We also refer to [58] for further details 
related to the topic. 

We consider the Schwartz space of rapidly decaying functions defined as 


) := J/G 


Vo, (3 e Nq, sup \x°'di 3 f{x)\ < oo 


For any / e 5(IR"'), denoting the space variable x e and the frequency vari¬ 
able ^ e K", the Fourier transform and the inverse Fourier transform are defined, 
respectively, as 

(3.1.1) /(^) := ^/(^) := r /(x)e-2— 
and 

(3.1.2) fix) = 

jRri 

Another useful notion is the one of principal value, namely we consider the definition 


(3.1.3) 



ujx) - ujy) 

\x-y\n+2s y 


lim 


f 


ujx) - ujy) 

\x-y\r^+2s 


Notice indeed that the integrand above is singular when y is in a neighborhood 
of x, and this singularity is, in general, not integrable (in the sense of Lebesgue): 
indeed notice that, near x, we have that uix) — u(y) behaves at the first order 
like Vu(x) • {x — y), hence the integral above behaves at the first order like 


Vu(x) • (x — 
\x — 


whose absolute value gives an infinite integral near x (unless either Vu(x) = 0 
or s < 1/2). 

The idea of the definition in (3.1.3) is that the term in (3.1.4) averages out in a 
neighborhood of x by symmetry, since the term is odd with respect to x, and so it 
does not contribute to the integral if we perform it in a symmetric way. In a sense. 
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the principal value in (3.1.3) kills the first order of the function at the numerator, 
which produces a linear growth, and focuses on the second order remainders. 


The notation in (3.1.3) allows us to write (2.0.1) in the following more compact 
form: 


(-A)"u(a;) = 


C{n, s) r 2u(x) — u{x + y) — u{x — y) 


i 


dy 


2 

2 e->0 

2 £—>0 
Cin, s) lim 


I 

I. 

I. 


|y|n+2s 

2u(x) — u(x + y) — u(x — y) 


u(x) — u(x + y) 

."\Be 


|y|n+2s 

u{x) — u{r]) 
u{x) — u{r]) 


dy + 


I 


dy 

u{x) — u{x — y) 


drj + 


J 


AS, 

u(x)-u(C) 
ABeW k-CI" + '^ 


dy 


dC 


dr], 


\x - ry|"+2s 

where the changes of variable rj \= x + y and (^ := x — y were used, i.e. 


(3.1.5) (-A)"u(t) = (17(71, s)P.V. r \ 

I 


u{x) - u{y) 


■y\ 


n+2s 


-dy. 


The simplification above also explains why it was convenient to write (2.0.1) with 
the factor 2 dividing C{n, s). Notice that the expression in (2.0.1) does not require 
the P.V. formulation since, for instance, taking u e L®(IR”) and locally (7^, using a 
Taylor expansion of it in Si, one observes that 



u{x + y)- u{x -y)\ , 

|y|n+2s 


< I|m||l‘=°(R") 
^ I|w||l«>(r") 



dy, 


and the integrals above provide a finite quantity. 


Formula (3.1.5) has also a stimulating analogy with the classical Laplacian. 
Namely, the classical Laplacian (up to normalizing constants) is the measure of the 
infinitesimal displacement of a function in average (this is the “elastic” property of 
harmonic functions, whose value at a given point tends to revert to the average in 
a ball). Indeed, by canceling the odd contributions, and using that 

\x-y\^dy=y\\ {xk-ykfdy = n\ {x^-y^)'^dy, 

f^_-^JBj.{x) JBr{x) 

for any i e {1,..., n}. 
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we see that 


(3.1.6) 


lim u{x) 


r—tO r 


\Br{x)\ 


u{y) 

JBr-ix) 


dy 


— 2IR I Wj/) “ “(*)) 

r^o r^\Br[x)\ 

If 1 

— n +2 IR I \/u{x) ■ (x-y) +-D^u{x){x-y) ■ (x-y) 

r^o r"+^|Si| 2 


r(x) 

0{\x-y\^)dy 


1 " f 


= lim — - 

r—>-0 2r"'+^ 




^iMx) {Xi - yif dy 


= lim — ^ 

2ti 7 *'^+^ 

= - CnAu{x), 


1 "■ r 

\Bi\ Jbux) 


for some C„ > 0. In this spirit, when we compare the latter formula with (3.1.5), 
we can think that the fractional Laplacian corresponds to a weighted average of the 
function’s oscillation. While the average in (3.1.6) is localized in the vicinity of a 
point X, the one in (3.1.5) occurs in the whole space (though it decays at infinity). 
Also, the spacial homogeneity of the average in (3.1.6) has an extra factor that is 
proportional to the space variables to the power —2, while the corresponding power 
in the average in (3.1.5) is —2s (and this is consistent for s —> 1). 


Furthermore, for u e 5(K") the fractional Laplace operator can be expressed 
in Fourier frequency variables multiplied by (27r|^|)^®, as stated in the following 
lemma. 

Lemma 3.1.1. We have that 

(3.1.7) (-A)M^)=-^-'((2T|e|)^*S(a). 

Roughly speaking, formula (3.1.7) characterizes the fractional Laplace operator 
in the Fourier space, by taking the s-power of the multiplier associated to the 
classical Laplacian operator. Indeed, by using the inverse Fourier transform, one 
has that 

-Au{x) = -A{T-\u)){x) = -A r u(C)e2’^“f d£, 

J]R^ 

= f (27r|e|)22(^)e2-“«d^ = 7^-i((2T|e|)^2(C)), 

which gives that the classical Laplacian acts in a Fourier space as a multiplier 
of (27r|^|)^. From this and Lemma 3.1.1 it also follows that the classical Laplacian 
is the limit case of the fractional one, namely for any u e 5(M") 

lim(—A)®u = —Au and also lim(—A )'^m = u. 

s—>-l s—»-0 

Let us now prove that indeed the two formulations (2.0.1) and (3.1.7) are equiv¬ 
alent. 
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Proof of Lemma 3.1.1. Consider identity (2.0.1) and apply the Fourier trans¬ 
form to obtain 




C(n, s) f (2u{x) - u{x + y)- u{x - y)j 


(3.1.8) 


2 

C{n, s) 


I 

f uiO 

Je” 

C{n,s)u{^) r 
Je^ 

.riablt 

m--= f 


\y\ 


n+2s 


dy 


2 _ g27ri{ i/ _ g-27ri{ y 


|y|n+2. 

1 — cos(27r^ • y) 


dy 


\y\ 


n+2s 


dy. 


Now, we use the change of variable z = |^|y and obtain that 

1 — cos(27r^ • y) 


I 2 /I” 


dy 


JR'. 


1 — COS ^ • z 


|n+2s 


■ dz. 


JR- \Zl 

Now we use that J is rotationally invariant. More precisely, we consider a rotation 
R that sends ei = (1,0, ...,0) into ^/|^|, that is Rei = ^/|^|, and we call its 
transpose. Then, by using the change of variables ui = R?"z we have that 

2 s f 1 — cos(27ri?ei • z) 


m = icp* f 

JR'' 


■\n-\-2s 


dz 


2 s f 1 — cos(27ri?’^z • ei) 


Jk 

JR- 


IR'^Z 


|n+2s 


dz 


1 — cos(27ra;i) 


- duj. 


\^\n+2s 

Changing variables oj = 2 ttui (we still write w as a variable of integration), we obtain 
that 

cos Wi 


(3.1.9) 


J(C) = (2^|^|)^^ f 

Je"^ I 


,|n+2s 


■ duj. 


Notice that this latter integral is finite. Indeed, integrating outside the ball Bi we 
have that 


r 

Jr"\Si I 


cos (jJi \ 


duj ^ 


J 


|^|n+2s 


< 00 , 


|w|”+2« 

while inside the ball we can use the Taylor expansion of the cosine function and 
observe that 


I, 


1 — cos oji I 

|^|n+2s 

Therefore, by taking 


dio ^ 


I 


duj ^ 


(3.1.10) C{n,s) := ( I 

\ Jr’ 

it follows from (3.1.9) that 

m = 


l^\n+2s 

1 — COSWI 


I 


dco 


|^|n+2s 




dw 


|n+2s—2 


< 00 . 


w 


(SttICI) 


2s 


C{n, s) 

By inserting this into (3.1.8), we obtain that 


J- 


((-A)®u(x)) = C{n,s)u{^)J{C) = (27r|^|)2"zi(^), 
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which concludes the proof. □ 

Notice that the renormalization constant C{n,s) introduced in (2.0.1) is now 
computed in (3.1.10). 


Another approach to the fractional Laplacian comes from the theory of semi¬ 
groups (or, equivalently from the fractional calculus arising in subordination iden¬ 
tities). This technique is classical (see [143]), but it has also been efficiently used 
in recent research papers (see for instance [52, 134, 34[). Roughly speaking, the 
main idea underneath the semigroup approach comes from the following explicit 
formulas for the Euler’s function: for any A > 0, one uses an integration by parts 
and the substitution r = At to see that 
-sr(-s) = r(i - s) 


r+oo 

i 

-f 


T ^ dr 


d 

dr 


J '+oo 

■ 

0 


_ —S—1/ —T 

r (e 


1) dr 


= — s\ 


^+00 

'I ' 


J. — S — 1 ( ^ — 




1) dt, 


that is 

(3.1.11) = 

r(-s) Jo 

Then one applies formally this identity to A := —A. Of course, this formal step 
needs to be justified, but if things go well one obtains 

r(-s) Jo 

that is (interpreting 1 as the identity operator) 


1 c+°o 

(3.1.12) (—A)®u(a;) = ——r — u{x)) dt. 

r(-s) Jo 

Formally, if U{x,t) := e*^u{x), we have that U{x,0) = u(x) and 

dfU = -^(e‘^M(x)) = Ae*^u{x) = At/, 
dt 

that is U{x,t) = e*^u{x) can be interpreted as the solution of the heat equation 
with initial datum u. We indeed point out that these formal computations can be 
justified: 


Lemma 3.1.2. Formula (3.1.12) holds true. That is, if u e iS(]R") and U = 
U{x,t) is the solution of the heat equation 


then 

(3.1.13) 


dtU = AU 

^\t=o = 


m t > 0, 


{-Ayu{x) 


1 r+” . 

r^Jo ^ 


{U{x, t) — u{x)) dt. 
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Proof. From Theorem 1 on page 47 in [70] we know that U is obtained by 
Gaussian convolution with unit mass, i.e. 


(3.1.14) 


U{x,t)= G{x-y,t)u{y)dy = G{y,t) u{x - y) dy, 


where G{x,t) ■= (47rt)-"/2e-PI /G*). 

As a consequence, using the substitution r := |yp/(4t), 

r+00 

t) — u{x)) dt 

Jo 

^ + 00 ~ n 

= t~^~^G{y,t){uix-y)-u{x))dy dt 

Jo LJm^ 

^ + 00 ” r 

= (u(x — y) — u{x)) dy 

Jo LJr^ 

^ + 00 ” r 

T”/^(7r|j/p)“”/^|y|"2®(4r)®+^e“^ (^(a; - 1 /) - M(a;)) dy 

Jo LJr^ 

p + X) r r 

\-n/2 

Jo LJm 


Jo 

= 


"+S-1 + y) + u(x -y)- 


\y\ 


n+2s 


dy 


dr 

47-2 

dr. 


Now we notice that 


Jo ■ 


r+s— 


e '^(ir = r — + s), 


so we obtain that 


p+oo 

t~^~^{U{x,t) — u{x)) dt 
Jo 


_ 22s-l.jj.-'n./2p 


+ s) J 

^ Jr^ 


i(x + y) + u(x — y) — 2u(x) 


|y| 


n-|-2s 


dy 


22s ^-n/2p ^ 


(-A)^u(x). 


G{n, s) 

This proves (3.1.13), by choosing C{n,s) appropriately. And, as a matter of fact, 
gives the explicit value of the constant C{n, s) as 

22«r(f+ s) _ 22'*sr(f+ s) 


(3.1.15) 


C(n, s) = —- 


7r"/2r(—s) 7r"/2r(l — s) ’ 
where we have used again that r(l — s) = —sT{—s), for any s e (0,1). 


□ 


It is worth pointing out that the renormalization constant G(n, s) introduced 
in (2.0.1) has now been explicitly computed in (3.1.15). Notice that the choices 
of C{n,s) in (3.1.10) and (3.1.15) must agree (since we have computed the frac¬ 
tional Laplacian in two different ways): for a computation that shows that the 
quantity in (3.1.10) coincides with the one in (3.1.15), see Theorem 3.9 in [23]. For 
completeness, we give below a direct proof that the settings in (3.1.10) and (3.1.15) 
are the same, by using Fourier methods and (3.1.11). 


Lemma 3.1.3. For any n e N, n > 1, and s e (0,1), we have that 


(3.1.16) 


J 

Jr' 


1 — cos(27ra;i) 


,|n-|-2s 


du) = 


7r?+2'’r(P 


s) 


sr(f+ s) 
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Equivalently, we have that 


(3.1.17) 


J - 


— cosa;i , ttS r(l — 5 ) 
duj = 


,|n+2s 


22«sr(f + s)' 


Proof. Of course, formula (3.1.16) is equivalent to (3.1.17) (after the substi¬ 
tution Cj := 2710;). Strictly speaking, in Lemma 3.1.1 (compare (2.0.1), (3.1.7), 
and (3.1.10)) we have proved that 
(3.1.18) 


1 


JR' 


1 — COSCJi 


.|n+2s 


du! 


J 

JR- 


2u{x) — u{x y) — u{x — y) 


\y\ 


n+2s 


dy = E ^((27r|C|)2®u(^)). 


Similarly, by means of Lemma 3.1.2 (compare (2.0.1), (3.1.13) and (3.1.15)) we 
know that 


(3.1.19) 


2^® ^sr(^-l-s) C 2u{x) — u{x + y) — u{x — y) 
7r"/2r(l — s) 


I. 


dy 


|y|n+2s 

1 C+® 

——- t-^-\U{x,t)-u{x))dt. 

r(-s) Jo 

Moreover (see (3.1.14)), we have that U{x,t) := Ft * u{x), where 
rt(x) := G{x,t) = (47rt)-”/2e-PI'/(4‘). 

We recall that the Fourier transform of a Gaussian is a Gaussian itself, namely 

j-(e--l"l') = 

therefore, for any fixed t > 0, using the substitution y := xjJiTTt, 


-^r*(a = 


(47rt) 


t)"/^ L 




\ 

JR'‘ 


- I y I ^ c - 2 Ti-i y ■ () 


e e 

2 + 1 c |2 


dy 


= g-471- t|C| 


As a consequence 


F{U{x, t) — u(x)) = E(rt * u(x) — u(x)) 

= E(r ,. u)(0 - S(0 = i^m) - i)u(0 


- 1 ) 2 ( 0 - 


We multiply by t ^ ^ and integrate over t > 0, and we obtain 




G + OO ^ + 00 

t~’^~^{U{x,t) — u{x)) dt = ^-s-l(g-47r t|5l —l)dtu{^) 

Jo Jo 


= r(-s) (4^^100^2(0, 


thanks to (3.1.11) (used here with A := Itt^IO^). By taking the inverse Fourier 
transform, we have 

^ + 00 

t-^-^{U{x,t) - u{x)) dt = F(-s) (27r)2"J-0l0^'*2(O). 

Jo 
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We insert this information into (3.1.19) and we get 
22 s-i gp r 2u{x) — u{x + v) — u{x — y) 


\ 

JR- 


\y\ 


n+2s 


7r"/2r(l _ g) 

Hence, recalling (3.1.18), 

22 s-i g p (ra ^ 2u{x) — u{x + y) — u{x — y) 


dy=(27r)2^.F-i(|^P^2(a). 


7 r"/2r(l — s) 


J 

JR' 


dy 


1 


1 — COSWI 

, I, ,|n+2s 
jR" 1^1 

which gives the desired result. 


•■j 

JR 


did 


J 


\y\n+2s 

2u{x) — u{x + y) — u{x — y) 


|y| 


n+2s 


dy, 


□ 


For the sake of completeness, a different proof of Lemma 3.1.3 will be given 
in Appendix A.2. There, to prove Lemma 3.1.3, we will use the theory of special 
functions rather than the fractional Laplacian. For other approaches to the proof of 
Lemma 3.1.3 see also the recent PhD dissertations [76] (and related [77]) and [93] . 


3.2. Fractional Sobolev Inequality and Generalized Coarea Formula 

Fractional Sobolev spaces enjoy quite a number of important functional inequal¬ 
ities. It is almost impossible to list here all the results and the possible applications, 
therefore we will only present two important inequalities which have a simple and 
nice proof, namely the fractional Sobolev Inequality and the Generalized Coarea 
Formula. 

The fractional Sobolev Inequality can be written as follows: 


Theorem 3.2.1. For any s e (0,1), p e (l, and u e C'®(M"’), 


(3.2.1) 



|M(a;) - u{y)Y‘ 

\x - y|"+*P 



1 

P 


for some C > 0, depending only on n and p. 


Proof. We follow here the very nice proof given in [117] (where, in turn, the 
proof is attributed to Hai'm Brezis). We fix r > 0, a > 0, P > 0 and x e K". Then, 
for any y e M", 

|M(a:)| < \u{x) - u{y)\ + \u{y)\, 
and so, integrating over Bj.{x), we obtain 


l^r-l |w(a:)| 




}Br{x) 


\u{x) - u{y)\dy + \ 

J J 


{y)\dy 


)BA^) \x-y\‘^ 

|u(a:) - u{y)\ 


S r“[ l‘‘W-"fa)L „+ [ 

Jba^) f “ yr Jbaoo) 


BAx) 
u{y)\dy. 


|u(j/)| dy 
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Now we choose a := and we make use of the Holder Inequality (with expo¬ 
nents p and and with exponents and to obtain 


|i?r||u(a;)| ^ r p 


I 


|u(a:) - u{y)\ 


^ T P 


Br-{x) \x- y\ p 

n±£P f f lu(x) - u(y)lP 


VL.(X) \x-y\^+^P 

+ \u{y)\^ dy 


dy + 

dy 


f \u{y)\ 

Jb^(x) 

1 


1 


dy 


Br-(x) 


dy 


a(p-l) + sp 


dy 


Br{x) 


^ Cr 


\u{x) - u{y)\P 


"dy 


(L.(.) la; - 2/1"+"^ 

+Cr |u(i/)|^ 

\JbAx) J 


for some C > 0. So, we divide by r" and we possibly rename C. In this way, we 
obtain 

|M(a;) - u{y)Y’ 


l{x)\ ^ Cr" 


(I. 


JBAx) \x-y\^+^p 
That is, using the short notation 


dy \ + r 


f Hy)\' 

JB^{x) 


.-ap (j^y 


a:= J 

JlR 

and /? := 


\u{x) - u{y)\P 
\x - y|”+®P 

\u{y)\^ dy 

jRri 


dy 


we have that 


i(a;) I ^ Cr" p -I- r p (3 pp ^ , 


hence, raising both terms at the appropriate power and renaming C 


(3.2.2) 

We take now 


np nsp / 1 ri n — sp \ gp 

|ii(x)| '^-^p ^ Cv^-^p + r~p /3 


ar 

n — sp 


With this setting, we have that r p j3 pp is equal to a p. Accordingly, possibly 
renaming C, we infer from (3.2.2) that 


ttur "““p 


^ Cap" 


= C 


j 


|u(a;) - u{y)\P 


dy 


\x - J/|"+®P 

for some C > 0, and so, integrating over x e M", 

Wjx) - u{y)\P 
\x - y\^+^P' 


f Hy)\-' 

Jh" 


.-ap dy 


r im(x)|"-pp dx ^ c [ r r 

JRp \Jr" Jr' 


dx dy 


f Hy)\-^ 

jRri 


■-sp dy 
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This, after a simplification, gives (3.2.1). 


□ 


What follows is the Generalized Co-area Formula of [139] (the link with the 
classical Co-area Formula will be indeed more evident in terms of the fractional 
perimeter functional discussed in Chapter 6) . 


Theorem 3.2.2. For any s e (0,1) and any measurable function u : 12 — > [0,1], 

1 f f f‘(f f 

2 Jn Jo [a: — ?/|"+® Jq u(x)>t} Jfi/en, u(2/)sst} ~ / 

Proof. We claim that for any x, y B Vt 
c ^ 

(3.2.3) |u(x) - u{y)\ = J (x{u>t}(x) X{u^t}iv) + X{u^t}(x) X{u>t} (y)) dt. 

To prove this, we fix x and y in 17, and by possibly exchanging them, we can suppose 
that u{x) > u{y). Then, we define 

Fit) ■= X{u>t}ix) X{u^t}iy) + X{«s:t}(a;) X{«>t}(y)- 


By construction 


therefore 


Fit) = 


0 if t < uiy) and t > u(a;), 
1 if uiy) < uix), 


rl ru{x) 

(pit) dt = dt = uix) — u{y), 

Jo a«(y) 


which proves (3.2.3). 

So, multiplying by the singular kernel and integrating (3.2.3) over 17 x 17, we 
obtain that 

\uix) - uiy)\ 


X — 


■ dx dy 


JJ' 

rY r r X{u>t}ix) X{u^t}iy) + X{u^t}ix) X{u>t}iy) 

Jo \JnJn |a;-?/|"+® J 

L (L„L,„t^+LJ, 


= 2 


0 yj{u>t} 1^ y\ 

dx dy 


J J J 


\x - y| 


n+s 


J {u^t} J{u>t} 
dt, 


dx dy 
\x - 


as desired. 


□ 


3.3. Maximum Principle and Harnack Inequality 

The Harnack Inequality and the Maximum Principle for harmonic functions are 
classical topics in elliptic regularity theory. Namely, in the classical case, if a non¬ 
negative function is harmonic in Bi and r e (0,1), then its minimum and maximum 
in Br must always be comparable (in particular, the function cannot touch the level 
zero in B^). 

It is worth pointing out that the fractional counterpart of these facts is, in 
general, false, as this next simple result shows (see [95]): 
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Theorem 3.3.1. There exists a bounded function u which is s-harmonic in Bi, 
non-negative in Bi, but such that infu = 0 . 


Sketch of the proof. The main idea is that we are able to take the datum 
of u outside Bi in a suitable way as to “bend down” the function inside Bi until 
it reaches the level zero. Namely, let M > 0 and we take um to be the function 
satisfying 


(3.3.1) 


(—A)®um = 0 

^ Um = 1 ~ ^ 
UM = ^ 


in Bi, 
in B 3 \B 2 , 

in M"\((S 3 \.B 2 ) uBi). 


When M = 0, the function um is identically 1. When M > 0, we expect um 
to bend down, since the fact that the fractional Laplacian vanishes in Bi forces 
the second order quotient to vanish in average (recall ( 2 . 0 . 1 ), or the equivalent 
formulation in (3.1.5)). Indeed, we claim that there exists M* > 0 such that mm, ^ 0 
in Bi with inf = 0. Then, the result of Theorem 3.3.1 would be reached by 

Bi 

taking u := um^- 

To check the existence of such M*, we show that inf um ooasM^+oo. 

Bi 

Indeed, we argue by contradiction and suppose this cannot happen. Then, for 
any M > 0 , we would have that 


(3.3.2) 

for some fixed a e K. We set 


inf Um ^ —a, 

Bi 


VM ‘ = 


Um + M — 1 


Then, by (3.3.1), 

{-AYvm = 0 in Bi, 

VM = 0 in B3\B2, 

Um = 1 in K"\((B 3 \B 2 ) Cl Bi). 

Also, by (3.3.2), for any x e Bi, 

-a + M-1 


vm{x) ^ 


M 


By taking limits, one obtains that vm approaches a function Voo that satisfies 

(-A)^Uoo = 0 in Bi, 

Uoo = 0 in B 3 \B 2 , 

Uoo = 1 in K"\((B3\B2) c Bi) 

and, for any x e Bi, 

Vccix) > 1. 

In particular the maximum of Uqo is attained at some point a:* e Bi, with Uc»(a;*) 5= 
1. Accordingly, 


'vj 

JR'' 


0 = P. 

> p. 


Woo (a;*) - Vas{y) 


\x^_y\n+2s 
1-0 

iBiXB 


W. f — 

JB3\B2 F* 


y\ 


n+2s 


dy > P.V. ( 

Ji 

dy > 0 , 


WQo(a:«) - Vosiy) 
b3\b2 k* - 


dy 
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which is a contradiction. 


□ 


The example provided by Theorem 3.3.1 is not the end of the story concerning 
the Harnack Inequality in the fractional setting. On the one hand, Theorem 3.3.1 is 
just a particular case of the very dramatic effect that the datum at infinity may have 
on the fractional Laplacian (a striking example of this phenomenon will be given 
in Section 3.5). On the other hand, the Harnack Inequality and the Maximum 
Principle hold true if, for instance, the sign of the function u is controlled in the 
whole of R”. 

We refer to [ ' ^ , 130, 95, 78] and to the references therein for a detailed 
introduction to the fractional Harnack Inequality, and to [56] for general estimates 
of this type. 

Just to point out the validity of a global Maximum Principle, we state in detail 
the following simple result: 


Theorem 3.3.2. If {—A)^u ^ 0 in Bi and 0 in W^\Bi, then u > 0 in Bi. 


Proof. Suppose by contradiction that the minimal point cc* e Bi satisfies u(a;*) < 
0. Then u{x^) is a minimum in R” (since u is positive outside Hi), if y e i ?2 we 
have that 2u(a;*) — + y) — u{xt — y) ^ 0. On the other hand, in R”\H 2 we 

have that a:* ± y e R"\Hi, hence u(a:* ± y) > 0. We thus have 


0 ^ 




< 


i 

JW 

J 

JW‘ 

J 

JR' 


2 u(x*) - M(a:* + y) - u(a;* - y) 


dy 


|y|n+2s 

2 u(x*) - u{x^ + y) - - y) 


'\B2 


'\B2 

This leads to a contradiction. 


2u(xt,} 


\y\ 

dy < 0. 


n+2s 


dy 


□ 


Similarly to Theorem 3.3.2, one can prove a Strong Maximum Principle, such 
as: 

Theorem 3.3.3. If (—A)'*u > 0 m Hi and 0 in R"\Hi, then u > 0 in Hi, 
unless u vanishes identieally. 

Proof. We observe that we already know that u > 0 in the whole of R”, 
thanks to Theorem 3.3.2. Hence, if u is not strictly positive, there exists Xq e Hi 
such that u{xo) = 0. This gives that 

2u{xo) - u{xo + y) - u{xo -y) , f u{xo + y) + u{xo - y) , 

--''"“■J..- 

Now both u{xo + y) and u{xo — y) are non-negative, hence the latter integral is less 
than or equal to zero, and so it must vanish identically, proving that u also vanishes 
identically. □ 

A simple version of a Harnack-type inequality in the fractional setting can be 
also obtained as follows: 

Proposition 3.3.4. Assume that (—A)'’n ^ 0 in H 2 , with u ^ 0 in the whole 
o/R". Then 

m(0 ) > c u(x) dx, 

JBi 



for a suitable c > 0. 
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Proof. Let P e with r(a;) e [0,1] for any x e M", and r(0) = 1. 

We fix e > 0, to be taken arbitrarily small at the end of this proof and set 

(3.3.3) T ] := u(0) + e > 0. 

We define ra(x) := 2?7r(a;) — a. Notice that if a > 2rj, then ra(a;) ^ 2r] — a < 0 ^ 
u{x) in the whole of M", hence the set {Pa < u in K”} is not empty, and we can 
define 

0 , 1 = := inf {Pa < u in K"}. 
aeR 

By construction 

(3.3.4) a* ^ 277. 

If a < ?7 then Pa(0) = 2rj — a > rj > u(0), therefore 

(3.3.5) a* > r ]. 

Notice that 


(3.3.6) Pq^ ^ u in the whole of R”. 

We claim that 

(3.3.7) there exists Xq g B 1/2 such that Pa,|,(a;o) = u(xo). 

To prove this, we suppose by contradiction that u > Pa,^ in Bij 2 , i-c. 

^ := min(M - Pq^) > 0. 

-Bi/2 


Also, if X G we have that 

u{x) — PoH=(a;) = u{x) — 2riT{x) + a* = u{x) + 0 ^= > a* > 77 , 

thanks to (3.3.5). As a consequence, for any x e K", 

u{x) - Pa,i= (a:) ^ min{/ 7 , 77 } =: > 0 . 

So, if we define ou := 0 ,,= — {^^/2), we have that aj < 0 ,,= and 

u{x) - Pa„(a;) = u{x) - Pa,, (x) - ^ > y > 0 - 

This is in contradiction with the definition of 0 = 1 = and so it proves (3.3.7). 

From (3.3.7) we have that xo G Bi/ 2 , hence (—A)® 7 i(xo) 3=0. Also |(—A)'’Pa„ (x)| 
277 1(—A)'’P(a;)| ^ Cr], for any x e R”, therefore, recalling (3.3.6) and (3.3.7), 

C 77 > (-Arr^Jxo) - (-Aru(xo) 

[ra,, (a;o) - Pa,, (a:o + v)] - [^(xo) - u{xo + y)] 


> 


C'(n,s)P.V. r 

Jr' 

C'(n,s)P.V. r 

Jr' 

C'(n,s)P.V. r 

JB 


|y|n+2s 

u{xo + y)- Pa*(a:o + v) , 
- 

ujxQ + 7/) - Pa„(xo + y) 

,|n+2s 


dy 


dy. 


tsic-xo) \y\'' 

Notice now that if y e Bi{—xo), then \y\ ^ |xo| + 1 < 2, thus we obtain 
C{n, s) 




2n+2s 


I 


Bi{-xo) 


[u(a:o + y)- PoH=(a;o + y)] dy. 
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Notice now that ra^(a:) = 2 ? 7 r(a;) — ^ 77 , due to (3.3.5), therefore we conclude 

that 

CV > ^}n+2s ( f u{xo + y) dy - r]\BM . 

^ \dBi{-xo) j 

That is, using the change of variable x := xq + y, recalling (3.3.3) and renaming the 
constants, we have 

C{u{Q) + e) = Crj ^ f u{x)dx, 

JBi 

hence the desired result follows by sending e ^ 0 . □ 


3.4. An s-harmonic function 


We provide here an explicit example of a function that is s-harmonic on the 
positive line K+ := (0,-l-oo). Namely, we prove the following result: 


Theorem 3.4.1. For any a: e K, let Ws{x) := x^ = max{x,0}®. Then 
{-Ayws{x) = 


-Cs\X\ 

0 


if X < 0, 
if X > 0, 


for a suitable eonstant Cg > 0. 



Figure 2. An s-harmonic function 


At a first glance, it may be quite surprising that the function x’f is s-harmonic 
in ( 0 ,-too), since such function is not smooth (but only continuous) uniformly up 
to the boundary, so let us try to give some heuristic explanations for it. 

We try to understand why the function x’f is s-harmonic in, say, the inter¬ 
val (0,1) when s e (0,1]. From the discussion in Section 2.2, we know that the 
s-harmonic function in ( 0 , 1 ) that agrees with outside ( 0 , 1 ) coincides with the 
expected value of a payoff, subject to a random walk (the random walk is classical 
when s = 1 and it presents jumps when s e (0,1)). If s = 1 and we start from the 
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middle of the interval, we have the same probability of being moved randomly to 
the left and to the right. This means that we have the same probability of exiting 
the interval (0,1) to the right (and so ending the process at x = 1, which gives 1 as 
payoff) or to the left (and so ending the process at x = 0, which gives 0 as payoff). 
Therefore the expected value starting at x = 1/2 is exactly the average between 0 
and 1, which is 1/2. Similarly, if we start the process at the point x = 1/4, we 
have the same probability of reaching the point 0 on the left and the point 1/2 to 
the right. Since we know that the payoff at x = 0 is 0 and the expected value of 
the payoff at x = 1/2 is 1/2, we deduce in this case that the expected value for the 
process starting at 1/4 is the average between 0 and 1/2, that is exactly 1/4. We 
can repeat this argument over and over, and obtain the (rather obvious) fact that 
the linear function is indeed harmonic in the classical sense. 

The argument above, which seems either trivial or unnecessarily complicated 
in the classical case, can be adapted when s e (0,1) and it can give a qualitative 
picture of the corresponding s-harmonic function. Let us start again the random 
walk, this time with jumps, at x = 1/2: in presence of jumps, we have the same 
probability of reaching the left interval (—oo, 0] and the right interval [1, +cx)). Now, 
the payoff at (—oo, 0] is 0, while the payoff at [1, +oo) is bigger than 1. This implies 
that the expected value at x = 1/2 is the average between 0 and something bigger 
than 1, which produces a value larger than 1/2. When repeating this argument over 
and over, we obtain a concavity property enjoyed by the s-harmonic functions in 
this case (the exact values prescribed in [1, -too) are not essential here, it would be 
enough that these values were monotone increasing and larger than 1). 



Figure 3. A payoff model: case s = 1 and s e (0,1) 


In a sense, therefore, this concavity properties and loss of Lipschitz regularity 
near the minimal boundary values is typical of nonlocal diffusion and it is due to the 
possibility of “reaching far away points”, which may increase the expected payoff. 


Now we present a complete, elementary proof of Theorem 3.4.1. The proof 
originated from some pleasant discussions with Fernando Soria and it is based on 
some rather surprising integral cancellations. The reader who wishes to skip this 
proof can go directly to Section 3.3 on page 20. Moreover, a shorter, but technically 
more advanced proof, is presented in Appendix A.l. 

Here, we start with some preliminary computations. 


Lemma 3.4.2. For any s e (0,1) 

(1-f t)"-f (1 - t)" - 2 


J 

JO 


fl + 2s 


dt + 


I 


(1 -f ty 

^l + 2s 


dt = 
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Proof. Fixed e > 0, we integrate by parts: 
(l + t)® + (l-t)"-2 


I 


^l + 2s 


■ dt 


(3.4.1) 


1 


2s 

(l + e)" + (l-e)*-2 




dt 


- 2 " + 2 


+ 2 


1 r (1 + ty-^ - {1 -1) 


I 


\s-l 


l2s 


■ dt 


= ^ [o(i) - 2^ + 2] + i J (1 + ty-h-^^ dt-j (1 - ty-y-^^ dt^, 

with o(l) infinitesimal as c \ 0. Moreover, by changing variable t := t/{l — t), that 
is t := f/(l + i), we have that 


J 'l r + 00 

{1 - ty-y-^^ dt = {i + iy-y 

e Je/(1—e) 


' dt. 


Inserting this into (3.4.1) (and writing t instead of t as variable of integration), we 

obtain 

(3.4.2) 

(i + t)" + (i-t)"-2 


I 


^l+2s 


■ dt 


.^[„(1)_2.h.2] + 1 

= i[o(l)-2- + 2] + l 

Now we remark that 

re/{l-e) 


pi p + GO 

{1 + ty-y-"^^ dt - {1 + ty-y-'^^ dt 

Je/(1—e) 


pe/(l-e) p+GO 

J {1 + ty-y-^^ dt - j {i + ty-y-^^dt 


re/(i-e) ^e/(l-E) 

J (i + t)"-H-2*dt s: J {1 + ey-y-^^ dt = - e)-yi + sy-y 


therefore 


^e/(l-£) 

{i + ty-y-^^dt = o. 

Us 

So, by passing to the limit in (3.4.2), we get 


lim 
£\.0 . 


(3.4.3) r’ 

JO 


(1 + *)- + (1-T-2,,.-2- + 2 1 


^l + 2s 


2s 


r+cx) 


dt. 


Now, integrating by parts we see that 


2 ! + 


•dt = 


£+00 . 

Ji dt 

ns p+GO 

= -—+ t-^-^yi + ty dt. 

2s 


By plugging this into (3.4.3) we obtain that 


J’ 


{1 + ty + {i-ty -2 - 2 ®+ 2 2 


s r+00 


^1+2£ 

which gives the desired result. 


dt = 


2s 2s 


I 


+ —- t-^-^yi + tydt 


□ 
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From Lemma 3.4.2 we deduce the following (somehow unexpected) cancellation 
property: 

COROLARY 3.4.3. Let Ws be as in the statement of Theorem 3.4-. 1. Then 

= 0 . 

Proof. The function < (1 + t)® + (1 — t)® — 2 is even, therefore 


j: 


^ (i+t)®+(i-t)® 


|i| 


l+2s 


■ dt = 2 


I 


^ (l + f)® + (l-t)®-2 


^l+2s 


* dt. 


Moreover, by changing variable t := —t, we have that 


£ 


(l-t)®-2 


l£ 


l+2s 


dt = 


r+oo 

J1 


(1 + iy - 2 

ll + 2s 


dt. 


Therefore 


?ils(l + t) + Ws(l -t) - 2Ws(l) 


— 00 
-1 


|i| 


l+2s 


dt 


{i-ty-2^^ (i + f)® + (i-t)®-2 


l+2s 


■ dt + 


i; 


= 2 


I 


^ (l + t)® + (l-t)®-2 


^l + 2s 


dt + 2 


\t\l + 2s 

+” (1 + t)® - 2 


= 2 


= 2 


J’ 


{i-^ty + {i-ty-2 


i 


ll+2s 

+” dt 

^l+2s 


dt + 


^l + 2s 


dt + 

dt 


r+OD 

Jl 


(1 + 0®-2 


|l + 2s 


dt 


^-^dt-2r^ 

i.l + 2s i.l + 2s 


where Lemma 3.4.2 was used in the last line. Since 

dt 1 


£ 


we obtain that 


I 


^l + 2s 2s’ 

4(;s(l + t) + Ws{l -t) - 2wsil) 


|i| 


l + 2s 


dt = 0, 


that proves the desired claim. 


□ 


The counterpart of Corollary 3.4.3 is given by the following simple observation: 
Lemma 3.4.4. Let Wg be as in the statement of Theorem 3.4. L Then 

-i-AYwsi-l) > 0. 

Proof. We have that 

Ws(-1 + 0 + Ws(-1 -t) - 2ws{-l) = (-1 + t)+ + (-1 -0+^0 
and not identically zero, which implies the desired result. □ 

We have now all the elements to proceed to the proof of Theorem 3.4.1. 
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Proof of Theorem 3.4.1. We let a e {+1,-1} denote the sign of a fixed 
X e IR\{0}. We claim that 

Ws(ct(1 + i)) + Ws(cr(l - t)) - 2ws{a) 


(3.4.4) 


— 00 
+ 00 


|l + 2s 


dt 


,(a + t) + Ws(a — t) — 2ws{<j) 
|l + 2s ' 


dt. 


I-OO W 

Indeed, the formula above is obvious when a; > 0 (i.e. a = 1), so we suppose x < 0 
(i.e. cr = — 1) and we change variable r := —t, to see that, in this case. 


Ws(ct( 1 + t)) + Ws(cr(l - t)) - 2ws{cr) 


— 00 
+ 00 


|l + 2s 


dt 


— 00 
+ 00 


— 00 
+ 00 


g(—1 Ws{ — 1 + t) — 2Ws{(T) 

Ws{-1 + r) + iCs(-l -t) - 2ws{<7) 


|l + 2s 


dt 

dr 


Ws{a + r) + Ws{(T — t) — 2ws{a-) 


|l + 2s 


dr, 


thus checking (3.4.4). 

Now we observe that, for any r e M, 

w^sda^k) = (|a;|r)+ = |a:|V* = |x|®u;s(r). 

That is 

Ws{xr) = Ws{(j\x\r) = \x\^Ws{aT). 

So we change variable y = tx and we obtain that 


Ws{x + y) + Ws{x -y) - 2wsix) 


— 00 
+ 00 


\y\i+2s dy 

Ws{x{l + t)) + 'Ws{x{l — t)) — 2Ws{x) 


\xnt\ 


2s| + |l+2s 


dt 


= X 


Ws(cr(l + t)) + Ws(cr(l - t)) - 2ws(cr) 


— 00 
+ 00 


jtjWs 

Ws(a- + t) + Ws{(T — t) — 2ws{cr) 


dt 


\t\ 


l + 2s 


dt, 


where (3.4.4) was used in the last line. This says that 


{-Ayws{x) = 


' {—A)^Ws{—l) if X < 0, 


\x\ '*(—A)®Ws(l) if X > 0, 

hence the result in Theorem 3.4.1 follows from Corollary 3.4.3 and Lemma 3.4.4. □ 


3.5. All functions are locally s-harmonic up to a small error 

Here we will show that s-harmonic functions can locally approximate any given 
function, without any geometric constraints. This fact is rather surprising and it is 
a purely nonlocal feature, in the sense that it has no classical counterpart. Indeed, 
in the classical setting, harmonic functions are quite rigid, for instance they cannot 
have a strict local maximum, and therefore cannot approximate a function with 
a strict local maximum. The nonlocal picture is, conversely, completely different. 
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as the oscillation of a function “from far” can make the function locally harmonic, 
almost independently from its local behavior. 

We want to give here some hints on the proof of this approximation result: 


Theorem 3.5.1. Let k e N be fixed. Then for any f e C^{Bi) and any e > 0 
there exists R > 0 and u e n C'®(K"') sueh that 


(3.5.1) 

and 


{—A)^u{x) = 0 in Bi 

u = 0 in 

11/ ~ ^ 


Sketch of the proof. For the sake of convenience, we divide the proof into 
three steps. Also, for simplicity, we give the sketch of the proof in the one¬ 
dimensional case. See [63] for the entire and more general proof. 

Step 1. Reducing to monomials 

Let A: e N be fixed. We use first of all the Stone-Weierstrass Theorem and we have 
that for any £ > 0 and any / e 1]) there exists a polynomial P such that 


11/ ^llc'“(Bi) ^ 

Hence it is enough to prove Theorem 3.5.1 for polynomials. Then, by linearity, it 

N 


is enough to prove it for monomials. Indeed, if P{x) = ^ CmX'^ and one finds an 
s-harmonic function Um such that 


m=0 


\\Um - X 




< 


N 

then by taking m := ^ CmUm we have that 

m = l 


N 

I'M — P||(3.fc(g^) ^ 2 |Cm|||'am— X ||cfc(57)^e- 
m = l 

Notice that the function u is still s-harmonic, since the fractional Laplacian is a 
linear operator. 


Step 2. Spanning the derivatives 

We prove the existence of an s-harmonic function in i?i, vanishing outside a compact 
set and with arbitrarily large number of derivatives prescribed. That is, we show 
that for any m G N there exist i? > r > 0, a point x G K and a function u such that 

(—A)®u = 0 in (x — r, X -I- r), 

(3.5.2) V j /> 

u = 0 outside {x — R,x + R), 


and 

(3.5.3) 


D^u{x) = 0 for any j G {0,..., m — 1}, 
D'^u{x) = 1. 


To prove this, we argue by contradiction. 
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We consider Z to be the set of all pairs (m, x') of s-harmonic functions in a 
neighborhood of x, and points a; e K satisfying (3.5.2). To any pair, we associate 
the vector 

{u{x),Du{x ),..., D^u{x)) e 

and take V to be the vector space spanned by this construction, i.e. 

y := I [u{x),Du{x ),..., D'^u{x)), for (u, x) e z|. 

Notice indeed that 

(3.5.4) is a linear space. 

Indeed, let Vi, V 2 e V and oi, 02 e K. Then, for any i e {1, 2}, we have that T) = 
(ui{xi), Dui{xi),..., D"^Ui{xi)), for some {ui,Xi) G Z, i.e. Ui is s-harmonic in [xi — 
Ti, Xi -I- Xi) and vanishes outside {xi — Ri, Xi + Ri), for some Ri ^ Xi > 0. We set 

11 , 3 ( 0 ;) := aiUi(x + Xi) -I- a 2 U 2 {x + X 2 )- 

By construction, U 3 is s-harmonic in (—x^yXs), and it vanishes outside (—i? 3 ,i? 3 ), 
with X 3 := min{ri,r 2 } and R 3 := max{i?i, i? 2 }, therefore ( 103 ,0) G Z. Moreover 

D^u^{x) = aiD^ui{x + Xi) + a 2 D^U 2 {x + X 2 ) 

and thus 

fllVi + 02 V 2 

= ai(ui(a;i), £)ui(o:i),..., D'^ui{xi)) + a 2 {u 2 {x 2 ), Du 2 {x 2 )y ■ ■ .yD'^U 2 {x 2 )) 
This establishes (3.5.4). 

Now, to complete the proof of Step 2, it is enough to show that 

(3.5.5) F = ]R'"+^ 

Indeed, if (3.5.5) holds true, then in particular (0,..., 0,1) G V, which is the desired 
claim in Step 2. 

To prove (3.5.5), we argue by contradiction: if not, by (3.5.4), we have that V 
is a proper subspace of and so it lies in a hyperplane. 

Hence there exists a vector c = (cq, ..., Cm) e such that 

Vc s.t. c-C = 0}. 

That is, taking (uyx) e Z, the vector c = (cq, ... ,Cm) is orthogonal to any vector 
{u{x)yDu{x)y . . . , D'^u{x))y namely 

Yi CjD^u{x) = 0 . 

jsim 

To find a contradiction, we now choose an appropriate s-harmonic function u and 
we evaluate it at an appropriate point x. As a matter of fact, a good candidate for 
the s-harmonic function is a:^, as we know from Theorem 3.4.1: strictly speaking, 
this function is not allowed here, since it is not compactly supported, but let us 
say that one can construct a compactly supported s-harmonic function with the 
same behavior near the origin. With this slight caveat set aside, we compute for a 
(possibly small) x in ( 0 , 1 ): 

= s(s — 1 )... (s — j + 


D^x' 
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and multiplying the sum with (for a; 7 ^ 0 ) we have that 

Cjs{s — 1 ) .. .{s — j + = 0. 

j^m 

But since s e (0,1) the product s(s — 1)... (s — j + 1) never vanishes. Hence 
the polynomial is identically null if and only if cj = 0 for any j, and we reach 
a contradiction. This completes the proof of the existence of a function u that 
satisfies (3.5.2) and (3.5.3). 

Step 3. Rescaling argument and completion of the proof 

By Step 2, for any m 6 N we are able to construct a locally s-harmonic function u 
such that u{x) = a;'" + 0(a:™+^) near the origin (up to a translation). By considering 
the blow-up 

ux{x) = ^ =x"‘ + AO(x™+i) 

we have that for A small, ua is arbitrarily close to the monomial a;"*. As stated in 
Step 1, this concludes the proof of Theorem 3.5.1. □ 

It is worth pointing out that the flexibility of s-harmonic functions given by 
Theorem 3.5.1 may have concrete consequences. For instance, as a byproduct of 
Theorem 3.5.1, one has that a biological population with nonlocal dispersive atti¬ 
tudes can better locally adapt to a given distribution of resources (see e.g. Theo¬ 
rem 1.2 in [105]). Namely, nonlocal biological species may efficiently use distant 
resources and they can fit to the resources available nearby by consuming them 
(almost) completely, thus making more difficult for a different competing species to 
come into place. 


3.6. A function with constant fractional Laplacian on the ball 

We complete this chapter with the explicit computation of the fractional Lapla¬ 
cian of the function U(x) = (1 — |a;p) + . In Bi we have that 

{-^yU{x) = C{n,s)'^B{s,l- s), 

where C{n, s) is introduced in (3.1.10) and B is the special Beta function (see section 
6.2 in [4]). Just to give an idea of how such computation can be obtained, with 
small modifications respect to [69, 68[ we go through the proof of the result. The 
reader can And the more general result, i.e. for U{x) = (1 — for p > —1, in 

the above mentioned [69, 68[. 

Let us take u: K ^ M as u{x) = (1 — |a;p)^. Consider the regional fractional 
Laplacian restricted to (—1,1) 

N uix) — uiy) , 

Cu(x) ■'!> 

and compute its value at zero. By symmetry we have that 

/:u(o) = 2 lim £ ^ J ^ dy- 
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Changing the variable uj = and integrating by parts we get that 
Cu{Q) = 2 lim J dy - J (1 - dy^ 

Je^ 


1 .. 


=-h lim 


5 \ S 

2 s 28/1 ^ 2 \s 




=-h lim 

S £->-0 


g-2«_g-2«(i_^2 


+ J U! ®(1 —w)® ^dw^. 


Using the integral representation of the Gamma function (see [4], formula 6.2.1), 


i.e. 


B{a,P)= r 

Jo 

it yields that 

£u{0) = — s,s) -. 

s 

For X e Bi we use the change of variables oj = • We obtain that 


Cu{x) = 




-1 1 ^ 
2 \s 


= (l-x^)"*P.U. J 

(«.i) 


3._y|l + 2s 

^ {l-ujxf^-^-{l-ujy{l-ujx)-^ 


2 s + l 


duj 




2 s + l 


I 


| 2 s+l 


duo + 


dw 


i: 


(I-UJX)^’^-^ - 1 


| 2 s + l 


duj 




where we have recognized the regional fractional Laplacian and denoted 


/(x) P.V. 


In J{x) we have that 


J{x) = P. 


= lim 

er^O 


'KI. 

r 


^ (1 — ^ 




2 s+l 


duj — 




-l- 2 s 


d(jj 


(1 + OJxY^ ^ + (1 — U)x) 


2 s-l 


| 2 s + l 


‘duj — 2 


f 


-l- 2 s 
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With the change of variable t = — 

00 

J{x) = - + lirn ^ J |^(t + dt — - -^ 

( 3 . 6 . 2 ) 1 (1 + + (1 - 1 (l + ga :)^^ + ( l -£ a :)^^-2 

S 2s 2s e^o £2. 

_ 1 _ (1 + a;) 2 ^ + (1 - x) 2 « 
s 2s 

To compute I{x), with a Taylor expansion of (1 — oox)~^ at 0 we have that 


lix) = P.V. 

J-i 


|2s+l 


-(1 — oo^y duo. 


The odd part of the sum vanishes by symmetry, and so 


2 l.„ 

^ ’ e^O w2s + l ^ ’ 


= — 2 lim 
£^0 


GO pi 

k=l 


.2fc-2s-l 


(1 — uj^y doo. 


We change the variable t = and integrate by parts to obtain 


00 ^1 

i{x) = - lim y x^’^ - ty dt, 


= y lim 

^ £^0 


k = l 


yi-ey^ 


k — s 


k — s 


r y-yi-ty-Ut 


For k ^ 1, the limit for £ that goes to zero is null, and using the integral represen¬ 
tation of the Beta function, we have that 


d(x) = y 


-,2k 


k = l 


—s 
k — s 


B{k + 1 — s, s). 


We use the Pochhammer symbol defined as 

/o R ON / N _ J1 for A: = 0, 

^ ^ [ <z(<Z -f 1 ) • • • ((? -f fc — 1 ) for fc > 0 

and with some manipulations, we get 




k- 


{k-s)T{k+l) 

(-s)r(A: - s)r(s) 


k\ 

= B{1 — s, s) 


{-s)k 

kl 


I{x) = B{1 - s,s) y 


,2fc 


k=l 


i-s)k 
k\ ■ 


And so 
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By the definition of the hypergeometric function (see e.g. page 211 in [112]) we 
obtain 


00 2k 

I{x) = - B{1- s,s) + B{1- s,s)Yji.-s)k\r 

k=0 


= B{1 — s, s) 



1 1 

2 ’ 2 ’"^ 



Now, by 15.1.8 in [4] we have that 


F 




and therefore 


/(a:) = i3(l — s, s) ^(1 — — 1^ 

Inserting this and (3.6.2) into (3.6.1) we obtain 


(3.6.4) Cu{x) = B{1 — s, s) — (1 — x^) 

Now we write the fractional Laplacian of u as 


_(l + ^) 2 . + (l_^) 2 « 


(—A)®M(a;) 
C(l,s) 


Cu{x) + J 


u(x) 


[a; _ y|l + 2s 
I —1 


dy + 


f 


2 s 


u(x) 


yi 


l + 2s 


dy 


Cu(.i:) + (1 - x^yfj \x-y\ ^ + J \x-y\ ^ dy] 


= Cu(x) + (1 - + 

Inserting (3.6.4) into the computation, we obtain 
(3.6.5) (— A)®u(a;) = 0(1, s)B{l — s, s). 

To pass to the n-dimensional case, without loss of generality and up to rotations, 
we consider x = (0, 0,..., Xn) with x„ > 0. We change into polar coordinates 
X — y = th, with h e dBi and t ^ 0. We have that 


(3.6.6) 

(-A)-fi(x) 

C{n, s) 


= P.V. 


J 

J j 

JdBi \ Jl 


(i-ixi^)^-(i-i2/n 


2\s 


dy 


X — J/|"+2s 
{1 - \x\'^y - {1 - \x + ht\y 


\t\ 


l + 2s 


dt]dH'^-yh). 


Changing the variable t = —\x\hn + r-y/l/inXp — |xp + 1, we notice that 


1 — |x + Wp = (1 — T^)(l — + \hnx\'^) 
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and so 


p.v. r 

Jr 


(1 — \x\^Y — (1 — Icc + htf) 


2\s 


' dt 


=P.V. 


’■'■J 

Ji 


|i|l + 2s 

(1 - \x\‘^y - (1 - T2)®(|/i„a:p - |a;p + 1)" 


=P.V. 


Jr 


(-A)®u 


1 - 


\x\hn + T-yJlhnXl^ - \x\^ + 1 
\x\‘^hl 


l + 2s 


^/\hYx\^^^]xY^^ dr 


\hnx\^ - Ixp + 1 


-( 1 -T^)^ 


•\hr: 


- |a ;|2 + 1 

\x\K 


l+2s 


dr 


^y\hnx\'^ - |a:|2 + 1 


C(l,s) 

=B{1 - s,s), 

where the last equality follows from identity (3.6.5). Hence from (3.6.6) we conclude 
that 

(-A)«f^(x) = C(n,s)H(l-s,s)‘:J. 




















CHAPTER 4 


Extension problems 


We dedicate this part of the book to the harmonic extension of the fractional 
Laplacian. We present at first two applications, the water wave model and the 
Peierls-Nabarro model related to crystal dislocations, making clear how the exten¬ 
sion problem appears in these models. We conclude this part by discussing^ in detail 
the extension problem via the Fourier transform. 

The harmonic extension of the fractional Laplacian in the framework consid¬ 
ered here is due to Luis Caffarelli and Luis Silvestre (we refer to [31] for details). 
We also recall that this extension procedure was obtained by S. A. Molcanov and 
E. Ostrovskii in [109] by probabilistic methods (roughly speaking “embedding” a 
long jump random walk in M” into a classical random walk in one dimension more, 
see Figure 4). 



Figure 4. The random walk with jumps in K" can be seen as 
a classical random walk in 


The idea of this extension procedure is that the nonlocal operator (—A)® acting on 
functions defined on K" may be reduced to a local operator, acting on functions 
defined in the higher-dimensional half-space := M” x (0, -too). Indeed, take 


^Though we do not develop this approach here, it is worth mentioning that extended problems 
arise naturally also from the probabilistic interpretation described in Section 2. Roughly speaking, 
a stochastic process with jumps in R" can often be seen as the “trace” of a classical stochastic 
process in R" x [0,-too) (i.e., each time that the classical stochastic process in R" x [0,-too) 
hits R" X {0} it induces a jump process over R”). Similarly, stochastic process with jumps may 
also be seen as classical processes at discrete, random, time steps. 
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U: ^ K such that U{x, 0) = u{x) in M", solution to the equation 

div(2/i-2*V;7(a;,2/)) = 0 in 
Then up to constants one has that 

- lim (y^~'^^dyU{x,y)) = {-Ayu{x). 

4.1. Water wave model 

Let us consider the half space = M" x (0, +oo) endowed with the coor¬ 

dinates X e K" and y e (0,-too). We show that the half-Laplacian (namely when 
s = 1/2) arises when looking for a harmonic function in with given data 

on M" X {y = 0}. Thus, let us consider the following local Dirichlet-to-Neumann 
problem: 

Ja17 = 0 inK”+\ 

|?7(x, 0) = u(a:) forcceM”. 

The function U is the harmonic extension of u, we write U = Eu, and define the 
operator C as 

(4.1.1) Cu{x) ■.= —dyU{x,0). 

We claim that 

(4.1.2) C = 
in other words 

= -A.. 

Indeed, by using the fact that E{Cu) = —dyU (that can be proved, for instance, by 
using the Poisson kernel representation for the solution), we obtain that 

C^u{x) = C[Cu){x) 

= — dyE(^Cu){x,0) 

= - ^y{ - ^yU){x,0) 

= {dyyU + A,U - A,u){x,0) 

= AU{x,0) — Au{x) 

= — Au{x), 

which concludes the proof of (4.1.2). 

One remark in the above calculation lies in the choice of the sign of the square 
root of the operator. Namely, if we set Cu{x) := dyU{x,0), the same computation 
as above would give that C? = —A. In a sense, there is no surprise that a quadratic 
equation offers indeed two possible solutions. But a natural question is how to 
choose the “right” one. 

There are several reasons to justify the sign convention in (4.1.1). One reason is 
given by spectral theory, that makes the (fractional) Laplacian a negative definite 
operator. Let us discuss a purely geometric justification, in the simpler n = 1- 
dimensional case. We wonder how the solution of problem 


(4.1.3) 


(-A)®m =1 in (-1,1), 

u = 0 in IR\(—1,1) 



4.1. WATER WAVE MODEL 
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should look like in the extended variable y. First of all, by Maximum Principle 
(recall Theorems 3.3.2 and 3.3.3), we have that u is positive^ when x e (—1,1) 
(since this is an s-superharmonic function, with zero data outside). 

Then the harmonic extension [/ in y > 0 of a function u which is positive 
in (—1,1) and vanishes outside (—1,1) should have the shape of an elastic membrane 
over the halfplane that is constrained to the graph of u on the trace {y = 0}. 



We give a picture of this function U in Figure 5. Notice from the picture that 
dyU{x,0) is negative, for any x e (—1,1). Since (—A)®'u(ai) is positive, we de¬ 
duce that, to make our picture consistent with the maximum principle, we need to 
take the sign of C opposite to that of dyU{x,0). This gives a geometric justifica¬ 
tion of (4.1.1), which is only based on maximum principles (and on “how classical 
harmonic functions look like”). 

Application to the water waves. 

We show now that the operator C arises in the theory of water waves of irrota- 
tional, incompressible, inviscid fluids in the small amplitude, long wave regime. 

Consider a particle moving in the sea, which is, for us, the space K" x (0,1), 
where the bottom of the sea is set at level 1 and the surface at level 0 (see Figure 
6 ). The velocity of the particle is v: M" x (0,1) ^ and we write v{x,y) = 

(vj;{x,y),Vy{x,y)), where v^- K" x (0,1) M" is the horizontal component and 
Vy: K” x (0,1) ^ K is the vertical component. We are interested in the vertical 
velocity of the water at the surface of the sea. 


In our model, the water is incompressible, thus div u = 0 in K" x (0,1). Furthermore, 
on the bottom of sea (since water cannot penetrate into the sand), the velocity has 
only a non-null horizontal component, hence Vy{x,l) = 0. Also, in our model we 
assume that there are no vortices: at a mathematical level, this gives that v is 
irrotational, thus we may write it as the gradient of a function U: ^ M. 


^As a matter of fact, the solution of (4.1.3) is explicit and it is given by (1 — up to 

dimensional constants. See [69] for a list of functions whose fractional Laplacian can be explic¬ 
itly computed (unfortunately, differently from the classical cases, explicit computations in the 
fractional setting are available only for very few functions). 
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V 


y = 1 






Figure 6 . The water waves model 


This says that the vertical component of the velocity at the surface of the sea is 
Vy{x,0) = dyU{x,0). We are led to the problem 


'AU = 0 inK!;:+i 

^ dyU{x, 1) = 0 in K", 

U{x,0)=u{x) in K”. 


(4.1.4) 


Let £ be, as before, the operator Cu{x) := —dyU{x,0). We solve the problem 
(4.1.4) by using the Fourier transform and, up to a normalization factor, we obtain 
that 



Notice that for large frequencies this operator is asymptotic to the square root 
of the Laplacian: 



The operator £ in the two-dimensional case has an interesting property, that is 
in analogy to a conjecture of De Giorgi (the forthcoming Section 5.2 will give further 
details about it): more precisely, one considers entire, bounded, smooth, monotone 
solutions of the equation £u = f{u) for given /, and proves that the solution only 
depends on one variable. More precisely: 

Theorem 4.1.1. Let f e C'^(K) and u be a bounded smooth solution of 


£u = f{u) in 
dx^u >0 in 


Then there exist a direction uj e and a function uq • K —> K such that, for any 


X G 


u{x) = uo(x ■ w). 
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See Corollary 2 in [53] for a proof of Theorem 4.1.1 and to Theorem 1 in [53] 
for a more general result (in higher dimension). 


4.2. Crystal dislocation 


A crystal is a material whose atoms are displayed in a regular way. Due to some 
impurities in the material or to an external stress, some atoms may move from their 
rest positions. The system reacts to small modifications by pushing back towards 
the equilibrium. Nevertheless, slightly larger modifications may lead to plastic 
deformations. Indeed, if an atom dislocation is of the order of the periodicity size 
of the crystal, it can be perfectly compatible with the behavior of the material at a 
large scale, and it can lead to a permanent modification. 

Suitably superposed atom dislocations may also produce macroscopic deforma¬ 
tions of the material, and the atom dislocations may be moved by a suitable external 
force, which may be more effective if it happens to be compatible with the periodic 
structure of the crystal. 

These simple considerations may be framed into a mathematical setting, and 
they also have concrete applications in many industrial branches (for instance, in 
the production of a soda can, in order to change the shape of an aluminium sheet, 
it is reasonable to believe that applying the right force to it can be simpler and less 
expensive than melting the metal). 

It is also quite popular (see e.g. [103]) to describe the atom dislocation motion 
in crystals in analogy with the movement of caterpillar (roughly speaking, it is less 
expensive for the caterpillar to produce a defect in the alignment of its body and to 
dislocate this displacement, rather then rigidly translating his body on the ground). 


The mathematical framework of crystal dislocation presented here is related to 
the Peierls-Nabarro model, that is a hybrid model in which a discrete dislocation 
occurring along a slide line is incorporated in a continuum medium. The total 
energy in the Peierls-Nabarro model combines the elastic energy of the material in 
reaction to the single dislocations, and the potential energy of the misfit along the 
glide plane. The main result is that, at a macroscopic level, dislocations tend to 
concentrate at single points, following the natural periodicity of the crystal. 

To introduce a mathematical framework for crystal dislocation, first, we “slice” 
the crystal with a plane. The mathematical setting will be then, by symmetry 
arguments, the half-plane = {{x,y) e s.t. y ^ 0} and the glide line will be 
the cc-axis. In a crystalline structure, the atoms display periodically. Namely, the 
atoms on the x-axis have the preference of occupying integer sites. If atoms move 
out of their rest position due to a misfit, the material will have an elastic reaction, 
trying to restore the crystalline configuration. The tendency is to move back the 
atoms to their original positions, or to recreate, by translations, the natural periodic 
configuration. This effect may be modeled by defining v^{x) := u(a;, 0) to be the 
discrepancy between the position of the atom x and its rest position. Then, the 
misfit energy is 


(4.2.1) 
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Figure 7. Crystal dislocation 


where W is a smooth periodic potential, normalized in such a way that W{u+1) = 
W{u) for any it e K and 0 = hF(0) < W{u) for any u e (0,1). We also assume that 
the minimum of W is nondegenerate, i.e. W"(0) > 0. 

We consider the dislocation function v{x,y) on the half-plane The elastic 
energy of this model is given by 


(4.2.2) 



\/v{x,y) 


2 

dx dy. 


The total energy of the system is therefore 


(4.2.3) F{v) ■.= £{v) + M{v°) =- 


1 r \7v{x,y) dxdy+ f W(v{x,0)') dx. 

2 JbI Jk ^ '' 


Namely, the total energy of the system is the superposition of the energy in (4.2.1), 
which tends to settle all the atoms in their rest position (or in another position 
equivalent to it from the point of view of the periodic crystal), and the energy 
in (4.2.2), which is the elastic energy of the material itself. 


Notice that some approximations have been performed in this construction. 
For instance, the atom dislocation changes the structure of the crystal itself: to 
write (4.2.1), one is making the assumption that the dislocations of the single atoms 
do not destroy the periodicity of the crystal at a large scale, and it is indeed this 
“permanent” periodic structure that produces the potential W. 

Moreover, in (4.2.2), we are supposing that a “horizontal” atom displacement 
along the line {y = 0} causes a horizontal displacement at {y = e} as well. Of 
course, in real life, if an atom at {y = 0} moves, say, to the right, an atom at 
level {y = e} is dragged to the right as well, but also slightly downwards towards 
the slip line {y = 0}. Thus, in (4.2.2) we are neglecting this “vertical” displacement. 
This approximation is nevertheless reasonable since, on the one hand, one expects 
the vertical displacement to be negligible with respect to the horizontal one and, on 
the other hand, the vertical periodic structure of the crystal tends to avoid vertical 
displacements of the atoms outside the periodicity range (from the mathematical 
point of view, we notice that taking into account vertical displacements would make 
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the dislocation function vectorial, which would produce a system of equations, rather 
than one single equation for the system). 

Also, the initial assumption of slicing the crystal is based on some degree of 
simplification, since this comes to studying dislocation curves in spaces which are 
“transversal” to the slice plane. 

In any case, we will take these (reasonable, after all) simplifying assumptions 
for granted, we will study their mathematical consequences and see how the results 
obtained agree with the physical experience. 


To find the Euler-Lagrange equation associated to (4.2.3), let us consider a 
perturbation (j) e (^“(IR^), with ip{x) := ^(a;, 0) and let u be a minimizer. Then 


—J'(u + e(j)) 


= 0 , 


£ = 0 


which gives 


f 

Jr- 


Vu • V(j)dxdy + 


r W'{v°)ip 

Jr 


dx = 0. 


Consider at first the case in which supp(/>n5IR^ = 0, thus v? = 0. By the Divergence 
Theorem we obtain that 

r (j) Av dxdy = 0 for any (j) e (7® (IR^), 

Jr2 

If supp0 n dlR^ A 0 then we have that 

0=1 div(^V v) dx dy + j W'(v^)(pdx 

Jr2 Jr 

r W'(v°)y?< 

Jr 


thus Av = 0 in 


J 

Ja 


,dv 
h— dx + 


r (p^ dx + r w'{v^)<p 

Jr Jr 


dx 

dx 


for an arbitrary ip e C'^(IR) therefore — (a:,0) = W'{v^{x)) for a; e IR. Hence the 

dy 

critical points of T are solutions of the problem 


Av{x, y) = 0 for a: e ]R and y > 0, 

< v{x, 0) = v°{x) for a: e ]R, 

dyu(x,0) = IT'^u(a;,0)^ for a: e R 

and up to a normalization constant, recalling (4.1.1) and (4.1.2), we have that 

—V—Au(a;,0) = W' (y{x,Q)), for any a: G R. 

The corresponding parabolic evolution equation is dtv{x,Q) = —^—Av (x,0) - 
W'(v(x,0)). 


After this discussion, one is lead to consider the more general case of the frac¬ 
tional Laplacian (—A)® for any s e (0,1) (not only the half Laplacian), and the 
corresponding parabolic equation 

dfV = —(—A)®u — W'(y) + O', 
where ct is a (small) external stress. 
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If we take the lattice of size e and rescale v and a as 

/ t x\ 2 s ( 

and a = e 

then the rescaled function satisfies 

(4.2.4) dtv^ = - ( — (—-■*" ™ (0, +oo) x M 

with the initial condition 


Ve{0,x) = v^{x) for a: e K. 

To suitably choose the initial condition , we introduce the basic layer^ solution u, 
that is, the unique solution of the problem 


|-(-A)®u(a;) = IT'(m) in M, 

> 0 and u{—co) = 0,u(0) = 1/2,^(+ 00 ) = 1. 


For the existence of such solution and its main properties see [114] and [26]. Fur¬ 
thermore, the solution decays polynomially at ±00 (see [61] and [60]), namely 


(4.2.6) 


i{x) — H{x) 


1 

2 slF"(0) |a/ 


X 

l+2s 


< 


c 


where 1 ? > 2s and H is the Heaviside step function 


for any x e K", 


H{x) 


1, X ^ 0 
0, X < 0. 


We take the initial condition of the solution of (4.2.4) to be the superposition of 
transitions all occurring with the same orientation, i.e. we set 

2s ^ / 0 \ 

(4.2.7) ^^(^’0) := :j^^a(0,x) + 2u(^^-^j, 

where x°,..., x^ are N fixed points. 


The main result in this setting is that the solution approaches, as e —> 0, the 
superposition of step functions. The discontinuities of the limit function occur at 
some points {xi{t))^_^ which move accordingly to the following^ dynamical 
system 


(4.2.8) 


Xj = 7 


a{t,Xi) + 2 




2s\Xi — Xj + i 


Xj(0) = x°, 


in (0, - 1 - 00 ), 


where 

(4.2.9) q=Q^(«')2^ 


More precisely, the main result obtained here is the following. 


3As a matter of fact, the solution of (4.2.5) coincides with the one of a one-dimensional 
fractional Allen-Cahn equation, that will be discussed in further detail in the forthcoming Section 
5.1. 

^ The system of ordinary differential equations in (4.2.8) has been extensively studied in [81]. 
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Xi 


Xn 


Xn 


X 


N 


Figure 8. The initial datum when £ —> 0 


Theorem 4.2.1. There exists a unique viscosity solution of 

dtVe = - ^ — {—Ayv^ -2s ^ ®"’ 

2 s ^ / 0 \ 

^e(0,x) = ^^„( 0 ) Cr( 0 ,a;) + ^ ' j V 6 


such that 
(4.2.10) 


N 


Vi^vyt,x) = Y^H[x- Xi{t)), 


2=1 


where {xi{t))^^^ ^ is solution to (4.2.8). 


We refer to [88] for the case s = 2 ’ [61] for the case ® ^ 2 ’ [60] for 

the case s < - (in these papers, it is also carefully stated in which sense the limit 
in (4.2.10) holds true). 

We would like to give now a formal (not rigorous) justification of the ODE 
system in (4.2.8) that drives the motion of the transition layers. 

Justification of ODE system (4.2.8). We assume for simplicity that the 
external stress a is null. We use the notation ~ to denote the equality up to 
negligible terms in e. Also, we denote 

ix-xyt)' 

Ui[t, X) := u' 

and, with a slight abuse of notation 


i'(t,x) := 


,fx- xyt) 
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By (4.2.6) we have that the layer solution is approximated by 


(4.2.11) 


Ui{t, x) ^ H 


X - Xi{t) 


’{x - Xi{t)) 


2sW"{Q)\x-Xi{t) 


l + 2s ■ 


We use the assumption that the solution is well approximated by the sum of N 
transitions and write 


N ^ _ (A 

Ve{t,x) ^ ^Ui{t,x) = p^ 


For that 


N 


dtVe{t,x) = - - Y Ui{t,x)Xi{t) 

2 = 1 

and, since the basic layer solution u is the solution of (4.2.5), we have that 

N 

-(-A)"u, ^ -Yj{-^yu^{t,x) 

2=1 

2=1 

1 ^ 

= 3lS»"l 


2=1 

N 


X - Xi{t) 

e 


2=1 

Now, returning to the parabolic equation (4.2.4) we have that 

1 AT 1 / ^ ^ \ 

(4.2.12) - - y|u'(t,x)i*(f) = XI W'{u^{t,x)) - W(^Y,u,{t,x)^ V 

^ 2 = 1 2=1 2 = 1 ^ 

Fix an integer k between 1 and N, multiply (4.2.12) by u'i.(t,x) and integrate over 
K. We obtain 


(4.2.13) 


1 ^ f 

u'yt,x)u'f.{t,x)dx 

^ Jk 


^2s + l 


N ^ 

§1 


r ^ 

W' x))u'j.{t, x) dx — j W' Ui{t, x)^ u'k^, x) dx 


We compute the left hand side of (4.2.13). First, we take the fc**' term of the sum 
(i.e. we consider the case i = k). By using the change of variables 


(4.2.14) 


y ■= 


Xk{t) 
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we have that 


(4.2.15) 



{u't,f{t,x)dx = 


- ^Xk{t) j dx 

-Xk{t) f {u'f{y)dy 

Jm 

Xk{t) 

5 

7 


where 7 is defined by (4.2.9). 

Then, we consider the term of the sum on the left hand side of (4.2.13). By 
performing again the substitution (4.2.14), we see that this term is 


-Xi{t) r u'^{t,x)u'i,{t,x) dx 

e Jb 



dx 


^ 0 , 


where, for the last equivalence we have used that for e small. 


u y + 


Xk{t) - Xi{t) 


is asymptotic to m'(±oo) = 0 . 

We consider the first member on the right hand side of the identity (4.2.13), 
and, as before, take the term of the sum. We do the substitution (4.2.14) and 
have that 


1 

e 


J 

Jr 


x) dx 


W'{u{y))u'{y)dy 

JR 

+00 

W{u{y)) 

— 00 

W(l) - W{0) = 0 


by the periodicity of W. Now we use (4.2.11), the periodicity of W' and we perform 
a Taylor expansion, noticing that W'{Q) = 0. We see that 


W {u^{t,x)) 



/ X - a:^(t) \ _ [x - x,{t)) \ 

V e ) 2sW"{0)\x-Xi{t)\^^‘^") 

e^^[x — Xi{t)) \ 

2sW"{0)\x - ) 


-e^^{x-x^{t)) 
2 s|a: — Xi(t)\^~^'^^ 


Therefore, the i**' term of the sum on the right hand side of the identity (4.2.13) for 
i ¥= k, by using the above approximation and doing one more time the substitution 
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(4.2.14), for e small becomes 
(4.2.16) 

- I W'[ui{t,x))uj^{t,x) dx 
e Jb 


1 r [x — Xi{t)) ,fx — Xk{t) 


I. 


2s\x — Xi{t) 


l + 2s ' 


dx 


j 

JR 


^2s 


{ey + Xkit) - Xi{t)) , 


2s\ey + Xk{t) - Xi{t) 
- x,{t)) 


l + 2s ' 


{y)dy 


l + 2s 


2 s|a;fe(f) - Xi{t) 
e'^'’{xk{t) - Xt{t)) 


f u\y) 

Jk 


dy 


|l + 2s • 


2s\xkit) - Xi{t)\ 

We also observe that, for e small, the second member on the right hand side of the 
identity (4.2.13), by using the change of variables (4.2.14), reads 

If ^ 

-J W'(^^Ui{t,x)ju'kit,x)dx 

= - j W (uk{t,x) +'Y^Ui{t,x)'ju'k{t,x) dx 

^ *JR A-i-U 


i^k 


For e small, ul y 


= r W'(u{y) + Yi'^{y + 
Xkit) - Xiit)\ , 


Xkit) - Xiit) 


u'iy) dy. 


is asymptotic either to ui+co) = 1 for Xk > Xi, 
or to u(—cx)) = 0 for Xk < Xi. By using the periodicity of W, it follows that 

N 


J W'(^f]'u-iit,x)^u'kit,x)dx = ^ W'(uiy)^u'iy)dy = Wil)-WiO) = 0, 


again by the asymptotic behavior of u. Concluding, by inserting the results (4.2.15) 
and (4.2.16) into (4.2.13) we get that 

Xkjt) _ yi Xkjt) - Xijt) 

7 i^k M^kit) - x^it)\^'^'^'’ 

which ends the justification of the system (4.2.8). □ 


We recall that, till now, in Theorem 4.2.1 we considered the initial data as a 
superposition of transitions all occurring with the same orientation (see (4.2.7)), i.e. 
the initial dislocation is a monotone function (all the atoms are initially moved to 
the right). 

Of course, for concrete applications, it is interesting to consider also the case 
in which the atoms may dislocate in both directions, i.e. the transitions can occur 
with different orientations (the atoms may be initially displaced to the left or to the 
right of their equilibrium position). 

To model the different orientations of the dislocations, we introduce a parameter 
e { — 1,1} (roughly speaking = 1 corresponds to a dislocation to the right 
and = — 1 to a dislocation to the left). 

The main result in this case is the following (see [115]): 
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Theorem 4.2.2. There exists a viscosity solution of 

dtVe = -^ ( - + O'e) (0> +°0) X 

" 2s ■^ / 0 \ 

^^e(0,x) = -^^^;^cr(0,a;) + 2 forxeR 

such that 

N 

Vv^v^{t,x) = H (^^^{x - x^(t))y 

i=l 

where {xiftf )^ is solution to 


(4.2.17) 


X* = + 2 (0,+oo), 

x,(0) = x°. 


We observe that Theorem 4.2.2 reduces to Theorem 4.2.1 when = ••• = 
= 1. In fact, the case discussed in Theorem 4.2.2 is richer than the one in 
Theorem 4.2.1, since, in the case of different initial orientations, collisions can occur, 
i.e. it may happen that XiiTf) = Xi+i(Tc) for some f e {1,..., TV — 1} at a collision 
time Tc- 

For instance, in the case N = 2, for = 1 and ^2 = ~1 (two initial dislocations 
with different orientations) we have that 

afli+ 2 s 

if cr s: 0 then < ° , 

( 2 s + 1)7 

if 00 < ( 2 s||cr||oo)" 5 ; then ^ —- ° ,, ,, 

7(1 - 2s0o||cr||oo) 

where 0o := X 2 — xj is the initial distance between the dislocated atoms. That is, if 
either the external force has the right sign, or the initial distance is suitably small 
with respect to the external force, then the dislocation time is finite, and collisions 
occur in a finite time (on the other hand, when these conditions are violated, there 
are examples in which collisions do not occur). 


This and more general cases of collisions, with precise estimates on the collision 
times, are discussed in detail in [115]. 


An interesting feature of the system is that the dislocation function does 
not annihilate at the collision time. More precisely, in the appropriate scale, we 
have that at the collision time vanishes outside the collision points, but it still 
preserves a non-negligible asymptotic contribution exactly at the collision points. 
A formal statement is the following (see [115]): 

Theorem 4.2.3. Let N = 2 and assume that a collision occurs. Let Xc be the 
collision point, namely Xc = XiiTf) = X 2 (Tc). Then 

(4.2.18) lim \\iB.Ve{t,x) = 0 for any x A Xc, 

t—yTc e—>-0 

but 

(4.2.19) limsupwe(t, Xc) > 1. 
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Formulas (4.2.18) and (4.2.19) describe what happens in the crystal at the 
collision time. On the one hand, formula (4.2.18) states that at any point that 
is not the collision point and at a large scale, the system relaxes at the collision 
time. On the other hand, formula (4.2.19) states that the behavior at the collision 
points at the collision time is quite “singular”. Namely, the system does not relax 
immediately (in the appropriate scale). As a matter of fact, in related numerical 
simulations (see e.g. [ 1 ]) one may notice that the dislocation function may persists 
after collision and, in higher dimensions, further collisions may change the topology 
of the dislocation curves. 

What happens is that a slightly larger time is needed before the system relaxes 
exponentially fast: a detailed description of this relaxation phenomenon is presented 
in [ 116 ]. For instance, in the case N = 2, the dislocation function decays to zero 
exponentially fast, slightly after collision, as given by the following result: 

Theorem 4.2.4. Let iV = 2, = 1, ^2 = o-nd let be the solution given 

by Theorem 4-2.2, with cr = 0. Then there exist eg > 0, c > 0, > Tc and > 0 
satisfying 


lim Tj = Tc 


and lim o, = 0 

e^O 

such that for any e < cq we have 


(4.2.20) |ne(f,a;)| ^ , for any x G M and t > T^. 


The estimate in (4.2.20) states, roughly speaking, that at a suitable time T^ 
(only slightly bigger than the collision time Tc) the dislocation function gets below 
a small threshold p^, and later it decays exponentially fast (the constant of this 
exponential becomes large when e is small). 

The reader may compare Theorem 4.2.3 and 4.2.4 and notice that different 
asymptotics are considered by the two results. A result similar to Theorem 4.2.4 
holds for a larger number of dislocated atoms. For instance, in the case of three 
atoms with alternate dislocations, one has that, slightly after collision, the disloca¬ 
tion function decays exponentially fast to the basic layer solution. More precisely 
(see again [ 116 ]), we have that: 

Theorem 4.2.5. Let = 3, = ^3 = 1, ^2 = and let be the solution 

given by Theorem 4 - 2 . 2 , with a = 0. Then there exist eg > 0, c > 0, Tf,Tf > Tc 
and Pc > 0 satisfying 


lim T} = lim Tf = Tc, 


and lim p, = 0 
e^O 


and points and Zc satisfying 


lim jzc - yd = 0 
e^O 


such that for any e < Cq we have 


(4.2.21) 



for any x G M and t ^ T} , 


and 


(4.2.22) 



for any x G K and t ^ T^, 
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where u is the basic layer solution introduced in (4.2.5). 

Roughly speaking, formulas (4.2.21) and (4.2.22) say that for times T}, just 
slightly bigger than the collision time T^, the dislocation function Ug gets trapped 
between two basic layer solutions (centered at points and z^), up to a small error. 
The error gets eventually to zero, exponentially fast in time, and the two basic layer 
solutions which trap Ue get closer and closer to each other as e goes to zero (that is, 
the distance between y^ and goes to zero with e). 

We refer once more to [ 116 ] for a series of figures describing in details the 
results of Theorems 4.2.4 and 4.2.5. We observe that the results presented in The¬ 
orems 4.2.1, 4.2.2, 4.2.3, 4.2.4 and 4.2.5 describe the crystal at different space and 
time scale. As a matter of fact, the mathematical study of a crystal typically 
goes from an atomic description (such as a classical discrete model presented by 
Frenkel-Kontorova and Prandtl-Tomlinson) to a macroscopic scale in which a plas¬ 
tic deformation occurs. 

In the theory discussed here, we join this atomic and macroscopic scales by 
a series of intermediate scales, such as a microscopic scale, in which the Peierls- 
Nabarro model is introduced, a mesoscopic scale, in which we studied the dynamics 
of the dislocations (in particular. Theorems 4.2.1 and 4.2.2), in order to obtain at 
the end a macroscopic theory leading to the relaxation of the model to a permanent 
deformation (as given in Theorems 4.2.4 and 4.2.5 , while Theorem4.2.3 somehow 
describes the further intermediate features between these schematic scalings). 

4.3. An approach to the extension problem via the Fourier transform 

We will discuss here the extension operator of the fractional Laplacian via the 
Fourier transform approach (see [31] and [134] for other approaches and further 
results and also [82], in which a different extension formula is obtained in the 
framework of the Heisenberg groups). 

Some readers may find the details of this part rather technical: if so, she or he 
can jump directly to Section 5 on page 59, without affecting the subsequent reading. 

We fix at first a few pieces of notation. We denote points in := K" x 

(0, -too) as X = {x, y), with x e K" and y > 0. When taking gradients in we 

write Vx = C^xjSy), where is the gradient in M". Also, in we will often 

take the Fourier transform in the variable x only, for fixed y > 0. We also set 

a := 1 - 2s e (-1,1). 

We will consider the fractional Sobolev space 7J^(IR") defined as the set of 
functions u that satisfy 

+ [u]g < -1-00, 

where 

Mg := 

For any u e Hjoc ((^i +'*))> we consider the functional 
(4.3.1) G(u) := J t‘^(^\u{t)f + \u'(t)\^^ dt. 

By Theorem 4 of [128], we know that the functional G attains its minimum among 
all the functions u e ((0, -too)) n G°([0, -too)) with u(0) = 1. We call g such 
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minimizer and 

(4.3.2) Cjj := G{g) = min G(u). 

“=”^1 o c « ° ° ( p • + » 

u{0) = l 

The main result of this section is the following. 

Theorem 4.3.1. Let u e iS(K”) and let 

(4.3.3) U{x,y) := J-i(2(0ff(|e|y)). 

Then 


(4.3.4) div(?/“V17) = 0 

for any X = {x,y) e In addition, 


(4.3.5) 


-y%t/ =Gj(-A)*m 


{!/=0} 


in M", both in the sense of distributions and as a pointwise limit. 


In order to prove Theorem 4.3.1, we need to make some preliminary computa¬ 
tions. At first, let us recall a few useful properties of the minimizer function g of 
the operator G introduced in (4.3.1). 

We know from formula (4.5) in [128] that 

(4.3.6) O^g^l, 
and from formula (2.6) in [128] that 

(4.3.7) g' ^ 0. 

We also cite formula (4.3) in [128], according to which g is a solution of 

(4.3.8) g''{t) + at-^g'{t) = g{t) 

for any t > 0, and formula (4.4) in [128], according to which 

(4.3.9) lim t°'g'(t) = — Ch. 

Now, for any V e Wjq^^(K""'"^) we set 

y-\XxV{X)\^dX. 

Notice that \V^a is well defined (possibly infinite) on such space. Also, one can 
compute [y]a explicitly in the following interesting case: 

Lemma 4.3.2. Let ip e iS(K") and 

(4.3.10) U{x,y)-.= X-^{iP{09my))- 


Then 

(4.3.11) 


[Ufa = Ca Wa. 
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Proof. By (4.3.6), for any fixed y > 0, the function ^ '0(^) g(|Cl2/) belongs 

to and so we may consider its (inverse) Fourier transform. This says that 

the definition of U is well posed. 

By the inverse Fourier transform definition (3.1.2), we have that 


.t7(ai,2/) = V, r m9my)e^^-^d^ 

JR^ 

^(*CV'(6ff(ICl2/))(2;)- 


= J- 


Thus, by Plancherel Theorem, 

f \\/^U{x,y)f dx = { \iip{ 09 i\^\y)\^ d^. 

JR^ jRr^ 

Integrating over y > 0, we obtain that 

r y^\^,u{x)\^dx = \ IClXe|'fr”yl5(ieiy)|'rfy 

Jr"+^ Jk" LJo 


(4.3.12) 


r r 

i^r-lV'(ei' tiswr 

^ + 00 r 

t^\9{t)fdt-\ |eP^|i/-(6|'dC 

Jo Je^ 

r+oo „ 

bpfc i“|5(0| dt. 

Jo 


d^ 


d^ 


Let us now prove that the following identity is well posed 
(4.3.13) dyU{x,y) = X-U\^\m9'{ 


For this, we observe that 

(4.3.14) |5'(t)|^C„t-“. 

To check this, we define 7 (t) := \ 9 '(t)\. From (4.3.7) and (4.3.8), we obtain that 

lit) = { 9 " it) + at-^g'it)) = -t“ff(t) < 0. 

Hence 

jit) ^ lim 7 (t) = lim r“| 5 f'(r)| = Cj, 

r-.0+ T^0+ 

where formula (4.3.9) was used in the last identity, and this establishes (4.3.14). 

From (4.3.14) we have that |^| |V’(ei IV'(ei 6 L^(ffi”), 

and so (4.3.13) follows. 

Therefore, by (4.3.13) and the Plancherel Theorem, 
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‘ dy 


dt 


d^ 

di 


Integrating over y > 0 we obtain 

f y''\dyU{x,y)\^ dx = f f y“|5'(IC|y) 

Jh"+^ Jh" LJo 

= f i^r-iV'(oi' 

LJo 

^ + 00 r 

Jo Jr^ 

r+oo „ 

= bpfc dt. 

Jo 

By summing this with (4.3.12), and recalling (4.3.2), we obtain the desired result 
[[/]^ = Cjj [V'Ig- This concludes the proof of the Lemma. □ 

Now, given u e Lj(,j,(IR"'), we consider the space Xu of all the functions V e 
^ioc(®+''"^) that, for any x e K", the map y V{x,y) is in (^^([O,+oo)), 
with V (x, 0) = u{x) for any x e M". Then the problem of minimizing [ • ]a over Xu 
has a somehow explicit solution. 

Lemma 4.3.3. Assume that u e S{M.^). Then 


(4.3.15) 


mm [l/]„ = [U]^ = Ctj [u]g, 


(4.3.16) 


Uix,y) := 


Proof. We remark that (4.3.16) is simply (4.3.10) with ip := u, and by Lemma 
4.3.2 we have that 

ml = 

Furthermore, we claim that 

(4.3.17) UeXu. 

In order to prove this, we first observe that 


(4.3.18) 


\9{T) -9{t)\ ^ 


2s 


To check this, without loss of generality, we may suppose that T ^ t ^ 0. Hence, 
by (4.3.7) and (4.3.14), 

rT 


|5(T)-g(t)| ^ \ 9 '{r)\dr 

‘r 


^ Cj J r-“ dr 
_ Cj (ri-“ - ti-“) 


1 — a 

that is (4.3.18). 

Then, by (4.3.18), for any y, y e (0, +oo), we see that 


5(IC|y)-5(ICIy) < 


2s 
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Accordingly, 


\U{x,y) - U{x,y)\ = 


< 


< 


uiO (sdCly) -5(1^1 y)) 

[ (sdciy)-5(ICIy)) 

-f 

JR^ 


di 


2s 


and this implies (4.3.17). 

Thanks to (4.3.17) and (4.3.11), in order to complete the proof of (4.3.15), it 
suffices to show that, for any V e A„, we have that 

(4.3.19) [F]2 [U]l 


To prove this, let us take V e Xu- Without loss of generality, since [t/]a < +oo 
thanks to (4.3.11), we may suppose that \V^a < +oo- Hence, fixed a.e. y > 0, we 
have that 

f \VxV{x,y)\'^ dx < i/“ f \VxV{x,y)\^ dx < +co, 

JR^ JR^ 


hence the map x g \VxV{x,y)\ belongs to L^(E^). Therefore, by Plancherel Theo¬ 
rem, 


(4.3.20) r \X,V{x,y)\^dx= f I J-(V,H(ai, y)) (0 

J]R»^ JR^ ' 


d^. 


Now by the Fourier transform definition (3.1.1) 


X{X,V{x,y)){0 


hence (4.3.20) becomes 


f V,V{x,y)e-^^<dx 

r i^V{x,y) dx 

Jr" 

i^X{V{x,y)){^), 


(4.3.21) r \X,V{x,y)\^dx= f \J^{V{x,y))m^ d^ 

Jr" Jr" 

On the other hand 

X{dyV{x,y)){0 = dyX{V{x,y)){0 
and thus, by Plancherel Theorem, 

f \dyV{x,y)\^dx= \ \X{dyV{x,y)){ 0 \^d^=\ \dyF{V{x,y)){ 0 ? d^. 

Jr" Jr" Jr" 

We sum up this latter result with identity (4.3.21) and we use the notation ^(^, y) := 
F{y{x,y)){^) to conclude that 


(4.3.22) r \XxV{x,y)\^dx= \ m,yW + \dyH^,y)\" d^. 

Je'^ Je”' 

Accordingly, integrating over y > 0, we deduce that 

(4.3.23) [P]2= \ ^^j^{\em^,y)\^ + \Mi^,y)\")d^dy. 

*'E , 
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Let US first consider the integration over y, for any fixed ^ e K"’\{0}, that we now 
omit from the notation when this does not generate any confusion. We set h{y) := 
We have that h'(y) = therefore, using the 

substitution t = |^| y, we obtain 


J 

JO 


3 24) = dt 

= \e^G{h). 

Now, for any A e M, we show that 

(4.3.25) min w(0) = XG{w) = A^Cjj. 

1«€ ((0.+00)) n C 0 ([ 0, + CX))) 

Indeed, when A = 0, the trivial function is an allowed competitor and 0(0) = 0, 
which gives (4.3.25) in this case. If, on the other hand, A 7 ^ 0, given w as above 
with w(0) = A we set w\{x) := X~^w{x). Hence we see that wa( 0) = 1 and 
thus G{w) = X'^ G{wx) ^ X^ G{g) = A^Cjj, due to the minimality of g. This 
proves (4.3.25). From (4.3.25) and the fact that 

h{o) = ^{^,0) = :f{v{x,o)){ 0 = m, 

we obtain that 

G(/i)>Q|S(e)|^ 

As a consequence, we get from (4.3.24) that 

r+oo 


j 

Jo 


+\dy(l){ty)\ )dy^ C^\i\ 


Integrating over ^ e IR"\{0} we obtain that 

f y“(lCI^ + \dy<l}{^,y)\^) d^dy > Cj [2] 

j]gn+i \ / 

Hence, by (4.3.23), 

[^]a > Q [2]^, 

which proves (4.3.19), and so (4.3.15). 

We can now prove the main result of this section. 


2 

G- 


□ 


Proof of Theorem 4.3.1. Formula (4.3.4) follows from the minimality prop¬ 
erty in (4.3.15), by writing that [1/]^ ^ [1/ -I- eip]^ for any ip smooth and compactly 
supported inside and any e e K. 

Now we take ip e G®(M"') (notice that its support may now hit {y = 0}). We 
define '■= u + ep, and as in (4.3.3), with 2 replaced by 2^ (notice that (4.3.3) 
is nothing but (4.3.16)), hence we will be able to exploit Lemma 4.3.3. 

We also set 

p^{x,y) := (piO 9{\^\y)) ■ 

We observe that 

(4.3.26) (^,(a.,0) = J-i((?(C)g(0)) = (<^(e)) = ^(x) 
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and that 


As a consequence 


U, = U + eJ^-^{0{O9{m) = U + e^,. 


[U,]l = ml + 2e r dX + o(e). 

Jr^+1 

Hence, using (4.3.4), (4.3.26) and the Divergence Theorem, 


(4.3.27) 


ml = ml + 2e f div m y'^XxU) dX + o(e) 
Jh"+i ^ ' 


= ml-‘^^\ (py'^dyUdx + oie). 

JR"x{0} 

Moreover, from Plancherel Theorem, and the fact that the image of ip is in the reals, 
mh = [u]g + 2e r |CP'*u(C) piO d^ + o(e) 

jRri 

= [m]g + 2e J m (x) p{x) dx + o(e) 

= [2]g + 2e {—Ayu{x) p{x) dx + o(e). 

JR" 

By comparing this with (4.3.27) and recalling (4.3.15) we obtain that 

ml-‘2-<^\ py^'dyU dx +o{e) 

JR"x{0} 

= mt 
= GyunTc 

= (7jj[u]G + 2Cjje J (—A)®u(x) (p{x) dx + o(e) 


and so 


= [U]l + 2C^e{ {—Xy u ip dx + o{e) 

Jr" 

— I py°‘dyUdx = C^ I {—Ayupdx, 
JR"xfOl JR" 


for any p 6 (^“(M"), that is the distributional formulation of (4.3.5). 
Furthermore, by (4.3.3), we have that 

y‘^dyU{x, y) = m (id uiO 5(l?l2/)) = ^ (id^-^ 2(d (ld?/)“ <7(l^l2/)) • 

Hence, by (4.3.9), we obtain 

\^ydyU{x,y) = -QJ-i(|d'-“2(d) 


= -c.r 


■'(iep*2(d) 


= - (-A)®m(x), 

that is the pointwise limit formulation of (4.3.5). This concludes the proof of The¬ 
orem 4.3.1. □ 





CHAPTER 5 


Nonlocal phase transitions 


Now, we consider a nonlocal phase transition model, in particnlar described 
by the Allen-Cahn eqnation. A fractional analogne of a conjectnre of De Giorgi, 
that deals with possible one-dimensional symmetry of entire solutions, naturally 
arises from treating this model, and will be consequently presented. There is a 
very interesting connection with nonlocal minimal surfaces, that will be studied in 
Chapter 6. 

We introduce briefly the classical case^. The Allen-Cahn equation has various 
applications, for instance, in the study of interfaces (both in gases and solids), in 
the theory of superconductors and superfluids or in cosmology. We deal here with a 
two-phase transition model, in which a fluid can reach two pure phases (say 1 and 
— 1) forming an interface of separation. The aim is to describe the pattern and the 
separation of the two phases. 

The formation of the interface is driven by a variational principle. Let u{x) be 
the function describing the state of the fluid at position a; in a bounded region fl. 
As a first guess, the phase separation can be modeled via the minimization of the 
energy 



where W is a double-well potential. More precisely, W: [—1,1] ^ [0, oo) such that 


(5.0.28) 


WeC^ ([-1,1]), W{±1) = 0, W > 0 in (-1,1), 
± 1) = 0 and W"{±1) > 0. 


. The classical example is 

- 11 ^ 

(5.0.29) W{u) := ^ ^ W 

On the other hand, the functional in £q produces an ambiguous outcome, since any 
function u that attains only the values +1 is a minimizer for the energy. That is, the 
energy functional in £q alone cannot detect any geometric feature of the interface. 

To avoid this, one is led to consider an additional energy term that penalizes 
the formation of unnecessary interfaces. The typical energy functional provided by 
this procedure has the form 

(5.0.30) £{u) := f W(^u{x)) dx + — f \Vu{x)\'^dx. 

JQ 2 Jf2 

In this way, the potential energy that forces the pure phases is compensated by 
a small term, that is due to the elastic effect of the reaction of the particles. As 
a curiosity, we point out that in the classical mechanics framework, the analogue 


^We would like to thank Alberto Farina who, during a summer-school in Cortona (2014), 
gave a beautiful introduction on phase transitions in the classical case. 
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of (5.0.30) is a Lagrangian action of a particle, with n = 1, x representing a time 
coordinate and u(x) the position of the particle at time x. In this framework 
the term involving the square of the derivative of u has the physical meaning of a 
kinetic energy. With a slight abuse of notation, we will keep referring to the gradient 
term in (5.0.30) as a kinetic energy. Perhaps a more appropriate term would be 
elastic energy, but in concrete applications also the potential may arise from elastic 
reactions, therefore the only purpose of these names in our framework is to underline 
the fact that (5.0.30) occurs as a superposition of two terms, a potential one, which 
only depends on u, and one, which will be called kinetic, which only depends on 
the variation of u (and which, in principle, possesses no real “kinetic” feature). 

The energy minimizers will be smooth functions, taking values between —1 and 
1, forming layers of interfaces of e-width. If we send e ^ 0, the transition layer will 
tend to a minimal surface. To better explain this, consider the energy 

(5.0.31) J{u) = J i I Vup + W{u) dx, 

whose minimizers solve the Allen-Cahn equation 
(5.0.32) - Au-t W'(u) = 0. 

In particular, for the explicit potential in (5.0.29), equation (5.0.32) reduces (up to 
normalizations constants) to 

(5.0.33) —Au = u — u^. 

In this setting, the behavior of u in large domains reflects into the behavior of the 
rescaled function Ug{x) = u(-) in Bi. Namely, the minimizers of J in are the 
minimizers of Jg in Bi, where Jg is the rescaled energy functional 

(5.0.34) Je{u) = f +-W(zi) da;. 

Jsi ^ £ 

We notice then that 

Je{u)^ ( V'2VP(m) |Vu|da; 

Jsi 

which, using the Co-area Formula, gives 

Je{u) > J ■sj2W{t) ({u = t}) dt. 

The above formula may suggest that the minimizers of Jg have the tendency to 
minimize the (n — l)-dimensional measure of their level sets. It turns out that 
indeed the level sets of the minimizers of Jg converge to a minimal surface as e —> 0: 
for more details see, for instance, [121] and the references therein. 

In this setting, a famous De Giorgi conjecture comes into place. In the late 70’s, 
De Giorgi conjectured that entire, smooth, monotone (in one direction), bounded 
solutions of (5.0.33) in the whole of K" are necessarily one-dimensional, i.e., there 
exist Lo e and uq • ® ^ ® such that 

u{x) = uo(w • x) for any x G K". 

In other words, the conjecture above asks if the level sets of the entire, smooth, 
monotone (in one direction), bounded solutions are necessarily hyperplanes, at least 
in dimension n ^ 8. 
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One may wonder why the number eight has a relevance in the problem above. 
A possible explanation for this is given by the Bernstein Theorem, as we now try 
to describe. 

The Bernstein problem asks on whether or not all minimal graphs (i.e. surfaces 
that locally minimize the perimeter and that are graphs in a given direction) in K" 
must be necessarily affine. This is indeed true in dimensions n at most eight. On 
the other hand, in dimension n > 9 there are global minimal graphs that are not 
hyperplanes (see e.g. [87]). 

The link between the problem of Bernstein and the conjecture of De Giorgi 
could be suggested by the fact that minimizers approach minimal surfaces in the 
limit. In a sense, if one is able to prove that the limit interface is a hyperplane 
and that this rigidity property gets inherited by the level sets of the minimizers Ug 
(which lie nearby such limit hyperplane), then, by scaling back, one obtains that the 
level sets of u are also hyperplanes. Of course, this link between the two problems, 
as stated here, is only heuristic, and much work is needed to deeply understand the 
connections between the problem of Bernstein and the conjecture of De Giorgi. We 
refer to [74] for a more detailed introduction to this topic. 

We recall that this conjecture by De Giorgi was proved for n ^ 3, see [86, 12, 5] . 
Also, the case 4 ^ n ^ 8 with the additional assumption that 

(5.0.35) lim u(a;', a:„) = ±1, for any a;'e 

Xn—* + OD 

was proved in [120]. 

For n > 9 a counterexample can be found in [55]. Notice that, if the above 
limit is uniform (and the De Giorgi conjecture with this additional assumption is 
known as the Gibbons conjecture), the result extends to all possible n (see for 
instance [73, 74] for further details). 

The goal of the next part of this book is then to discuss an analogue of these 
questions for the nonlocal case and present related results. 

5.1. The fractional Allen-Cahn equation 

The extension of the Allen-Gahn equation in (5.0.32) from a local to a nonlocal 
setting has theoretical interest and concrete applications. Indeed, the study of long 
range interactions naturally leads to the analysis of phase transitions and interfaces 
of nonlocal type. 

Given an open domain D c K" and the double well potential W (as in (5.0.29)), 
our goal here is to study the fractional Allen-Gahn equation 

(—A)®u -I-W'(m) = 0 in D, 

for s e (0,1) (when s = 1, this equation reduces to (5.0.32)). The solutions are the 
critical points of the nonlocal energy 

(5.1.1) £{u,n):= ^^W{u{x))dx+^ JJ dx dy, 

R2"\(nC)2 

up to normalization constants that we omitted for simplicity. The reader can com¬ 
pare (5.1.1) with (5.0.30). Namely, in (5.1.1) the kinetic energy is modified, in 
order to take into account long range interactions. That is, the new kinetic en¬ 
ergy still depends on the variation of the phase parameter. But, in this case, far 
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away changes in phase may influence each other (though the influence is weaker and 
weaker towards inflnity). 

Notice that in the nonlocal framework, we prescribe the function on x fl'' and 
consider the kinetic energy on the remaining regions (see Figure 9). The prescription 
of values in fl'' x Vp reflects into the fact that the domain of integration of the 
kinetic integral in (5.1.1) is M^”\(n‘')^. Indeed, this is perfectly compatible with 
the local case in (5.0.30), where the domain of integration of the kinetic term was 
simply To see this compatibility, one may think that the domain of integration 
of the kinetic energy is simply the complement of the set in which the values of 
the functions are prescribed. In the local case of (5.0.30), the values are prescribed 
on or, one may say, in then the domain of integration of the kinetic energy 
is the complement of fl*", which is simply il. In analogy with that, in the nonlocal 
case of (5.1.1), the values are prescribed on x = (fl*")^, i.e. outside fl for 
both the variables x and y. Then, the kinetic integral is set on the complement 
of {Tff, which is indeed 

Of course, the potential energy has local features, both in the local and in 
the nonlocal case, since in our model the nonlocality only occurs in the kinetic 
interaction, therefore the potential integrals are set over O both in (5.0.30) and 
in (5.1.1). 

For the sake of shortness, given disjoint sets A, B K" we introduce the 
notation 


<a,b):= r r 

JA JB 


\u{x)-u{y)\^ 

n+2s 


lAJB \x-y\ 
and we write the new kinetic energy in (5.1.1) as 


dx dy, 


(5.1.2) 


/C(u,0) 


1 

2 


M(n, n) + u{vl, dP). 



Let us define the energy minimizers and provide a density estimate for the mini- 
mizers. 
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Definition 5.1.1. The function u is a minimizer for the energy £ in Bn if 
£{u,Bn) ^ £{v,Bn) for any v such that u = v outside Bn- 

The energy of the minimizers satisfy the following uniform bound property on 
large balls. 


Theorem 5.1.2. Let u be a minimizer in Bn +2 for a large R, say i? > 1. Then 


(5.1.3) 

More precisely, 


=Q. 

i?—*- + 00 


£{u,Bn) < 


if sG(il), 
logi? if s = |, 

if s e (o, . 

Here, C is a positive constant depending only on n, s and W. 


Notice that for s e 




> i?” These estimates are optimal (we 


refer to [125] for further details). 


Proof. We introduce at first some auxiliary functions. Let 

'4’{x) := —1 + 2 min |(|a;| — R— l) + ,l|, v{x) := min |u(a;), '0(a;)|, 

d{x) := max|(i? + 1 — |a;|), l|. 



Figure 10. The functions ip, v and d 


Then, for |a; — ^ d{x) we have that 

(5.1.4) 
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Indeed, if |a;| ^ R, then d{x) = i? + 1 — |a;| and 

\y\ < k - 2/1 + |a;| ^ d{x) + |a;| i? + 1, 

thus = tpiy) = 0 the inequality is trivial. Else, if |a;| > R, then d{x) = 1, 
and so the inequality is assured by the Lipschitz continuity of tp (with 2 as the 
Lipschitz constant). 

Also, we prove that we have the following estimates for the function d: 


(5.1.5) 


L 


d{x)~ 



'CR^-^ 

if 

s e 

(A) 


CR^-^\ogR 

if 

s = 

1 

2 ’ 



if 

s e 

(A) 


To prove this, we observe that in the ring Bh+tXBh, we have d{x) = 1. Therefore, 
the contribution to the integral in (5.1.5) that comes from the ring is 

bounded by the measure of the ring, and so it is of order namely 


(5.1.6) 


£ 


d{x)-^^dx = \Br+2\Br\ ^ CR^-\ 

^Bfi+2\BR 

for some C > 0. We point out that this order is always negligible with respect to 
the right hand side of (5.1.5). 

Therefore, to complete the proof of (5.1.5), it only remains to estimate the 
contribution to the integral coming from Br. 

For this, we use polar coordinates and perform the change of variables t = 
p/{R + !)• In this way, we obtain that 

''R ^n—1 


r d{x) dx=c r 

J Bji J 0 


- 1 - 


dp 


R+i 

= C(i?+l)"-2M 

Jo 

ri“HTT 

^ C (i? + 1)"-2H (i_i) 

Jo 


-2s 


dt 


dt, 


for some C > 0. Now we observe that 

r el 


1 1 
R+1 

a-ty 

Jo 


' dt ^ ^ 


r (i-t)-2 

Jo 

-log(l - 1) 


'■dt = C 


^ logi? 


if 

if 


s e 


s = 


2 ’ 

if s e (" i 


1 

2 ’ 


The latter two formulas and (5.1.6) imply (5.1.5). 

Now, we define the set 

A\= {v = Ip} 

and notice that Br+i c A c Br+ 2 - We prove that for any x e A and any y e A^ 
(5.1.7) 111 ( 2 ;) - v{y)\ s: max||ii( 2 :) - u{y)\, \ip{x) - i/>(y)||. 

Indeed, for a; e A and y e A^ vje have that 

v(x) = ipix) ^ u{x) and v{y) = u{y) ^ ip{y), 
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therefore 

v(x) — v{y) ^ u{x) — u{y) and v{y) — v(x) ^ 'ip{y) — 
which establishes (5.1.7). This leads to 
(5.1.8) v{A,A^)!iu{A,A^) + ip{A,A^). 

Notice now that 

£{u,Br+2) < £{v,Br+2) 

since u is a minimizer in i3_R+2 and v = u outside i?R+ 2 - We have that 

£{u,Bfij^ 2 ) = 2 ^(^R+ 2 ,^R+ 2 ) + u(i?i^+ 2 , i?£+ 2 ) +1" W{u) 

J BFi+2 

~ 2 

+ 2 ^ii:+ 2 \ 2 l) + u{Bji^ 2 \A, 73^+2) 


da; 


+ r W(w)da;+ f W{u)dx. 

I o\A 


Sf£ + 2\^ 

Since u and v coincide on A^, by using the inequality (5.1.8) we obtain that 

0 ^ £{v,Bii+2) - £{u,Br+2) 

1 


~ 2^ 
1 

2 ^ 


{A, A) — ^^(A, A) + v{A, A^) — u{A, A^) + J ^Vb(i;) — W{u^ dx 

{A, A) - ^u{A, A) + ip{A, A^) + J (w{v) - W(u)) dx. 

Moreover, v = tjj on A and we have that 

l-u{A,A)+ r W{u) dx ^ tl;{A,A) + tlj{A,A^) + f W{ilj)dx = £{tp,A), 

^ JA ^ JA 

and therefore, since Sr+i ^ A^ 7?r_|_2, 

(5.1.9) ^ u(Br+i,Br+i) + r W{u)dx £{tp,Bii+ 2 )- 

We estimate now £, B^j^ 2 )■ For a fixed x 6 73^+2 we observe that 

\^{x) - i^iy)? 


J 

Je- 


I 


a;-i/|"+2* 

|V;(a;) -i/;(y)p 

\x—y\^d(x) y\ 


dy + 


I 


\i>{x) -'ip{y)\^^ 

\x—y\^d{x) y\ 


dy 


^ C 


(^I 


d(^X^ j\x—y\^d{x) 


\^-y\ 


—n—2s+2 


dy + 


I 


\x-y\^d{x) 


\x-y\ 


—n—2s 


dy 


where we have used (5.1.4) and the boundedness of ijj. Passing to polar coordinates, 
we have that 


J 

Je^ 


|V;(x) -i/;(y)p 
\x - y|"+2s 


dy^C 


1 

d{xy 


^d{x) 

Jo 


„-2s + l 


dp + 


rCC 

J d(x) 


dp 


= Cd{x)- 
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Recalling that = —1 on Rr+i and W{—1) = 0, we obtain that 


£(!/>, Rr+ 2)= f f 

< f d{x)- 

^ B R J.O 


\x - y|"+2s 

I 


dy dx + 


r wi^) 

Br-^2 


dx 


^dx + 


Br + 2\Bh+1 


W{'tp) dx. 


Therefore, making use of (5.1.5), 





if 

s e 

(1‘) 

(5.1.10) 

£{iP,Br+2) sS ^ 

CR"-ilogR 

CRri-2s 

if 

if 

s = 

s e 

(M) 


For what regards the right hand-side of inequality (5.1.9), we have that 


lu( 


(5.1.11) 


2 

+ r W{u)dx. 

JBr 

We prove now that 



(5.1.12) 

u(^R,S^+i)< f d{x) dx. 

For this, we observe that if a: e Hr, 
and y e i?R_,_ i , then 

then d{x) = i? + 1 — a; 


1 

W 

1 

w 

II 

1 

1—1 

+ 

A\ 


Therefore, by changing variables z = x — y and then passing to polar coordinates, 
we have that 


p-^^-Up 


) <4 j 

dx 


Br J 

r 

pOO 

^ C 

dx 



'Br 

'd{x) 

= c 

d{x)-^^ 

'Br 


dx. 


This establishes (5.1.12). 

Hence, by (5.1.5) and (5.1.12), we have that 


(5.1.13) w(Rr,R^+i) ^ r d{x)-^^dx^< 

J Br^o 




if 

s e 

(M) 


CR^-^logR 

if 

s = 

1 

2 ’ 


CRn-2s 

if 

s e 

(o.i) 


We also observe that, by adding to inequality (5.1.11), we obtain that 

)B„ 


2^ 


(Rr+ 1 , Rr+i) + r W{u) dx + u{B[i, B^_^i) 

{BR,Bji) + u{B]i,Bii+i\Bji) + r W{u) dx + u{B]i,B‘^^i) 

JBr 


1 

»2“' 


= £{u,Br). 
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This and (5.1.9) give that 

£{u,Bp:) £{%!}, Bii+ 2 ) + 

Combining this with the estimates in (5.1.10) and (5.1.13), we obtain the desired 
result. □ 


Another type of estimate can be given in terms of the level sets of the minimizers 
(see Theorem 1.4 in [125]). 

Theorem 5.1.3. Let u be a minimizer of £ in B^. Then for any 61,62 e (—1,1) 
such that 

u{0) > 6»i 

we have that there exist R and C > 0 such that 

|{u> 02} > CR^ 

if R ^ i?(0i,02)- The constant C > 0 depends only on n, s and W and i?(0i,02) Is 
a large constant that depends also on 9i and 02 - 

The statement of Theorem 5.1.3 says that the level sets of minimizers always 
occupy a portion of a large ball comparable to the ball itself. In particular, both 
phases occur in a large ball, and the portion of the ball occupied by each phase is 
comparable to the one occupied by the other. 

Of course, the simplest situation in which two phases split a ball in domains 
with comparable, and in fact equal, size is when all the level sets are hyperplanes. 
This question is related to a fractional version of a classical conjecture of De Giorgi 
and to nonlocal minimal surfaces, that we discuss in the following Section 5.2 and 
Chapter 6. 

Let us try now to give some details on the proof of the Theorem 5.1.3 in the 
particular case in which s is in the range (0,1/2). The more general proof for all 
s e (0,1) can be found in [125], where one uses some estimates on the Gagliardo 
norm. In our particular case we will make use of the Sobolev inequality that we 
introduced in (3.2.1). The interested reader can see [122] for a more exhaustive 
explanation of the upcoming proof. 


Proof of Theorem 5.1.3. Let us consider a smooth function w such that 
w = 1 on (we will take in sequel re to be a particular barrier for u), and define 

v{x) := min{it(a;),r(;(a;)}. 

Since |m| ^ 1, we have that v = u in B^. Calling D = (K" x K”) \ (i?)^ x B’^) we 
have from definition (5.1.2) that 

lC{u - V, Br) + K,{v, Br) - IC{u, Br) 

1 ff |(m- u)(a;) - (u- '(;)(?/)p + |u(x) - w(j/)p - |u(x) - m(?/)P 


il 

D 


\X - ?/|"+2s 

,|2 , 1,2 _ 


■ dx dy. 


We use the algebraic identity |a — + r — = 26(6 — a) with a = u(x) — u{y) 

and 6 = v{x) — v{y) to obtain that 

]C{u-v,Br) + 1C{v,Br) - 1C{u,Br) 

{{u - v)(x) -(u- v){y)) {v{y) - v(x)) 


JJ 

D 


\x-y\ 


n+2s 


■ dx dy. 
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Since u — v = 0 on we can extend the integral to the whole space M" x M”, 
hence 


IC{u-v,Br) + - IC{u,Br) 

((u — v){x) — (u - 


r 


v){y)) {v{y) - v(x)) 


-y\ 


n+2s 


dx dy. 


Then by changing variables and using the anti-symmetry of the integrals, we notice 
that 


ff ((R-r;)(a:) - {u-v){y)) {v{y) - v{x)) 

BjiXBn 

JJ |x-y|"+2« 

BuxBn 

_ ff 

JJ |x-y|"+2^ 


BrxBji 


= 2 


■* r 


{u - v)(ai) {v{y) - i;(x)) 


\x - y\ 


n+2s 


dx dy 


BrXBr 


and 


II 


{{u - v){x) -{u- v){y)) {v{y) - v(x)) 


dx dy 




+ 


II 


\X - ?/|"+2s 

{{u - v)(x) -{u- v)(y)) (v(y) - v(x)) 


B^xBr 


\x-y\ 


n+2s 


dx dy 


II 


(m - v)(x) (v(y) - v(x)) 


BrxB% 


\x - y\ 


n+2s 


dx dy 


= 2 


^ ff (»-■■)(!/) 

JJ |a:-y|"+2^ 

B%xBr 

JJ |a:-y|"+2^ 

BrxB% 
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Therefore 


IC{u-v,Bfi) + 1 C{v,Br) - K,{u,Bii) 

(u(x) - v(x)) {v{y) - v(x)) 

\X - 

v{y) - v{x) 


= 2 


JJ 


2 I r 

\JIF 


\x-y\ 


n+2s 


dx dy 
dy ) dx 


= 2 

^ 2 

= 2 


r (M(a;) - w(a:)) ( f dy \ dx 

JBRn{u>v=w} VJr" \x-yr+^-^ J 


X - 2 /|"+ 2 « 

I 


Si?n{u>to} 


(u — w){x) (—(—A)^w) (x) dx. 


Hence 


‘ Bjir\{u>w} 


lC{u — V, Bfj) 

^]C{u,Bii)-]C{v,Bii) + 2 r 

Ji 

By adding and subtracting the potential energy, we have that 
JC{u - V, Bn) 


(u — w) (—(—A)^w) dx. 


< 


‘ Bfir\{u>w} 

and since u is minimal in Bn, 


S{u, Bn) - £{,v, Br) + r Wiv) - W{u) 
Jbr 

+ 2 J 

JB 


dx 


(u — w) (—(—A)®w) dx 


(5.1.14) 


lC{u — v,Bn)^ j W{w) — W{u)dx 

J Bjin{u>w=v} 


i 


Bjir\{u>w} 


(u — w) (—(—A)®rt;) dx. 


We deduce from the properties in (5.0.28) of the double-well potential W that there 
exists a small constant c > 0 such that 

W{t) — W{r) > c(l -I- r){t — r) -I- c{t — r)^ when —1-l-c 

Wir) — W (t) ^ ——— when —l^r^t^l. 

c 

We fix the arbitrary constants 9i and 6 * 2 , take c small as here above. Let then 


:= minimi, 02,-1 + c}. 
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It follows that 
(5.1.15) 


I 


W{w) — W{u)dx 




BRr\{u>w} 

j W{w) — W{u)dx+ j 

J Bjir\{0-t>u>w} 

; (1 — w){u — w) dx — c 

JBjir<{6ir>u>w} J Biiri{6i, 

(1 + w) dx 

^ JBRn{u>max.{9ir,w}} 

^ — c r (1— w){u—w) dx+- r 

J ^ Jl 
Therefore, in (5.1.14) we obtain that 


W{w) — W{u)dx 


^ — c 
1 

+ - 


>u>w} 


{u — w)'^ dx 


B n {lA > max { 0 *, 10 } } 


(1 + w) dx. 


JC{u — V, Bft) ^ — c {1 — w){u — w) dx 

J BRr\{6i,>u>w} 


(5.1.16) 


1 


I 


^ d B]in{u>max{9.t, ,w}} 


(1 + w) dx 


+ 2 


I 


B pir\{u>w} 


{u — w) (—(—A)®u;) dx. 


We introduce now a useful barrier in the next Lemma (we just recall here Lemma 
3.1 in [125] - there the reader can find how this barrier is build): 

Lemma 5.1.4. Given any r 0 there exists a constant C > 1 (possibly depend¬ 
ing on n, s and t) such that: for any R ^ C there exists a rotationally symmetric 
function w e C (K", [—1 + 1]) with w = 1 in and such that for any 

X e Bfi one has that 

(5.1.17) + 

o 


(5.1.18) —(—A)®rt;(a;) ^ t(1 + w(x)). 

Taking w as the barrier introduced in the above Lemma, thanks to (5.1.16) and 
to the estimate in (5.1.18), we have that 


K,{u — V, Bfi) ^ — c {1 w){u — w) dx 

J BRr\{9ir>u>w} 

(1 + w) dx 


-f 

c Jl 


+ - 
c 

+ 2r 


Sfjo {ii>max{@* ,ii;}} 

(u — ?r;)(l + w) dx. 


I 


Bjin(u>w} 
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Let then t = |, and we are left with 


ic{u - v,B[i) ^ c r 

Ji 


B{u>ma.x.{6 ,io}} 
1 


(u — + w) dx 


I 


^ ai5i^n{u>max{0*,u>}} 


(1 + w) dx 


^ Cl 


L 


(1 + w) dx, 

with Cl depending on c (hence on W). Using again Lemma 5.1.4, in particular the 
right hand side inequality in (5.1.17), we have that 


1C{u-v,Br) ^ Cl - C r 

Ji 


B n {lA > max { 0 *, 10 } } 


(i?+l-|x|) 


-2s 


We set 

(5.1.19) V{R) :=\BRr,{u> e^}\ 
and the Co-Area formula then gives 

rR 

(5.1.20) IC{u-v,Br) ^C2 \ {R+l-t)-^^V'{t)dt, 

Jo 

where C 2 possibly depends on n, s, W. 

We use now the Sobolev inequality (3.2.1) for p = 2, applied to it — u (recalling 
that the support of it — u is a subset of Br) to obtain that 

|(ii-i;)(x) - (u-i;)(2/)p 


IC{u-v,Br) = IC{u-v,R^) = JJ 


(5.1.21) 

From (5.1.17) one has that 


\x - y\ 


n+2s 


■ dx dy 


> Clllt — up 2,1 = Cllli — up 2,. 


ii;(x)^C(i?-tl-|x|)-2"-l. 
We fix K large enough so as to have R > 2K and in Br-r 

w{x) ^ C(1 -t iF)-2" - 1 ^ -1 -t ^ ^ 
Therefore in Br-r n {it > 0*} we have that 


lit — ul ^u — w^u+1 — 


1 


> 


2 

1 + 0* 


Using definition (5.1.19), this leads to 


It — u 2 ^ 

(Br) 


I 


\u — v\ ^- 2 ^ dx 


:> 


Br 

1 + ^ 




dx 


n{u>S,} 


> CsViR-K)'- 


In (5.1.21) we thus have 


JC{u-v,BR)^C3V{R-Ky 
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and from (5.1.20) it follows that 

CiViR-K)"^ \ {R+l-t)-'^^V'{t)dt. 

Jo 

Let R^ 2K. Integrating the latter integral from p to ^ we have that 

r ~ 

Ci^V{p-K)'^ <C4 I ' V{R- dR 

2 Jp 

/ pR \ 

^ J \J (^+ 


dR 


3 p / ^ \ 

I' y'{t) (I' (R+i-ty^^dR] 

Jo 


dt 


3p 


1 - 2s 

l-2s 


dt. 


1 ^ (2/9)^ hence, 


I dt 


Since 1 — 2s > 0, one has for large p that ^ ~ ^ 

noticing that the function V is nondecreasing, 

P-V{p-K)^ ^ £ V'{t ), 

^ C^p^-^^V{2p). 

Therefore 

(5.1.22) p^^V{p-^2C^V{2p). 

Now we use an inductive argument as in Lemma 3.2 in [125], that we recall here: 

Lemma 5.1.5. Let cr, /r e (0, oo), v e (cr, co) and 7 , Rq, C e (1, co). 

Let V (0, 00 ) —> (0, 00 ) be a nondecreasing function. For any r e [i?o,oo), let 

a{r) := m 


log V (r) 1 

1, —^->. Suppose that V(Ro) > u and 

logr J 


r^air)V{r) •' ^ CV{'yr), 

for any r e [i?o, 00). Then there exist c G ( 0 , 1 ) and i?* G [i?o, 00), possibly depending 
on /i, "f, Rq,C such that 

V{r) > cr‘', 

for any r G [i?*, 00). 

For R large, one obtains from (5.1.22) and Lemma 5.1.5 that 

V{R) > coi?", 

for a suitable cq G (0,1). Let now 

9* := max{0i, 02 , —1 + c}. 

We have that 

\{u>e*}nBR\ + \{9.<u<e*}nBR\ 

= |{u > 0.} n Bn] 

= V{R) > coi?”. 


(5.1.23) 
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Moreover, from (5.1.2) we have that for some c > 0 

£{u,Bn) 

therefore 

> £{u, Br) > r W{u) dx 

^ {9i,<u<i9*]r\BR 

> inf W{t) |{0* < u < 9*} n Br\. 

■,9*') 

From this and (5.1.23) we have that 

coi?" < + |{u > 9*} n Br\, 


and finally 

\{u>9*}nBR\^CR”, 

with C possibly depending on n,s,W. This concludes the proof of Theorem 5.1.3 
in the case s e (0,1/2). □ 


5.2. A nonlocal version of a conjecture by De Giorgi 


In this section we consider the fractional counterpart of the conjecture by De 
Giorgi that was discussed before in the classical case. Namely, we consider the 
nonlocal Allen-Cahn equation 

-(-A)®u + lF(u) = 0 in M”, 


where IT is a double-well potential, and u is smooth, bounded and monotone in 
one direction, namely |m| ^ 1 and dx^u > 0. We wonder if it is also true, at least in 
low dimension, that u is one-dimensional. In this case, the conjecture was initially 
proved for n = 2 and s = ^ in [27]. In the case n = 2, for any s e (0,1), the result is 
proved using the harmonic extension of the fractional Laplacian in [26] and [132]. 

For n = 3, the proof can be found in [24] for s e ij ■ The conjecture is still open 

for n = 3 and s e |^0, and for n ^ 4. Also, the Gibbons conjecture (that is the 
De Giorgi conjecture with the additional condition that limit in (5.0.35) is uniform) 
is also true for any s e (0,1) and in any dimension n, see [75] . 

To keep the discussion as simple as possible, we focus here on the case n = 2 and 
any s e (0,1), providing an alternative proof that does not make use of the harmonic 
extension. This part is completely new and not available in the literature. The proof 
is indeed quite general and it will be further exploited in [43] . 

We define (as in (5.1.2)) the total energy of the system to be 


(5.2.1) 


£{utBr) — JCr 



W{u)dx, 


where the kinetic energy is 


(5.2.2) 


i JJ 

Qr 


|M(a:) — u{x)\^ 

[2, _ ^|n+2s 


dx dx, 


and Qr := = {Br x Br) u {Br x (M"\Bk)) u {{W^\Br) x Br). We 

recall that the kinetic energy can also be written as 


K-r{u) = -u{Br, Br) + u{Br, B^), 


( 5 . 2 . 3 ) 
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where for two sets A, B 


(5.2.4) 


iiA,B)= r r 

JA JB 


\u{x) — u{x)\'^ 
\x - x|"+2s 


dx dx. 


The main result of this section is the following. 

Theorem 5.2.1. Let u be a minimizer of the energy defined in (5.2.1) in any 
hall o/K^. Then u is 1-D, i.e. there exist ui e and uq : M ^ K such that 


u{x) = uo{oj ■ x) for any a; e 

The proof relies on the following estimate for the kinetic energy, that we prove 
by employing a domain deformation technique. 

Lemma 5.2.2. Let R > 1, cp e C®(i?i). Also, for any y e K", let 

(5.2.5) 4'fl,+ (y) := y + ei and 'i/R_{y):=y-(p(^^^ei. 

Then, for large R, the maps 4'^^+ and are dijfeomorphisms on M". Further¬ 

more, if we define ur^+{x) := u{'i’]^\{x)), we have that 


(5.2.6) 


K-r{ur,+ ) + ICr{ur-) - 2K,r{u) ^ —1Cr{u), 


for some C > 0. 


Proof. First of all, we compute the Jacobian of '^r^±. For this, we write 
to denote the i^^ component of the vector 4 '_r,+ = (4>ij_+i,--- _!_„,) and we 

observe that 


(5.2.7) 


d'i'R^+piy) 


d 




The latter term is bounded by 0{R~^), and this proves that + is a diffeomor- 
phism if R is large enough. 

For further reference, we point out that if Jr^± is the Jacobian determinant 
of 4 '_r,+, then the change of variable 


(5.2.8) 


X := 4'R,+(y), 


X ■■= 4'R,+ (y) 


gives that 


dxdx = jR^+{y) jR^+{y)dydy 

- (‘ ± +°(i 5 )) (‘ ± hMl )+ 

- ‘ ± yHii) - yAii) + °(i) 
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thanks to (5.2.7). Therefore 


(5.2.9) 


dx dx 


\x — ccj 


\x — 

|u(^'^^+(a;)) -m(T^^+(x))|" 


\u{y) — u{y)\'^ I ^R,±iy) ^R,±{y) 


\y_y\n + 2s 
1 


dx dx 


|y- 2 / 1 ^ 


1 + + of 

~ R ^^\r) ~ R ^^\r) V 



dydy. 


Now, for any y, y s K" we calculate 


(5.2.10) 


'^R,±{y) - ^R.+(y) 


{y-y)±[ 


= I2/-2/P + 


T’ 




ei 


10 / / y 


-^[j^)]iyi-yi)- 


Notice also that 


(5.2.11) 


^^r)-^\r 


< ■j^Mc^(R’^)\y-y\, 


hence (5.2.10) becomes 


'^R,±{y) - ^'i?..+ (y) 


\y-y\^ 


1 + 77 + 


where 


(5.2.12) 


v± ■■= 



yi) 



As a consequence 


'^R,±iy) - '^R,±iy) 


\y-yf^ 


= (1 + 77 +)' 


= 1 - 


n + 2 s 


77 + + C>(77|)- 
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We plug this information into (5.2.9) and use (5.2.12) to obtain 
|MR.+ (a;) - ur ^+{ x )\'^ 


. _ ^|n+2s 


|M(y) - uiy)\ 

\y-y\n+2s 

1 


1 - 


■ dx dx 
n + 2s 


2 + 

1 


1 ± ± + “^(d)) 


\u{y) - u(j/)|^ 


\y_y\n + 2s 


n + 2s 


ri±+ ±RMI)±7iMI) 


+0 


dydy. 


Using this and the fact that 

77+ + 77_ = 2 


-Hr 


= O 


(d)’ 


|y-y|2 \R‘^ 

thanks to (5.2.11), we obtain 

|uR,+(a;) - Ufl,+(i)p \uR-{x) - UR^_{x)\'^ 


+ 


\X - x|"+2'* 

\u{y)-u{y)\^ , 

\y_y\n + 2s 

Thus, if we integrate over Q-R we find that 


\x — X 


pp.\n-\-2s 


dx dx 


2 + C’(^) ) dydy. 


ICr{ur.+) +ICr{ur.+) = 2K.r{u) + 

Qr 

This establishes (5.2.6). 


□ 


Proof of Theorem 5.2.1. We organize this proof into four steps. 

Step 1. A geometrical consideration 

In order to prove that the level sets are flat, it suffices to prove that u is monotone 
in any direction. Indeed, if u is monotone in any direction, the level set {u = 0} is 
both convex and concave, thus it is flat. 


Step 2. Energy estimates 

Let tp G C^{Bi) such that </? = 1 in i?i/2) and let e = (1, 0). We define as in Lemma 
5.2.2 

'^R.+{y) ■= y + e and 'i>R_{y):=y-(p(^^'^e, 

which are diffeomorphisms for large R, and the functions ur^+{x) := u(5')^^_,_(x)). 
Notice that 

for ye B^^r 
for j/G Br/2. 


(5.2.13) 

(5.2.14) 


UR,+{y) = u{y) 
UR,+{y) = u{y - e) 
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By computing the potential energy, it is easy to see that 

j W{u]i,+ {x)) dx + I W{un.-{x)) dx — 2 j W{u{x))dx 

*jBji 

Using this and (5.2.6), we obtain the following estimate for the total energy 

(5.2.15) £{ur^+,Bii) + £{ur^_,Br) - 2£{u,Br) ^ ^£{u,Br). 

Also, since ur^+ = u in i?^, we have that 

£ (it, Br) ^ £ {ur^_ ,Br). 

This and (5.2.15) imply that 

(5.2.16) £{ur^+,Br) - £{u,Br) ^ ^£{u,Br). 

As a consequence of this estimate and (5.1.3), it follows that 

(5.2.17) lim (£iuR,+ ,BR)-£iu,BR)) =0. 

R^+oo V / 

Step 3. Monotonicity 

We claim that it is monotone. Suppose by contradiction that it is not monotone. 
That is, up to translation and dilation, we suppose that the value of it at the origin 
stays above the values of e and —e, with e := (1,0), i.e. 

ii(0) > u(e) and it(0) > it(—e). 

Take R to be large enough, say R > 8. Let now 

(5.2.18) i;ij(a;) := min {it(a;), itu_+(a;)} and it;fl;(x) := max {it(x), iti^^+(x)}. 

By (5.2.13) we have that vr = wr = u outside Br. Then, since u is a minimizer in 
Br and wr = it outside Br, we have that 

(5.2.19) £{wR,BR)^£iu,BR). 

Moreover, the sum of the energies of the minimum and the maximum is less than 
or equal to the sum of the original energies: this is obvious in the local case, since 
equality holds, and in the nonlocal case the proof is based on the inspection of the 
different integral contributions, see e.g. formula (38) in [114]. So we have that 

S{vr,Br) + £{wr,Br) s: £{u,Br) + £{ur^+,Br) 

hence, recalling (5.2.19), 

(5.2.20) £{vr, Br) £{ur^+, Br). 

We claim that vr is not identically neither u, nor ur^+. Indeed, since u(0) = 
UR^+{e) and u(—e) = ur,+(0) we have that 

vr(0) = min{u(0),iiH,+ (0)} = min{it(0),ii(-e)} 

= u{—e) = ur^+{0) < u(0) and 
VR{e) = min{u(e),itii:,+(e)} = min{u(e),it(0)} 

= u(e) < it(0) = it_R,+(e). 
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(5.2.21) 


lim (5 r = 0. 

R—►+0O 


By continuity of u and ur,+, we have that 

vr = ur^+ < u in a neighborhood of 0 and 
vr = u < ur,+ in a neighborhood of e. 

We focus our attention on the energy in the smaller ball i? 2 . We claim that vr is 
not minimal for £(-,B2). Indeed, if vr were minimal in B2, then on B2 both vr, 
and u would satisfy the same equation. However, vr ^ u in by definition and 
Vr = u in a neighborhood of e by the second statement in (5.2.21). The Strong 
Maximum Principle implies that they coincide everywhere, which contradicts the 
first line in (5.2.21). 

Hence vr is not a minimizer in H 2 . Let then be a minimizer of £(•, B 2 ), that 
agrees with vr outside the ball B2, and we define the positive quantity 

(5.2.22) SR:=£{vR,B2)-£iv%B2). 

We claim that 

(5.2.23) as R goes to infinity, 6r, remains bounded away from zero. 

To prove this, we assume by contradiction that 

(5.2.24) 

Consider u to be the translation of u, that is u(x) := u{x — e). Let also 

m(x) := min { 71 ( 0 ;), {((x)}. 

We notice that in Bjf /2 we have that u(x) = ur^+(x). This and (5.2.18) give that 

(5.2.25) m = Vr in Bji/2. 

Also, from (5.2.21) and (5.2.25), it follows that m cannot be identically neither u 
nor u, and 

m < u in a neighborhood of 0 and 
m = u in a neighborhood of e. 

Let z be a competitor for m in the ball B2, that agrees with m outside i? 2 . We take 
a cut-off function ip e (^“(IR") such that ip = 1 in 5^/4, ^ = 0 in B^^^. Let 

zr{x) ■■= ip{x)z{x) -f (1 - iP{x))vr{x). 

Then we have that zr = z on Bfij^ and 

(5.2.27) Zr = Vr On 

In addition, by (5.2.25), we have that z = m = vr in Bfi/2\B2- So, it follows that 

Zr{x) = iP{x)vr{x) + {1 - iP{x))vr{x) = Vr{x) = z{x) on Biii 2 \B 2 . 

This and (5.2.27) imply that zr = vr on B^- 

We summarize in the next lines these useful identities (see also Figure 11). 

in B2 


(5.2.26) 


in Bjii2\B2 
in Bjp\B[ii2 
in 


= U, m = VR, Z = Zr 
Ur^+ =U, V% = VR = m = Z = Zr 
= Vr = Zr, 771 = Z 


Ur,+ = U = Vr = Vr = Zr, 


m = z. 
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Figure 11. Energy estimates 


We compute now 


£{m ,B2) - £{z,B2) 

= £{m,B2) -£{vr,B2) + £{vr,B2) -£{zr,B2) + £{zr,B2) -£{z,B2). 


By the definition of Sr in (5.2.22), we have that 
(5.2.28) 

£{m ,B2) - £{z,B2) 

= £{m,B 2 ) - £{vr,B 2 ) + Sr + £{vr,B 2 ) -£{zr,B 2 ) +£{zr,B 2 ) - £{z,B 2 ). 

Using the formula for the kinetic energy given in (5.2.3) together with (5.2.4) we 
have that 

£{m ,^2) - £{vr,B2) 

= -m{B2,B2) + m{B2,B2) + f W[m{x))dx 
2 Jb2 

- \vr{B2,B2) - vr{B2,B2) - f W{vR{x))dx. 

^ J B 2 

Since m = vr on Br^2 (recall (5.2.25)), we obtain 

£{m ,^ 2 ) -£{vr,B 2 ) 

\m{x) - m(y)p - \m{x) - VR{y)\'^ 

Notice now that m and vr are bounded on K” (since so is u). Also, if x e B2 
and y e S ^/2 have that \x — y\ ^ \y\ — |a:| ^ |j/|/2 if R is large. Accordingly, 

£{m,B2) - £{vr,B 2) ^ C { dx { dy ^ 

Jb 2 dB%i \y\ 



(5.2.29) 
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up to renaming constants. Similarly, zu = z on 5^/2 we have the same bound 

(5.2.30) £{zr,B2)-£{z,B2) 

Furthermore, since is a minimizer for £(-,i? 2 ) and wj = zr outside of B 2 , we 
have that 

£{vI,B2)-£{zr,B2)^Q. 

Using this, (5.2.29) and (5.2.30) in (5.2.28), it follows that 
£{m, B 2 ) - £{z, B 2 ) ^ CR-^^ + 5r. 

Therefore, by sending R +00 and using again (5.2.24), we obtain that 

(5.2.31) £{m,B2) ii£{z,B2). 

We recall that z can be any competitor for m, that coincides with m outside of B 2 ■ 
Hence, formula (5.2.31) means that m is a minimizer for £(-,B 2 ). On the other 
hand, u is a minimizer of the energy in any ball. Then, both u and m satisfy the 
same equation in B 2 . Moreover, they coincide in a neighborhood of e, as stated in 
the second line of (5.2.26). By the Strong Maximum Principle, they have to coincide 
on B 2 ^ but this contradicts the first statement of (5.2.26). The proof of (5.2.23) is 
thus complete. 

Now, since = vr on from definition (5.2.22) we have that 
Sr = £{vr,Br) -£{vr,Br). 

Also, £{v’^,Br) > £(u,Br), thanks to the minimizing property of u. Using these 
pieces of information and inequality (5.2.20), it follows that 

Sr ^ £{ur,+,Br) - £{u,Br). 

Now, by sending R +00 and using (5.2.23), we have that 

lim £{ur^+,Br) - £{u,Br) > 0, 

R-^+oo 

which contradicts (5.2.17). This implies that indeed u is monotone, and this con¬ 
cludes the proof of this Step. 

Step 4. Conclusions 

In Step 3, we have proved that u is monotone, in any given direction e. Then, Step 
1 gives the desired result. This concludes the proof of Theorem 5.2.1. □ 

We remark that the exponent two in the energy estimate (5.2.6) is related to 
the expansions of order two and not to the dimension of the space. Indeed, the 
energy estimates hold for any n. However, the two power in the estimate (5.2.6) 
allows us to prove the fractional version of De Giorgi conjecture only in dimension 
two. In other words, the proof of Theorem 5.2.1 is not applicable for n > 2. One 
can verify this by checking the limit in (5.2.17) 

lim (£{ur^+ , Br) - £(u, Br)) = 0, 

R—y+00 \ / 

which was necessary for the Proof of Theorem 5.2.1 in the case n = 2. We know 
from Theorem 5.1.2 that 

= 0. 

R —►+00 R 
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Confronting this result with inequality (5.2.16) 

£{ur.+,Br) -£{u,Br) ^ ^£{u,Br), 

we see that we need to have n = 2 in order for the the limit in (5.2.17) to be zero. 

Of course, the one-dimensional symmetry property in Theorem 5.2.1 is inher¬ 
ited by the spatial homogeneity of the equation, which is translation and rotation 
invariant. In the case, for instance, in which the potential also depends on the 
space variable, the level sets of the (minimal) solutions may curve, in order to 
adapt themselves to the spatial inhomogeneity. 

Nevertheless, in the case of periodic dependence, it is possible to construct 
minimal solutions whose level sets are possibly not planar, but still remain at a 
bounded distance from any fixed hyperplane. As a typical result in this direction, 
we recall the following one (we refer to [45] for further details on the argument): 



Figure 12. Minimal solutions in periodic medium 


Theorem 5.2.3. Let > Q- > 0 and Q : K" —> [Q_,(5_|-]. Suppose 
that Q{x + k) = Q{x) for any k e Z". Let us consider, in any ball Br, the energy 
defined by 

£{u,Br) = ICr{u) + ^ r Q(x)(l-u^fdx, 

^ JBr 

where the kinetic energy K,r(u) is defined as in (5.2.2). 

Then, there exists a constant M > 0, such that, given any ui G dBi, there exists 
a minimal solution u^j of 

{-Ayuuj{x) = Q{x) {u„j{x) - ul,{x)) for any a; G K" 

for which the level sets ^ are contained in the strip {x G M" s.t. jo; • a;| ^ 

M}. 

Moreover, if ui is rotationally dependent, i.e. if there exists kg G Z" such thatui- 
ko = 0, then Uu> is periodic with respect to u, i.e. 

Uuj(x) = Uu){y) for any x, y e K” such that x — y = k and w • fc = 0. 





CHAPTER 6 


Nonlocal minimal surfaces 


In this chapter, we introduce nonlocal minimal surfaces and focus on two main 
results, a Bernstein type result in any dimension and the non-existence of nontrivial 
s-minimal cones in dimension 2. Moreover, some boundary properties will be dis¬ 
cussed at the end of this chapter. For a preliminary introduction to some properties 
of the nonlocal minimal surfaces, see [135], 

Let n cz K" be an open bounded domain, and E c M" be a measurable set, 
fixed outside 17. We will consider for s G (0,1/2) minimizers of the norm 

Notice that only the interactions between E and E^ contribute to the norm. 

In order to define the fractional perimeter of E in 17, we need to clarify the 
contribution of 17 to the iJ® norm here introduced. Namely, as E is fixed outside 
17, we aim at minimizing the “f7-contribution” to the norm among all measurable 
sets that “vary” inside 17. We consider thus interactions between E nfl and E^ and 
between E\[l and f7\E, neglecting the data that is fixed outside 17 and that does 
not contribute to the minimization of the norm (see Figure 13). 



Figure 13. Fractional Perimeter 
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We define the interaction I {A, B) of two disjoint subsets of M" as 

dx dy 


(6.0.32) 


i{A,B)-.= { r 

Ja Ji 

J 


AJB F- y|"+2^ 

Xa{x)xb{x) 


dx dy. 


\X - 

Then (see [29]), one defines the nonlocal s-perimeter functional of E in as 
(6.0.33) Per^iE, Pi) := I{E n fl, E^) + I{E\Vt, n\E). 

Equivalently, one may write 

Per^(£;, f2) = I{E n VL, n\E) + I{E n PI, E^\n) + I{E\n, n\E). 


Definition 6.0.4. Let PI be an open domain o/K". A measurable set E c K” 
is s-minimal in PI if Per^{E, Pi) is finite and if for any measurable set E such that 
E\Pl = E\Pl, we have that 

Pers{E,Pl) ^ Pera{E,Pi). 

A measurable set is s-minimal in M" if it is s-minimal in any ball Br, where r > 0. 

When s ^ the fractional perimeter Perg approaches the classical perimeter, 
see [22]. See also [46] for the precise limit in the class of functions with bounded 
variations, [35, 36] for a geometric approach towards regularity and [117, 6[ for an 
approach based on P-convergence. See also [136] for a different proof and Theorem 
2.22 in [104] and the references therein for related discussions. A simple, formal 
statement, up to renormalizing constants, is the following: 


Theorem 6.0.5. Let R> Q and E be a set with finite perimeter in Br. Then 
linj - s'jPera{E,Br) = Per {E,Br) 

for almost any r e (0, R). 


The behavior of Per^ as s —> 0 is slightly more involved. In principle, the limit 
as s ^ 0 of Perg is, at least locally, related to the Lebesgue measure (see e.g. [107]). 
Nevertheless, the situation is complicated by the terms coming from infinity, which, 
as s —> 0, become of greater and greater importance. More precisely, it is proved 
in [59] that, if Perg^(if, PI) is finite for some So e (0,1/2), and the limit 

( 6 . 0 . 34 ) & 1 „ 2 . 

exists, then 


(6.0.35) lim2sPerg(£;,D) = (jdBij -/?£;) \EnPl\+l3E jD/Ej. 

s —>0 

We remark that, using polar coordinates, 

0 ^/3b ^ lirn 2s r —= lirn 2s |5Bi| f dp = jdSi|, 

therefore /3 b e [0, |(3i3i|] plays the role of a convex interpolation parameter in the 
right hand-side of (6.0.35) (up to normalization constants). 

In this sense, formula (6.0.35) may be interpreted by saying that, as s —> 0, the 
s-perimeter concentrates itself on two terms that are “localized” in the domain PI, 
namely [if n flj and |n\i3|. Nevertheless, the proportion in which these two terms 
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count is given by a “strongly nonlocal” interpolation parameter, namely the quan¬ 
tity Pe in (6.0.34) which “keeps track” of the behavior of E at infinity. 

As a matter of fact, to see how Pe is influenced by the behavior of E at infinity, 

one can compute Pe for the particular case of a cone. For instance, if S c dBi, 
ISI 

with =: be [0,1], and E is the cone over S (that is E := {tp, pe t ^ 0}), 
we have that 

^ + 00 

/3E = lim2s|S| p-^-^Up=\E\ = b\dBil 

that is Pe gives in this case exactly the opening of the cone. 

We also remark that, in general, the limit in (6.0.34) may not exist, even for 
smooth sets: indeed, it is possible that the set E “oscillates” wildly at infinity, say 
from one cone to another one, leading to the non-existence of the limit in (6.0.34). 

Moreover, we point out that the existence of the limit in (6.0.34) is equivalent to 
the existence of the limit in (6.0.35), except in the very special case |if nSlI = |n\if|, 
in which the limit in (6.0.35) always exists. That is, the following alternative holds 
true: 

• if |A n n| 7 ^ then the limit in (6.0.34) exists if and only if the limit 

in (6.0.35) exists, 

• if jA n n| = |f2\£’|, then the limit in (6.0.35) always exists (even when the 
one in (6.0.34) does not exist), and 

lim2sPers(F;,n) = \dBi\ jA n f2| = \dBi\ |f2\F;|. 

s —>0 

The boundaries of s-minimal sets are referred to as nonlocal minimal surfaces. 
In [29] it is proved that s-minimizers satisfy a suitable integral equation (see in 
particular Theorem 5.1 in [29]), that is the Euler-Lagrange equation corresponding 
to the s-perimeter functional Perg. If E is s-minimal in and dE is smooth enough, 
this Euler-Lagrange equation can be written as 

( 6 . 0 . 36 ) 

for any xq G n dE. 

Therefore, in analogy with the case of the classical minimal surfaces, which have 
zero mean curvature, one defines the nonlocal mean curvature of E ai, xq e dE as 

(n n 0^3 TTS ^ \ r j 

( 6 . 0 . 37 ) 

In this way, equation (6.0.36) can be written as Elf, = 0 along dE. 

It is also suggestive to think that the function xe '■= Xe — Xe^ averages out 
to zero at the points on dE, if dE is smooth enough, since at these points the local 
contribution of E compensates the one of E^. Using this notation, one may take 
the liberty of writing 

1 r Xe{xo + y)+ Xe{xo -y) , 

2 V ^ 

If XE{xo + y)+XE{xo-y)-2xE{xo) 

2 J...- 

-(-A)'^Xi;((ro) 

C{n,s) ’ 


He{xo) 
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using the notation of (2.0.1). Using this suggestive representation, the Euler- 
Lagrange equation in (6.0.36) becomes 

(—A)®X£; = 0 along dE. 

We refer to [3] for further details on this argument. 

It is also worth recalling that the nonlocal perimeter functionals find appli¬ 
cations in motions of fronts by nonlocal mean curvature (see e.g. [33, 92, 41]), 
problems in which aggregating and disaggregating terms compete towards an equi¬ 
librium (see e.g. [79] and [57]) and nonlocal free boundary problems (see e.g. [30] 
and [62]). See also [107] and [139] for results related to this type of problems. 

In the classical case of the local perimeter functional, it is known that minimal 
surfaces are smooth in dimension n ^ 7. Moreover, if n > 8 minimal surfaces are 
smooth except on a small singular set of Hausdorff dimension n — 8. Furthermore, 
minimal surfaces that are graphs are called minimal graphs, and they reduce to 
hyperplanes if n ^ 8 (this is called the Bernstein property, which was also discussed 
at the beginning of the Chapter 5). If n > 9, there exist global minimal graphs that 
are not affine (see e.g. [87]). 

Differently from the classical case, the regularity theory for s-minimizers is still 
quite open. We present here some of the partial results obtained in this direction: 

Theorem 6.0.6. In the plane, s-minimal sets are smooth. More precisely: 

a) If E is an s-minimal set in Q cz then dE n is a -curve. 

b) Let E be s-minimal in LI cz K” and let c dE n LI denote its singular set. 
Then 'H'^{Ye) = 0 for any d > n — 8. 

See [124] for the proof of this results (as a matter of fact, in [124] only C^’“ 
regularity is proved, but then [9] proved that s-minimal sets with C^’“-boundary 
are automatically C®). Further regularity results of the s-minimal surfaces can be 

found in [36] . There, a regularity theory when s is near - is stated, as we see in 
the following Theorem: 

Theorem 6.0.7. There exists eg e ^0, such that if s ^ - — eg, then 

a) if n ^ 7, any s-minimal set is of class C®, 

b) if n = 8 any s-minimal surface is of class (7® except, at most, at countably many 
isolated points, 

c) any s-minimal surface is of class C” outside a closed set E of Hausdorff dimen¬ 
sion n — 8. 


6.1. Graphs and s-minimal surfaces 

We will focus the upcoming material on two interesting results related to graphs: 
a Bernstein type result, namely the property that an s-minimal graph in is 

flat (if no singular cones exist in dimension n) ; we will then prove that an s-minimal 
surface whose prescribed data is a subgraph, is itself a subgraph. 

The first result is the following theorem: 

Theorem 6.1.1. Let E = {{x,t) e M” x K s.f. t < u{x)} be an s-minimal graph, 
and assume there are no singular cones in dimension n (that is, if K, ez M" is an 
s-minimal cone, then 1C is a half-space). Then u is an affine function (thus E is a 
half-space). 
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To be able to prove Theorem 6.1.1, we recall some useful auxiliary results. In 
the following lemma we state a dimensional reduction result (see Theorem 10.1 in 
[29]). 

Lemma 6.1.2. Let E = F x W. Then if E is s-minimal if and only if F is 
s-minimal. 

We define then the blow-up and blow-down of the set E are, respectively 

E 

Eq := Ym\ Ej. and Eao '■= lim Er, where E^- = —. 

1 —>0 r^+00 r 

A first property of the blow-up of E is the following (see Lemma 3.1 in [80]). 
Lemma 6.1.3. If Eqq is affine, then so is E. 

We recall also a regularity result for the s-minimal surfaces (see [80] and [9] for 
details and proof). 

Lemma 6.1.4. Let E he s-minimal. Then: 

a) If E is Lipschitz, then E is C^’“. 

b) If E is then E is . 

We give here a sketch of the proof of Theorem 6.1.1 (see [80] for all the details). 

Sketch of the proof of Theorem 6.1.1. If E c is an s-minimal 

graph, then the blow-down E^o is an s-minimal cone (see Theorem 9.2 in [29] 
for the proof of this statement). By applying the dimensional reduction argument 
in Lemma 6.1.2 we obtain an s-minimal cone in dimension n. According to the 
assumption that no singular s-minimal cones exist in dimension n, it follows that 
necessarily Eqo can be singular only at the origin. 

We consider a bump function wq e C®(M, [0,1]) such that 

wo{t) = 0 in ^ - 00, u -too 

wo(t) = 1 in 
w{t) = Wod^D- 

The blow-down of E is 

Aqo = {{x',Xn+i) S.t. Xn+I < Uoo(a;')}- 
For a fixed cr e dBi , let 

Ft := {{x',Xn+i) s.t. x„+i ^ Uoo{x' -Y t0w{x')a) - t} 
be a family of sets, where t e (0,1) and 9 > 0. Then for 6 small, we have that 

(6.1.1) El is below Eao- 

Indeed, suppose by contradiction that this is not true. Then, there exists 9k ^ 0 
such that 

(6.1.2) Woo(a;fc + 9kw{x'k)(j) - 1 > Uaoix'k). 

But xf e suppru, which is compact, therefore x'^ := lim x'k belongs to the support 

fe—>. + 00 

of w, and w{x',„) is defined. Then, by sending k —> -too in (6.1.2) we have that 

Woo(a;'oo) - 1 > ^^(a^'oo)^ 
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which is a contradiction. This establishes (6.1.1). 

Now consider the smallest to e (0,1) for which Ft is below Ea^- Since Ea^ 
is a graph, then touches E^ from below in one point Xq = (xg, where 

Xq e suppw. Now, since Fqq is s-minimal, we have that the nonlocal mean curvature 
(defined in (6.0.37)) of the boundary is null. Also, since Etg is a diffeomorphism 
of Fqo we have that 

(6.1.3) 

and there is a region where Eqq and Etg are well separated by tg, thus 

|(Foo\Ft„) n {B3\B2) \ > cto, 
for some c > 0. Therefore, we see that 

^kg (P) = ^kg (P) - ^sip) > etg. 

This and (6.1.3) give that 9to > cto, for some c > 0 (up to renaming it). If 9 is 
small enough, this implies that tg = 0. 

In particular, we have proved that there exists 9 > 0 small enough such that, 
for any t e (0,1) and any a e dBi, we have that 

Uoo(x' + t9w(x')a) — t ^ Uoo(x'). 

This implies that 

Uoo(x' + t9w(x')a) — Uoo(x') 1 

r9 ^ 0’ 

hence, letting t ^ 0, we have that 

\/uoo{x')w{x')a ^ for any x e IR"'\{0}, and a e Bi. 

U 

We recall now that w = 1 in F 3 / 5 \i? 2/5 and a is arbitrary in dBi. Hence, it follows 
that 

|VMoo(a;)| ^ for any x e B3/3\B2/5- 

Therefore Mqo is globally Lipschitz. By the regularity statement in Lemma 6.1.4, 
we have that Uao is C®. This says that u is smooth also at the origin, hence (being 
a cone) it follows that E^g is necessarily a half-space. Then by Lemma 6.1.3, we 
conclude that E is a half-space as well. □ 

We introduce in the following theorem another interesting property related to 
s-minimal surfaces, in the case in which the fixed given data outside a domain is a 
subgraph. In that case, the s-minimal surface itself is a subgraph. Indeed: 

Theorem 6.1.5. Let Og be an open and bounded subset with boundary 

of class and let H := Og x K. Let E be an s-minimal set in fl. Assume that 

(6.1.4) E\Ll = {xn < u{x'), x' e IR"“^\ng} 

for some continuous function u: —> K. Then 

E n LI = {xn < v{x'), x' e fig} 
for some function v: ^ IR. 
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The reader can see [65], where this theorem and the related results are proved; 
here, we only state the preliminary results needed for our purposes and focus on 
the proof of Theorem 6.1.5. The proof relies on a sliding method, more precisely, 
we take a translation of E in the direction, and move it until it touches E from 
above. If the set £1 n is a subgraph, then, up to a set of measure 0, the contact 
between the translated E and E^ will be E itself. 

However, since we have no information on the regularity of the minimal sur¬ 
face, we need at first to “regularize” the set by introducing the notions of sup- 
convolution and subconvolution. With the aid of a useful result related to the 
sub/supconvolution of an s-minimal surface, we proceed then with the proof of the 
Theorem 6.1.5. 

The supconvolution of a set E c K"’ is given by 

4 := U M^). 

xeE 



In an equivalent way, the supconvolution can be written as 

£;“= IJ (E + v). 

veR'^ 


Indeed, we consider <5 > 0 and an arbitrary x e E. Let y e Bs{x) and we define 
V := y — X. Then 

|u| ^ |j/ — a;| ^ (5 and y = x + veE + v. 

Therefore Bs{x) ^ E + v for |u| ^ <5. In order to prove the inclusion in the opposite 
direction, one notices that taking y e E + v with |u| ^ S and defining x := y — v, it 
follows that 

\x-y\ = |u| ^ S. 

Moreover, xe{E + v)—v = E and the inclusion E + v e Bs{x) is proved. 

On the other hand, the subconvolution is defined as 

e'^s := M”\ (^{R^\E)l'^ . 
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Now, the idea is that the supconvolution of i? is a regularized version of E whose 
nonlocal minimal curvature is smaller than the one of E, i.e.: 


(6.1.5) 




for any x e dE^, where x := x — v e dE for some v e K" with |n| = 6. Then, by 
construction, the set E + v lies in Eg, and this implies (6.1.5). 

Similarly, one has that the opposite inequality holds for the subconvolution 
of E, namely 


( 6 . 1 . 6 ) 



(y) - (y) 

|a; - y|"+2s 


dy > 0, 


By (6.1.5) and (6.1.6), we obtain: 

Proposition 6.1.6. Let E be an s-minimal set in Ll. Let p G dEg and assume 
that Bs{p) £ LI. Assume also that Eg is touched from above by a translation of Eg, 
namely there exists lo e K" such that 

eI<^E\ + u: 


and 


Then 


pe{dEl)n{dE\+oj). 
e\ = e\ + cc. 


Proof of Theorem 6.1.5. One first remark is that the s-minimal set does 
not have spikes which go to infinity: more precisely, one shows that 

(6.1.7) Llo X (-CX), -M) c £; n O c flo X {-oo, M) 


for some M > 0. The proof of (6.1.7) can be performed by sliding horizontally a 
large ball, see [65] for details. 

After proving (6.1.7), one can deal with the core of the proof of Theorem 6.1.5. 
The idea is to slide E from above until it touches itself and analyze what happens 
at the contact points. For simplicity, we will assume here that the function u is 
uniformly continuous (if u is only continuous, the proof needs to be slightly modified 
since the subconvolution and supconvolution that we will perform may create new 
touching points at infinity). 

At this purpose, we consider Et = E + ten for t ^ 0. Notice that, by (6.1.7), if 
t > 2M, then E Et- Let then t be the smallest for which the inclusion E Et 
holds. We claim that t = 0. If this happens, one may consider 

V = inf{T s.t. {x,t) e E^} 


and, up to sets of measure 0, A n Oq is the subgraph of v. 

The proof is by contradiction, so let us assume that t > 0. According to (6.1.4), 
the set E\fl is a subgraph, hence the contact points between dE and dEt must lie 
in rio X M. Namely, only two possibilities may occur: the contact point is interior 
(it belongs to Qq x K), or it is at the boundary (on dflg ^ ®)- So, calling p the 
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contact point, one may have^ that 

(6.1.8) either p e Qq x R or 

(6.1.9) pednoxR. 

We deal with the first case in (6.1.8) (an example of this behavior is depicted 
in Figure 15). We consider Eg and Eg to be the supconvolution, respectively the 
subconvolution of E. We then slide the subconvolution until it touches the supcon¬ 
volution. More precisely, let r > 0 and we take a translation of the subconvolution. 
Eg + Ten- For r large, we have that e\ ^ E\ + ren and we consider Tg to be the 
smallest for which such inclusion holds. We have (since t is positive by assumption) 
that 

> ( > 0 . 

Moreover, for 5 small, the sets dEg and d^E'g + Tge„) have a contact point which, 
according to (6.1.8), lies in fig x K. Let pg be such a point, so we may write 

pg e (dE^g) n d(E^g + TgCn) and pg e fto x K. 

Then, for S small (notice that Bg{p) c fl), Proposition 6.1.6 yields that 

Eg = Eg TgCn- 

Considering S arbitrarily small, one obtains that 

E = E + Toe„, with tq > 0. 

But if is a subgraph outside of 12, and this provides a contradiction. Hence, the 
claim that i = 0 is proved. 


Let us see that we also obtain a contradiction when supposing that t > 0 and that 
the second case (6.1.9) holds. Let 

p = {p,Pn) and p e {dE) n {dEt). 

Now, if one takes sequences G dE and bk G dEt, both that tend to p as k goes to 
infinity, since E\[l is a subgraph and t > 0, necessarily ak,bk belong to 12. Hence 

(6.1.10) p G (dE) n 12 n (dEt) n 12. 

Thanks to Definition 2.3 in [29], one obtains that if is a variational subsolution in 
a neighborhood of p. In other words, if H c if n 12 and p e A, then 

0 > Per,(F;, 12) - Per,(F;\H, 12) = I {A, E‘^) - I {A, E\A) 


^As a matter of fact, the number of contact points may be higher than one, and even infinitely 
many contact points may arise. So, to be rigorous, one should distinguish the case in which all the 
contact points are interior and the case in which at least one contact point lies on the boundary. 

Moreover, since the surface may have vertical portions along the boundary of the domain, 
one needs to carefully define the notion of contact points (roughly speaking, one needs to take a 
definition for which the vertical portions which do not prevent the sliding are not in the contact 
set). 

Finally, in case the contact points are all interior, it is also useful to perform the sliding 
method in a slighltly reduced domain, in order to avoid that the supconvolution method produces 
new contact points at the boundary (which may arise from vertical portions of the surfaces). 

Since we do not aim to give a complete proof of Theorem 6.1.5 here, but just to give the main 
ideas and underline the additional difficulty, we refer to [ 65 ] for the full details of these arguments. 
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Figure 15. Sliding E until it touches itself at an interior point 


(we recall the definition of I in (6.0.32) and of the fractional perimeter PeVs in 
(6.0.33)). According to Theorem 5.1 in [29], this implies in a viscosity sense (i.e. if 
E is touched at p from outside by a ball), that 


( 6 . 1 . 11 ) 



XEjy) - Xr"\e(2/) 
\p - y|"+2s 


dy 5= 0. 


In order to obtain an estimate on the fractional mean curvature in the strong sense, 
we consider the translation of the point p as follows: 


Pt= p-te„ = {p',Pn -t) = {p\pn,t)- 

Since t > 0, one may have that either ^ u{p'), or pn,t u{p'). 

These two possibilities can be dealt with in a similar way, so we just continue 
with the proof in the case A u{p') (as is also exemplified in Figure 16). Taking 
r > 0 small, the set Br{p)\Xl is contained entirely in E or in its complement. 
Moreover, one has from [28] that dE n Br{p) is a 2 +®-graph in the direction 
of the normal to fl at p. That is: in Figure 16 the set E is C'^’ 2 +®, hence in the 
vicinity oi p = {p',Pn), it appears to be sufficiently smooth. 

So, let v{p) = (u'(p), u„(p)) be the normal in the interior direction, then up 
to a rotation and since fl is a cylinder (hence Vn{p) = 0), we can write u(p) = ei. 
Therefore, there exists a function ik of class C'^’ 2 +® such that pi = 'I'(p 2 ; ■ ■ ■ ,Pn) 
and, in the vicinity of p, we can write dE as the graph G = {xi = ^'( 2 : 2 ,..., Xn)}- 
Given (6.1.10), we deduce that there exists a sequence pk E G such that pk e XI 
and pk ^ p as k ^ cx). From this it follows that there exists a sequence of points 
Pk ^ p such that 

(6.1.12) dE in the vicinity of pk is a graph of class 


and 

(6.1.13) 



XE{y) - XEc{y) 


dy = 0. 


\pk - 2/|"+2" 
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Figure 16. Sliding E until it touches itself at a boundary point 


From (6.1.12) and (6.1.13), and using a pointwise version of the Euler-Lagrange 
equation (see [65] for details), we have that 

xsiy) - XEciy) 


I 


\p - y|’T’+2s 

Now, E Et for t strictly positive, hence 

I XEM-XE9iy) 


j 

JR 


dy = 0. 


dy > 0. 


\p - y|"+2s 

Moreover, we have that the set dEt n Br{p) must remain on one side of the graph 
G, namely one could have that 

Et n Br(p) C {xi ^ .. -.Xn)} Or 

Et n Br[p) 3 > 'i>{x2, ■ ■ ■,Xn)}. 

Given again (6.1.10), we deduce that there exists a sequence pk e dEt n 12 such that 
Pfc P as fc —> 00 and dEt n 12 in the vicinity of pk is touched by a surface lying in 
Et, of class G^’5+®. Then 

XEtiy) -XE^iy) 


i 


■dy ^ 0. 


\Pk - 

Hence, making use of a pointwise version of the Euler-Lagrange equation (see [65] 
for details), we obtain that 

XEtiy) -XE^iy) 


i 


■dy ^ 0. 


\p - y|"+2'< 

But this is a contradiction with (6.1.14), and this concludes the proof of Theo¬ 
rem 6.1.5. □ 


On the one hand, one may think that Theorem 6.1.5 has to be well-expected. 
On the other hand, it is far from being obvious, not only because the proof is not 
trivial, but also because the statement itself almost risks to be false, especially at 
the boundary. Indeed we will see in Theorem 6.3.2 that the graph property is close 
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to fail at the boundary of the domain, where the s-minimal surfaces may present 
vertical tangencies and stickiness phenomena (see Figure 23). 

6.2. Non-existence of singular cones in dimension 2 

We now prove the non-existence of singular s-minimal cones in dimension 2, as 
stated in the next result (from this, the more general statement in Theorem 6.0.6 
follows after a blow-up procedure): 

Theorem 6.2.1. If E is an s-minimal cone in then E is a half-plane. 

We remark that, as a combination of Theorems 6.1.1 and 6.2.1, we obtain the 
following result of Bernstein type: 

COROLARY 6.2.2. Let E = {{x,t) e K” x M s.t. t < u(ai)} be an s-minimal 
graph, and assume that ns {1,2}. Then u is an affine function. 



Figure 17. The cone K. 


Let us first consider a simple example, given by the cone in the plane 

K, := |(x, y) G s.t. i/^ > 

see Figure 17. 

Proposition 6.2.3. The cone 1C depicted in Figure 17 is not s-minimal in K^. 

Notice that, by symmetry, one can prove that JC satisfies (6.0.36) (possibly 
in the viscosity sense). On the other hand. Proposition 6.2.3 gives that K, is not 
s-minimal. This, in particular, provides an example of a set that satisfies the 
Euler-Lagrange equation in (6.0.36), but is not s-minimal (i.e., the Euler-Lagrange 
equation in (6.0.36) is implied by, but not necessarily equivalent to, the s-minimality 
property). 

Proof of Proposition 6.2.3. The proof of the non-minimality of JC is due 
to an original idea by Luis Caffarelli. 
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Suppose by contradiction that the cone JC is minimal in . We add to /C a small 
square adjacent to the origin (see Figure 18), and call 1C' the set obtained. Then 
/C and JC' have the same s-perimeter. This is due to the interactions considered in 
the s-perimeter functional and the unboundedness of the regions. We remark that 
in Figure 18 we represent bounded regions, of course, sets A, S, C, D, A', B', C and 
D' are actually unbounded. 



Figure 18. Interaction of M with A, B, C, B, A'^ B', C, D' 


Indeed, we notice that in the first image, the white square M interacts with the 
dark regions A, B,C, D, while in the second the now dark square M interacts with 
the regions A',B',C',D', and all the other interactions are unmodified. There¬ 
fore, the difference between the s-perimeter of 1C and that of 1C' consists only of 
the interactions I{A, M) + I{B, M) + I{C, M) + I{D, M) — I{A', M) — I{B', M) — 
I{C',M) - I{D',M). But Au B = A' yj B' and C u D = C u D' (since these 
sets are all unbounded), therefore the difference is null, and the s-perimeter of K. is 
equal to that of 1C'. Consequently, 1C' is also s-minimal, and therefore it satisfies the 
Euler-Lagrange equation in (6.0.36) at the origin. But this leads to a contradiction, 
since the the dark region now contributes more than the white one, namely 



XJC'jy) - XR'^wiy) 

\y\2+S 


dy > 0. 


Thus 1C cannot be s-minimal, and this concludes our proof. 


□ 


This geometric argument cannot be extended to a more general case (even, 
for instance, to a cone in made of many sectors, see Figure 19). As a matter 
of fact, the proof of Theorem 6.2.1 will be completely different than the one of 
Proposition 6.2.3 and it will rely on an appropriate domain perturbation argument. 

The proof of Theorem 6.2.1 that we present here is actually different than 
the original one in [124]. Indeed, in [124], the result was proved by using the 
harmonic extension for the fractional Laplacian. Here, the extension will not be 
used; furthermore, the proof follows the steps of Theorem 5.2.1 and we will recall 
here just the main ingredients. 
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Figure 19. Cone in 

Proof of Theorem 6.2.1. The idea of the proof is the following: if if c 
is an s-minimal cone, then let if be a perturbation of the set E which coincides 
with a translation of E in 5^/2 and with E itself outside Bn. Then the difference 
between the energies of E and E tends to 0 as i? ^ +oo. This implies that also the 
energy of if n if is arbitrarily close to the energy of E. On the other hand if E is 
not a half-plane, the set E E can be modified locally to decrease its energy by a 
fixed small amount and we reach a contradiction. 

The details of the proof go as follows. Let 

u '■= Xe — Xr 2 \e- 

From definition (5.2.4) we have that 

u{Bn, Bn) = 2i(if n if^, Bn\E) 

and 

u{Bn, Bi) = I {Bn n F, E^\Bn) + I{Bn\E, E\Bn), 

thus 

(6.2.1) Pers{E,Bn) = JCn{u), 

where JCn{u) is given in (5.2.2) and Pers(if, Bn) is the s-perimeter functional defined 
in (6.0.33). Then E is s-minimal if u is a minimizer of the energy ICn in any ball 
Bn, with R > 0. 

Now, we argue by contradiction, and suppose that E is an s-minimal cone 
different from the half-space. Up to rotations, we may suppose that a sector of E 
has an angle smaller than tt and is bisected by 62 . Thus there exists M ^ 1 and 
p e E n Bm on the e 2 -axis such that p ± ei e M.'^\E (see Figure 19). 

We take p e C^{Bi), such that (p{x) = 1 in i?i/ 2 . For R large (say R > 8M), 
we define 

'pR..+{y) ■=y + ‘^(d) ei- 

We point out that, for R large, 41^; + is a diffeomorphism on 
Furthermore, we define u]^{x) := rt(T)^^_|_(a:)). Then 

Un{y) = u{y - ei) for p e B 2 M 

and u^{y) = u{y) for p e ^\Bn- 
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We recall the estimate obtained in (5.2.6), that, combined with the minimality of 
M, gives 

But u is a minimizer in any ball, and by the energy estimate in Theorem 5.1.2 we 
have that 

1Cr{u+^) - JCr{u) ^ CR-^^. 

This implies that 

(6.2.2) ^lim^/Ci^(u+) - 1Cr{u) = 0. 

Let now 

Ufl(a;) := min{u(a:),Mj(a;)} and ^^^(a;) := max{u(a:),Mj(a;)}. 

We claim that vr is not identically u nor Indeed 

= Hp - ei) = iXE - Xh 2 \£;)(p - ei) = -1 and 

u{p) = iXE - Xr^\e){p) = 1- 

On the other hand, 

'^r(p + ®i) “ 1 

u{p+ei) = ixE - Xm\E){p + ei) = -1. 

By the continuity of u and uj, we obtain that 

(6.2.3) Vr = u'^ < u in a, neighborhood of p 
and 

(6.2.4) Vr = u < rtj in a neighborhood of p -I- ei. 

Now, by the minimality property of u, 

ICr{u) ^ ICr{vr). 

Moreover (see e.g. formula (38) in [114]), 

^r(l’r) + K.r{wr) ^ JCr{u) + K,r{u^). 

The latter two formulas give that 

(6.2.5) JCr{vr) ^ JCr{uI). 

We claim that 

(6.2.6) Vr, is not minimal for IC 2 M 

with respect to compact perturbations in i? 2 M- Indeed, assume by contradiction 
that Vr is minimal, then in i?2M both vr and u would satisfy the same equation. 
Recalling (6.2.4) and applying the Strong Maximum Principle, it follows that u = vr 
in B 2 M, which contradicts (6.2.3). This establishes (6.2.6). 

Now, we consider a minimizer of K, 2 M among the competitors that agree 
with Vr outside i? 2 M- Therefore, we can define 

5r := 1C2m{vr) — K^2m{ur)- 

In light of (6.2.6), we have that 6r > 0. 

The reader can now compare Step 3 in the proof of Theorem 5.2.1. There we 
proved that 

(6.2.7) Sr remains bounded away from zero as R ^ +00. 
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Furthermore, since and vr agree outside B 2 M we obtain that 

= ^r{vr). 

Using this, (6.2.5) and the minimality of u, we obtain that 

Sr = ICr{vr) - JCr{u*r) ^ ICr{u+) - 1Cr{u). 

Now we send R to infinity, recall (6.2.2) and (6.2.7), and we reach a contradiction. 
Thus, E is, a, half-space, and this concludes the proof of Theorem 6.2.1. □ 

As already mentioned, the regularity theory for s-minimal sets is still widely 
open. Little is known beyond Theorems 6.0.6 and 6.0.7, so it would be very interest¬ 
ing to further investigate the regularity of s-minimal surfaces in higher dimension 
and for small s. 


It is also interesting to recall that if the s-minimal surface E is a subgraph of 
some function u : ^ M (at least in the vicinity of some point xq = (xq, ^(xq)) g 

dE) then the Euler-Lagrange (6.0.36) can be written directly in terms of u. For 
instance (see formulas (49) and (50) in [9]), under appropriate smoothness assump¬ 
tions on u, formula (6.0.36) reduces to 


0 = 


J 

i 


-\e{xo + y)- Xe{xo + y) 

|y|n+2s 

u{.x'q + y') - u(x'o) 


\y'\ 


dy 

ay') 


\y'\ 


n—H-2s 


dy' + T(a:o), 


for suitable F and ik, and a cut-off function C supported in a neighborhood of x'q. 



Figure 20. A nonlocal catenoid 


Regarding the regularity problems of the s-minimal surfaces, let us mention the 
recent papers [49] and [50]. Among other very interesting results, it is proved there 
that suitable singular cones of symmetric type are unstable up to dimension 6 but 
become stable in dimension 7 for small s (these cones can be seen as the nonlocal 
analogue of the Lawson cones in the classical minimal surface theory, and the stabil¬ 
ity property is in principle weaker than minimality, since it deals with the positivity 
of the second order derivative of the functional). 
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This phenomenon may suggest the conjecture that the s-minimal surfaces may 
develop singularities in dimension 7 and higher when s is sufficiently small. 

In [ 50 ], interesting examples of surfaces with vanishing nonlocal mean curvature 
are provided for s sufficiently close to 1/2. Remarkably, the surfaces in [ 50 ] are 
the nonlocal analogues of the catenoids, but, differently from the classical case (in 
which catenoids grow logarithmically), they approach a singular cone at infinity, see 
Figure 20. 

Also, these nonlocal catenoids are highly unstable from the variational point of 
view, since they possess infinite Morse index (differently from the standard catenoid, 
which has Morse index equal to one, i.e. it is, roughly speaking, a minimizer in any 
functional direction with the exception of one). 

Moreover, in [ 50 ], there are also examples of surfaces with vanishing nonlocal 
mean curvature that can be seen as the nonlocal analogues of two parallel hyper¬ 
planes. Namely, for s sufficiently close to 1/2, there exists a surface of revolution 
made of two sheets which are the graph of a radial function / = ±/(r). When r 
is small, / is of the order of 1 -f (^ — s)r^, but for large r it becomes of the order 

of \ — s-r. That is, the two sheets “repel each other” and produce a linear growth 
at infinity. When s approaches 1/2 the two sheets are locally closer and closer to 
two parallel hyperplanes, see Figure 21. 

The construction above may be extended to build families of surfaces with 
vanishing nonlocal mean curvature that can be seen as the nonlocal analogue of 
k parallel hyperplanes, for any fc e N. These fc-sheet surfaces can be seen as the 
bifurcation, as s is close to 1/2, of the parallel hyperplanes {xn = Oi}, for i e 
{1,..., /c}, where the parameters Ui satisfy the constraints 

k 

(6.2.8) ai > • • • > afc, ^ Qi = 0 

2 = 1 

and the balancing relation 

(6.2.9) a, = 2 V - . 

n .: — n ^ 


It is actually quite interesting to observe that solutions of (6.2.9) correspond to 
(nondegenerate) critical points of the functional 

1 ^ 

E{ai,...,au) ■■= ^ (-l)*+Mog |a, - | 

i=l KSjSr. 

among all the fc-ples (oi,... ,afc) that satisfy (6.2.8). 

These bifurcation techniques rely on a careful expansion of the fractional perime¬ 
ter functional with respect to normal perturbations. That is, if if is a (smooth) set 
with vanishing fractional mean curvature, and h is a smooth and compactly sup¬ 
ported perturbation, one can define, for any t e M, 

Eh{t) := {x + th{x)i>{x), X e dE], 
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r 



Figure 21. A two-sheet surface with vanishing fractional 
mean curvature 


where i'{x) is the exterior normal of E at x. Then, the second variation of the 
perimeter of Eh(t) at t = 0 is (up to normalization constants) 


h{y) - h{x) 


■ |n+2s 


dE \x-y\ 

h{y) - h{x) 


\x - y\ 


n+2s 






-^{y) + h{x) r 

JdE 


{e{x) - v{y)) ■ v{x) 


\x - y|"+2s 
1 - v{x) ■ v{y ) 

SE |ai-y|-+2« 


dn^-\y) 


Jdj 


Notice that the latter integral is non-negative, since ^{x) ■ v{y) ^ 1. The quantity 
above, in dependence of the perturbation h, is called, in jargon, “Jacobi operator”. It 
encodes an important geometric information, and indeed, as s ^ 1/2, it approaches 
the classical operator 

^dEh + h, 

where AgE is the Laplace-Beltrami operator along the hypersurface dE and |Ag£:p 
is the sum of the squares of the principal curvatures. Other interesting sets that 


possess constant nonlocal mean curvature with the structure of onduloids have 
been recently constructed in [25] and [48]. This type of sets are periodic in a given 
direction and their construction has perturbative nature (indeed, the sets are close 
to a slab in the plane). 

It is interesting to remark that the planar objects constructed in [25] have no 
counterpart in the local framework, since hypersurfaces of constant classical mean 
curvature with an onduloidal structure only exist in K" with n > 3: once again, 
this is a typical nonlocal effect, in which the nonlocal mean curvature at a point is 
influenced by the global shape of the set. 


While unbounded sets with constant nonlocal mean curvature and interest¬ 
ing geometric features have been constructed in [50, 25], the case of smooth and 
bounded sets is always geometrically trivial. As a matter of fact, it has been recently 
proved independently in [25] and [44] that bounded sets with smooth boundary and 
constant mean curvature are necessarily balls (this is the analogue of a celebrated 
result by Alexandrov for surfaces of constant classical mean curvature). 
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6.3. Boundary regularity 

The boundary regularity of the nonlocal minimal surfaces is also a very inter¬ 
esting, and surprising, topic. Indeed, differently from the classical case, nonlocal 
minimal surfaces do not always attain boundary data in a continuous way (not 
even in low dimension). A possible boundary behavior is, on the contrary, a com¬ 
bination of stickiness to the boundary and smooth separation from the adjacent 
portions. Namely, the nonlocal minimal surfaces may have a portion that sticks at 
the boundary and that separates from it in a C'^’ 2 +'*-way. 



Figure 22. Stickiness properties of 
Theorem 6.3.1. 


As an example, we can consider, for any 5 > 0, the spherical cap 

Ks ■■= {Bi+s\Bi) n {xn < 0}, 

and obtain the following stickiness result: 

Theorem 6.3.1. There exists <5o > 0, depending on n and s, such that for 
any 6 e (0,(5o], we have that the s-minimal set in Bi that coincides with Ks out¬ 
side Bi is Ks itself. 

That is, the s-minimal set with datum Ks outside Bi is empty inside Bi. 

The stickiness property of Theorem 6.3.1 is depicted in Figure 22. 

Other stickiness examples occur at the sides of slabs in the plane. For instance, 
given M > 1, one can consider the s-minimal set Em in (—1, 1 ) x K with datum 
outside ( — 1, 1) X R given by the “jump” set Jm '■= Jm ^ Jm^ where 

Jm ■= (-00:-1] X (-00,-M) 
and := [1,-l-cx)) x (—cx),M). 

Then, if M is large enough, the minimal set Em sticks at the boundary of the slab: 

Theorem 6.3.2. There exist Mg > 0, Cg > 0, depending on s, such that 
if M ^ Mg then 

(6.3.1) [-1,1) X [CgM"^,M] c Eh 

(6.3.2) and (-1,1] x [-M, -CgM^^] c Em- 
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Figure 23. Stickiness properties of Theorem 6.3.2. 


The situation of Theorem 6.3.2 is described in Figure 23. 

We mention that the “strange” exponent in (6.3.1) and (6.3.2) is optimal. 

For the detailed proof of Theorems 6.3.1 and 6.3.2, and other results on the 
boundary behavior of nonlocal minimal surfaces, see [64]. Here, we limit ourselves 
to give some heuristic motivation and a sketch of the proofs. 

As a motivation for the (somehow unexpected) stickiness property at the bound¬ 
ary, one may look at Figure 22 and argue like this. In the classical case, correspond¬ 
ing to s = Iji, independently on the width 5, the set of minimal perimeter in Bi 
will always be the half-ball Bi n {xn < 0}. 

Now let us take s < 1/2. Then, the half-ball Bi n {xn < 0} cannot be an 
s-minimal set, since the nonlocal mean curvature, for instance, at the origin cannot 
vanish. Indeed, the origin “sees” the complement of the set in a larger proportion 
than the set itself. More precisely, in Bi (or even in Bi+s) the proportion of the set 
is the same as the one of the complement, but outside Bi^s the complement of the 
set is dominant. Therefore, to “compensate” this lack of balance, the s-minimal set 
for s < 1/2 has to bend a bit. Likely, the s-minimal set in this case will have the 
tendency to become slightly convex at the origin, so that, at least nearby, it sees 
a proportion of the set which is larger than the proportion of the complement (we 
recall that, in any case, the proportion of the complement will be larger at infinity, 
so the set needs to compensate at least near the origin). But when S is very small, 
it turns out that this compensation is not sufficient to obtain the desired balance 
between the set and its complement: therefore, the set has to “stick” to the half¬ 
sphere, in order to drop its constrain to satisfy a vanishing nonlocal mean curvature 
equation. 
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Of course some quantitative estimates are needed to make this argument work, 
so we describe the sketch of the rigorous proof of Theorem 6.3.1 as follows. 

Sketch of the proof of Theorem 6.3.1. First of all, one checks that for 
any fixed ry > 0, if <5 > 0 is small enough, we have that the interaction between Bi 
and Bi+s\Bi is smaller than rj. In particular, by comparing with a competitor that 
is empty in Bi, by minimality we obtain that 

(6.3.3) BevsiEs.Bi) 

where we have denoted by Es the s-minimal set in Bi that coincides with Ks 
outside Bi. 

Then, one checks that 

(6.3.4) the boundary of Es can only lie in a small neighborhood of dBi 
if 5 is sufficiently small. 

Indeed, if, by contradiction, there were points of dEs at distance larger than e 
from dBi, then one could find two balls of radius comparable to e, whose centers 
lie at distance larger than e/2 from dBi and at mutual distance smaller than e, and 
such that one ball is entirely contained in Bi n Es and the other ball is entirely 
contained in Bi\Es (this is due to a Clean Ball Condition, see Corollary 4.3 in [ 29 ]). 
As a consequence, Peis{Es, Bi) is bounded from below by the interaction of these 
two balls, which is at least of the order of Then, we obtain a contradiction 

with (6.3.3) (by choosing rj much smaller than and taking 6 sufficiently small). 

This proves (6.3.4). From this, it follows that 

(6.3.5) the whole set Es must lie in a small neighborhood of dBi. 

Indeed, if this were not so, by (6.3.4) the set Es must contain a ball of radius, 
say 1/2. Hence, Pers(iil 5 , Bi) is bounded from below by the interaction of this ball 
against {a;„ > 0}\i3i, which would produce a contribution of order one, which is in 
contradiction with (6.3.3). 

Having proved (6.3.5), one can use it to complete the proof of Theorem 6.3.1 
employing a geometric argument. Namely, one considers the ball Bp, which is out¬ 
side Es for small p > 0, in virtue of (6.3.5), and then enlarges p untill it touches dEs- 
If this contact occurs at some point p e Hi, then the nonlocal mean curvature of Es 
at p must be zero. But this cannot occur (indeed, we know by (6.3.5) that the 
contribution of Es to the nonlocal mean curvature can only come from a small 
neighborhood of dBi, and one can check, by estimating integrals, that this is not 
sufficient to compensate the outer terms in which the complement of Es is domi¬ 
nant). 

As a consequence, no touching point between Bp and dEs can occur in Hi, which 
shows that Es is void inside Hi and completes the proof of Theorem 6.3.1. □ 

As for the proof of Theorem 6.3.2, the main arguments are based on sliding a ball 
of suitably large radius till it touches the set, with careful quantitative estimates. 
Some of the details are as follows (we refer to [64] for the complete arguments). 

Sketch of the proof of Theorem 6.3.2. The first step is to prove a weaker 
form of stickiness as the one claimed in Theorem 6.3.2. Namely, one shows that 

(6.3.6) [-l,l)x[coM,M]c^Eli 

(6.3.7) and (-1,1] x [-M, -CoM] c Em, 
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for some Co e (0,1). Of course, the statements in (6.3.1) and (6.3.2) are stronger 
than the ones in (6.3.6) and (6.3.7) when M is large, since < 1, but we will 
then obtain them later in a second step. 

To prove (6.3.6), one takes balls of radius CqM and centered at {x2 = t}, for 
any t e [cqM, M], One slides these balls from left to right, till one touches BEm- 
When M is large enough (and Cq small enough) this contact point cannot lie 
in {|a;i| < 1}. This is due to the fact that at least the sliding ball lies outside Em, 
and the whole { 3:2 > M) lies outside Em as well. As a consequence, these contact 
points see a proportion of Em smaller than the proportion of the complement (it 
is true that the whole of {x 2 < —M} lies inside Em, but this contribution comes 
from further away than the ones just mentioned, provided that Cq is small enough). 
Therefore, contact points cannot satisfy a vanishing mean curvature equation and 
so they need to lie on the boundary of the domain (of course, careful quantitative 
estimates are necessary here, see [64] , but we hope to have given an intuitive sketch 
of the computations needed). 

In this way, one sees that all the portion [—1,1) x [cqM , M] is clean from the 
set Em and so (6.3.6) is established (and (6.3.7) can be proved similarly). 

Once (6.3.6) and (6.3.7) are established, one uses them to obtain the strongest 
form expressed in (6.3.1) and (6.3.2). For this, by (6.3.6) and (6.3.7), one has only to 

1+25 

take care of points in {\x 2 \ e [CoM ^+^=, CqM]}. For these points, one can use again 
a sliding method, but, instead of balls, one has to use suitable surfaces obtained by 
appropriate portions of balls and adapt the calculations in order to evaluate all the 
contributions arising in this way. 

The computations are not completely obvious (and once again we refer to [64] 
for full details), but the idea is, once again, that contact points that are in the set 

l + 2s 

{1x21 e [CqM 2 + 2 s j CqM]} cannot satisfy the balancing relation prescribed by the 
vanishing nonlocal mean curvature equation. □ 

The stickiness property discussed above also has an interesting consequence 
in terms of the “geometric stability” of the flat s-minimal surfaces. For instance, 
rather surprisingly, the flat lines in the plane are “geometrically unstable” nonlocal 
minimal surfaces, in the sense that an arbitrarily small and compactly supported 
perturbation can produce a stickiness phenomenon at the boundary of the domain. 
Of course, the smaller the perturbation, the smaller the stickiness phenomenon, but 
it is quite relevant that such a stickiness property can occur for arbitrarily small 
(and “nice”) perturbations. This means that s-minimal flat objects, in presence of a 
perturbation, may not only “bend” in the center of the domain, but rather “jump” 
at boundary points as well. 

To state this phenomenon in a mathematical framework, one can consider, for 
fixed i5 > 0 the planar sets 

H :=M. X (— 00 , 0 ), 

F_ := (-3,-2) X [0,<5) 
and F+ := (2,3) X [0,(5). 

One also Axes a set F which contains HuF_uF+ and denotes by E be the s-minimal 
set in (—1,1) X M among all the sets that coincide with F outside (—1,1) x K. 


Then, this set E sticks at the boundary of the domain, according to the next result: 
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Figure 24. The stickiness/instability property in Theo¬ 
rem 6.3.3, with f3 := 


Theorem 6.3.3. Fixeo > 0 arbitrarily small. Then, there exists Sq > 0, possibly 
depending on eg, such that, for any 6 e (0, ^g], 

E 3 (-1,1) X (-00,(5^]. 

The stickiness/instability property in Theorem 6.3.3 is depicted in Figure 24. 
We remark that Theorem 6.3.3 gives a rather precise quantification of the size of the 
stickiness in terms of the size of the perturbation: namely the size of the stickiness 
in Theorem 6.3.3 is larger than the size of the perturbation to the power f3 := 
for any eg > 0 arbitrarily small. Notice that /? ^ -too as s ^ 1/2, consistently with 
the fact that classical minimal surfaces do not stick at the boundary. 

The proof of Theorem 6.3.3 is based on the construction of suitable auxiliary 
barriers. These barriers are used to detach a portion of the set in a neighborhood of 
the origin and their construction relies on some compensations of nonlocal integral 
terms. In a sense, the building blocks of these barriers are “self-sustaining solutions” 
that can be seen as the geometric counterparts of the s-harmonic function dis¬ 
cussed in Section 3.4. 

Indeed, roughly speaking, like the function x\, these barriers “see” a proportion 
of the set in {xi < 0} larger than what is produced by their tangent plane, but a 
proportion smaller than that at infinity, due to their sublinear behavior. Once again, 
the computations needed to check such a balancing conditions are a bit involved, 
and we refer to [64] for the complete details. 


To conclude this chapter, we make a remark on the connection between so¬ 
lutions of the fractional Allen-Cahn equation and s-minimal surfaces. Namely, a 
suitably scaled version of the functional in (5.1.1) T-converges to either the classical 
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perimeter or the nonlocal perimeter functional, depending on the fractional param¬ 
eter s. The T-convergence is a type of convergence of functionals that is compatible 
with the minimization of the energy, and turns out to be very useful when dealing 
with variational problems indexed by a parameter. This notion was introduced by 
De Giorgi, see e.g. [51] for details. 

In the nonlocal case, some care is needed to introduce the “right” scaling of 
the functional, which comes from the dilation invariance of the space coordinates 
and possesses a nontrivial energy in the limit. For this, one takes first the rescaled 
energy functional 

Jg{u,fl) := fl)-f I W{u)dx, 

JQ 

where JC is the kinetic energy defined in (5.1.2). Then, one considers the functional 

' if sG (0,1/2), 


Fs{u,n) 


< |eloge| ^Je{u,U) if s = 1/2, 


ifsG(l/2, 1). 


The limit functional of as e ^ 0 depends on s. Namely, when s G (0,1/2), 
the limit functional is (up to dimensional constants that we neglect) the fractional 
perimeter, i.e. 


(6.3.8) F{u,n) 


JPers(F, n) if u|n = Xe — for some set F c fl 
( - 1-00 otherwise. 


On the other hand, when s G [1/2,1), the limit functional of F/ is (again, up to 
normalizing constants) the classical perimeter, namely 


(6.3.9) F{u,n) 


Per(F, O) if u|n = Xe — Xe^, for some set F c O 
- 1-00 otherwise. 
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That is, the following limit statement holds true: 

Theorem 6.3.4. Let s e (0, 1). Then, T-converges to F, as defined in 
either (6.3.8) or (6.3.9), depending on whether s G (0,1/2) or s G [1/2,1). 

For precise statements and further details, see [ 123 ]. Additionally, we remark 
that the level sets of the minimizers of the functional in (5.1.1), after a homogeneous 
scaling in the space variables, converge locally uniformly to minimizers either of the 
fractional perimeter (if s G (0,1/2)) or of the classical perimeter (if s G [1/2,1)): 
that is, the “functional” convergence stated in Theorem 6.3.4 has also a “geometric” 
counterpart: for this, see Corollary 1.7 in [ 125 ]. 

One can also interpret Theorem 6.3.4 by saying that a nonlocal phase transition 
possesses two parameters, e and s. When g —> 0, the limit interface approaches a 
minimal surface either in the fractional case (when s G (0,1/2)) or in the classical 
case (when s G [1/2,1)). This bifurcation at s = 1/2 somehow states that for lower 
values of s the nonlocal phase transition possesses a nonlocal interface in the limit, 
but for larger values of s the limit interface is characterized only by local features (in 
a sense, when s G (0,1/2) the “surface tension effect” is nonlocal, but for s G [1/2,1) 
this effect localizes). 

It is also interesting to compare Theorems 6.0.5 and 6.3.4, since the bifurcation 
at s = 1/2 detected by Theorem 6.3.4 is perfectly compatible with the limit behavior 
of the fractional perimeter, which reduces to the classical perimeter exactly for this 
value of s, as stated in Theorem 6.0.5. 



CHAPTER 7 


A nonlocal nonlinear stationary Schrodinger type 

equation 


The type of problems introduced in this chapter are connected to solitary solu¬ 
tions of nonlinear dispersive wave equations (such as the Benjamin-Ono equation, 
the Benjamin-Bona-Mahony equation and the fractional Schrodinger equation). In 
this chapter, only stationary equations are studied and we redirect the reader to 
[ 137 , 138 ] for the study of evolutionary type equations. 

Let n > 2 be the dimension of the reference space, s e (0,1) be the fractional 
parameter, and e > 0 be a small parameter. We consider the so-called fractional 
Sobolev exponent 

{ 2n 

-— for n ^ 3, or n = 2 and s e (0,1/2) 

77 / 

-l-cx) for n = 1 and s e (0,1/2] 


and introduce the following nonlocal nonlinear Schrodinger equation 


(7.0.10) 


je^^{-Ayu + u = uP inflcK" 
= 0 in R"\n, 


in the subcritical case p e (1, 2* — 1), namely when pell, - ). 

\ n — 2s J 

This equation arises in the study of the fractional Schrodinger equation when 
looking for standing waves. Namely, the fractional Schrodinger equation considers 
solutions T = T(x, t) : M" X K ^ C of 

(7.0.11) ihdt'f = + P)T, 


where s e (0,1), h is the reduced Planck constant and V = P(a;, t, [Tl) is a po¬ 
tential. This equation is of interest in quantum mechanics (see e.g. [102] and the 
appendix in [ 47 ] for details and physical motivations). Roughly speaking, the quan¬ 
tity |'k(a;, t)p dx represents the probability density of finding a quantum particle in 
the space region dx and at time t. 

The simplest solutions of (7.0.11) are the ones for which this probability density 
is independent of time, i.e. |T(a;,f)| = u(x) for some u : M" —> [0, -too). In this way, 
one can write T as u times a phase that oscillates (very rapidly) in time: that is 
one may look for solutions of (7.0.11) of the form 


^(x,t) := u(a:)e“*/^, 

for some frequency w e K. Choosing V = RdTj) = —= —uP~^, a substitu¬ 
tion into (7.0.11) gives that 

(^n^y-Ayu + uju- = h^y-Ay^^ - ihdt^ + pt = o, 

which is (7.0.10) (with the normalization convention a; := 1 and e := K). 
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The goal of this chapter is to construct solutions of problem (7.0.10) that con¬ 
centrate at interior points of the domain for sufficiently small values of e. We 
perform a blow-up of the domain, defined as 

a; e flj. 

We can also rescale the solution of (7.0.10) on flj, 

Ue{x) = u{ex). 



The problem (7.0.10) for then reads 


(7.0.12) 


j (—A)®it + u = uP in flj 

= 0 in IR”\r2g. 


When e ^ 0, the domain invades the whole of the space. Therefore, it is also 
natural to consider (as a first approximation) the equation on the entire space 

(7.0.13) {-Ayu + u = uP 


The first result that we need is that there exists an entire positive radial least energy 
solution w e i7®(M"') of (7.0.13), called the ground state solution. Here follow some 
relevant results on this. The interested reader can find their proofs in [ 84 ]. 

(1) The ground state solution w e i7'^(M") is unique up to translations. 

(2) The ground state solution w e i7®(M"’) is nondegenerate, i.e., the deriva¬ 
tives DiW are solutions to the linearized equation 

(7.0.14) {-Ayz + Z = pZP-^. 

(3) The ground state solution w e 77® (M”) decays polynomially at infinity, 
namely there exist two constants a, /3 > 0 such that 

a|x|-("+2®) ^ u{x) ^ /3|x|-("+2®). 

Unlike the fractional case, we remark that for the (classical) Laplacian, at infinity 
the ground state solution decays exponentially fast. We also refer to [ 83 ] for the 
one-dimensional case. 

The main theorem of this chapter establishes the existence of a solution that 
concentrates at interior points of the domain for sufficiently small values of e. 
This concentration phenomena is written in terms of the ground state solution w. 
Namely, the first approximation for the solution is exactly the ground state w, scaled 
and concentrated at an appropriate point ^ of the domain. More precisely, we have: 


Theorem 7.0.5. If e is sufficiently small, there exist a point f e H and a 
solution Ue of the problem (7.0.10) such that 



^ C'e”+2®, 


and dist{f, dn) ^ S > 0. Here, C and 5 are constants independent of e or U, and 
the function w is the ground state solution of problem (7.0.13). 


The concentration point ^ in Theorem 7.0.5 is influenced by the global geom¬ 
etry of the domain. On the one hand, when s = 1, the point f is the one that 
maximizes the distance from the boundary. On the other hand, when s e (0,1), 
such simple characterization of f does not hold anymore: in this case, ^ turns out 
to be asymptotically the maximum of a (complicated, but rather explicit) nonlocal 
functional: see [ 47 ] for more details. 
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We state here the basic idea of the proof of Theorem 7.0.5 (we refer again to 


[47] for more details). 

Sketch of the proof of Theorem 7.0.5. In this proof, we make use of the 
Lyapunov-Schmidt procedure. Namely, rather than looking for the solution in an 
infinite-dimensional functional space, one decomposes the problem into two orthog¬ 
onal subproblems. One of these problem is still infinite-dimensional, but it has 
the advantage to bifurcate from a known object (in this case, a translation of the 
ground state). Solving this auxiliary subproblem does not provide a true solution 
of the original problem, since a leftover in the orthogonal direction may remain. To 
kill this remainder term, one solves a second subproblem, which turns out to be 
finite-dimensional (in our case, this subproblem is set in K”, which corresponds to 
the action of the translations on the ground state). 

A structural advantage of the problem considered lies in its variational struc¬ 
ture. Indeed, equation (7.0.12) is the Euler-Lagrange equation of the energy func¬ 
tional 



(7.0.15) 


for any u e H§{n^) := {u e i7®(K") s.t. u = 0 a.e. in K"\r2£}. Therefore, the 
problem reduces to finding critical points of A. 

To this goal, we consider the ground state solution w and for any ^ e K" we 
let := w{x — For a given ^ e fig a first approximation for the solution of 
problem (7.0.12) can be taken as the solution of the linear problem 



(7.0.16) 


The actual solution will be obtained as a small perturbation of uj for a suitable 


point ^ = ^(e). We define the operator C := (—A)'^ -I- 1, where I is the identity and 
we notice that C has a unique fundamental solution that solves 


£r = 6o in K”. 

The Green function Gg of the operator C in fig satisfies 



CGt{x, y) = 5y{x) if T e fig, 

Gg(ai,j/) = 0 if T e IR"'\flg. 


(7.0.17) 


It is convenient to introduce the regular part of Gg, which is often called the Robin 
function. This function is defined by 


(7.0.18) 


He{x, y) := r(a; - y) - Gg(x, y) 


and it satisfies, for a fixed y e K” 


(7.0.19) 



Then 




and by (7.0.17) 
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The operator £ is self-adjoint and thanks to the above identity and to equation 
(7.0.16) it follows that 


u^{x) = Cu^{y)G^{x,y) dy 

Jn, 

= f w^Ay)G,{x,y)dy. 

Jn, 

So, we use (7.0.18) and we obtain that 

iZ 5 (t) = r inf (i/)r(x -y)dy-\ inf (j/)i7.(T, y) dy. 

Jn. Jn. 


Now we notice that, since is solution of (7.0.13) and T is the fundamental solution 
of £, we have that 


wAy)T{x-y)dy 

Jk" 



£w^{y)T{x - y) dy 
■Wf{y)Cr{x - y) dy 


w^{x). 


Therefore we have obtained that 

(7.0.20) U 5 (a:) = in 5 (a;) - f wl{y)T{x - y) dy - \ wAy)HAx,y) dy. 

jR"\n, Jn, 

Now we can insert (7.0.20) into the energy functional (7.0.15) and expand the errors 

in powers of e. For dist(^,(5ne) > - with S fixed and small, the energy of is a 

perturbation of the energy of the ground state in and one finds (see Theorem 4.1 
in [47]) that 

(7.0.21) lAu^) = I(w) + + o(e”+"^), 

where 

■He(0 := f f HAx,y)wAx)wAy)dxdy 
Jn, Jn, 

and / is the energy computed on the whole space M", namely 

(7.0.22) I{u) = - r (^(—A)'^it(x)-I-u(a;)')it(a:) da;- f u^~^^{x)dx. 

2 Jar. V J p -t 1 Jr„ 

In particular, £( 115 ) agrees with a constant (the term /(in)), plus a functional over 
a finite-dimensional space (the term HA^), which only depends on ^ e K"), plus a 
small error. 

We remark that the solution of equation (7.0.16) which can be obtained 
from (7.0.20) does not provide a solution for the original problem (7.0.12) (indeed, 
it only solves (7.0.16)): for this, we look for solutions of (7.0.12) as perturbations 
of u^, in the form 

(7.0.23) itj := + A- 


The perturbation functions A ar® considered among those vanishing outside fie and 
orthogonal to the space Z = Span(Zi,..., Z„), where Zi = are solutions of 


the linearized equation (7.0.14). This procedure somehow “removes the degeneracy”, 
namely we look for the corrector 1 /; in a set where the linearized operator is invertible. 
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This makes it possible, fixed any ^ e K", to find ip = such that the function u^, 
as defined in (7.0.23) solves the equation 

n 

(7.0.24) (-A)®M{ + uj = 2 CjZj in 

i = l 

That is, is solution of the original equation (7.0.12), up to an error that lies 
in the tangent space of the translations (this error is exactly the price that we 
pay in order to solve the corrector equation for ip on the orthogonal of the kernel, 
where the operator is nondegenerate). As a matter of fact (see Theorem 7.6 in [47] 
for details) one can see that the corrector ip = ip^ is of order £”+ 2 ^*. Therefore, 
one can compute leiu^) = A(w{ + V'^) ^ ^ higher order perturbation of Ie(u^). 
From (7.0.21), one obtains that 

(7.0.25) h{u^) = I{w) + ^77.(0 + 0(e”+4^), 

see Theorem 7.17 in [47] for details. 

Since this energy expansion now depends only on ^ G K", it is convenient to 
define the operator ^ K as 

JeiO 

This functional is often called the reduced energy functional. From (7.0.25), we 
conclude that 

(7.0.26) J,(0 = I{w) + + 0(e”+"^). 

The reduced energy J plays an important role in this framework since critical points 
of J correspond to true solutions of the original equation (7.0.12). More precisely 
(see Lemma 7.16 in [47]) one has that Ci = 0 for all i = 1,... ,n in (7.0.24) if and 
only if 

(7.0.27) 

In other words, when e approaches 0, to find concentration points, it is enough to 
find critical points of J, which is a finite-dimensional problem. Also, critical points 
for J come from critical points of "He, up to higher orders, thanks to (7.0.26). The 
issue is thus to prove that T-Le does possess critical points and that these critical 
points survive after the small error of size in fact, we show that "He possesses 

a minimum, which is stable for perturbations. For this, one needs a bound for the 
Robin function from above and below. To this goal, one builds a barrier function 
defined for ^ G fig and x G K" as 

j3^{x) := j r{z — ^)T{x — z) dz. 

Using this function in combination with suitable maximum principles, one obtains 
the existence of a constant c G (0,1) such that 

cH^{x,^) ^ I3^{x) < c~^H^{x,i), 

for any x G R" and any ^ e with dist(^, dfl^) > 1, see Lemma 2.1 in [47] . From 
this it follows that 


( 7 . 0 . 28 ) 


77,(0 - 
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for all points ^ G fie such that d e [5, ^/e]. So, one considers the domain of the 
points of rig that lie at distance more than 5/e from the boundary of By (7.0.28), 
we have that 

(7.0.29) for any ^ e 

Also, up to a translation, we may suppose that 0 G fl. Thus, 0 G and its distance 
from drig is of order 1/e (independently of (5). In particular, if 6 is small enough, 
we have that 0 lies in the interior of and (7.0.28) gives that 

By comparing this with (7.0.29), we see that 77e has an interior minimum in fie/. 
The value attained at this minimum is of order and the values attained at the 

boundary of are of order ^n+4s^ which is much larger than j£ g -g 

small enough. This says that the interior minimum of He in ^e,s is nondegenerate 
and it survives to any perturbation of order ^ g -g gjjjall enough. 

This and (7.0.26) imply that J has also an interior minimum at some point ^ 
in By construction, this point ^ satisfies (7.0.27), and so this completes the 
proof of Theorem 7.0.5. □ 

The variational argument in the proof above (see in particular (7.0.27)) has a 
classical and neat geometric interpretation. Namely, the “unperturbed” functional 
(i.e. the one with e = 0) has a very degenerate geometry, since it has a whole 
manifold of minimizers with the same energy: this manifold corresponds to the 
translation of the ground state w, namely it is of the form Mg := {w^, ^ G M"} and, 
therefore, it can be identified with K". 


Mo 



Figure 26. Geometric interpretation 


For topological arguments, this degenerate picture may constitute a serious 
obstacle to the existence of critical points for the “perturbed” functional (i.e. the 
one with e ¥= 0). As an obvious example, the reader may think of the function of 
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two variables /<. : —> M given by fe{x,y) := + ey. When e = 0, this function 

attains its minimum along the manifold {x = 0}, but all the critical points on this 
manifold are “destroyed” by the perturbation when e ^ 0 (indeed V/e(x, y) = {2x, e) 
never vanishes). 

In the situation described in the proof of Theorem 7.0.5, this pathology does not 
occur, thanks to the nondegeneracy provided in [84] . Indeed, by the nondegeneracy 
of the unperturbed critical manifold, when e ^ 0 one can construct a manifold, 
diffeomorphic to the original one (in our case of the form := {u^+tpiOj ? ^ ®"})j 
that enjoys the special feature of “almost annihilating” the gradient of the functional, 
up to vectors parallel to the original manifold Mq (this is the meaning of formula 
(7.0.24)). 

Then, if one finds a minimum of the functional constrained to M^, the theory of 
Lagrange multipliers (at least in finite dimension) would suggest that the gradient 
is normal to Mg. That is, the gradient of the functional is, simultaneously, parallel 
to Mq and orthogonal to Mg. But since Mg is diffeomorphically close to Mq, the 
only vector with this property is the null vector, hence this argument provides the 
desired critical point. 


We also recall that the fractional Schrodinger equation is related to a nonlocal 
canonical quantization, which in turn produces a nonlocal Uncertainty Principle. 
In the classical setting, one considers the momentum/position operators, which are 
defined in K", by 


(7.0.30) Pk ■■= -ihdk and Qfc := Xk 


for k e {1,... ,n}. Then, the Uncertainty Principle states that the operators P = 
(Pi,..., Pn) and Q = (Qi,..., Q„) do not commute (which makes it practically 
impossible to measure simultaneously both momentum and position). Indeed, in 
this case a simple computation shows that 


(7.0.31) [Q, P] := 2 [Qfc, Pk] = ihn. 

k = l 


The nonlocal analogue of this quantization may be formulated by introducing a 
nonlocal momentum, i.e. by replacing the operators in (7.0.30) by 


(7.0.32) 


P^ := —ih^dk{—^) 2 and Qfe := Xk- 
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In this case, using that the Fourier transform of the product is the convolution of 
the Fourier transforms, one has that 

{xk*g){0 = ^{xk^~^9{x)){0 

= f dxl dy 

= dx\ 

Zm J^n 

(7.0.33) = ^ dy dkg{y) 

ZTT Ju^Ti 

for any test function g. In addition, 

= {2nrh^^k ier-7- 

Therefore, given any test function if}, using this with f := ijj and / := Tfc'i/', and 
also (7.0.33) with g := P{P^tjj) and g := ip, we obtain that 

^{QkP^iPix) - P^QkHx)) 

= ^{xkPk'4’{x)- Pk{Xk^{x))) 

= Xk * PiPktpix)) - P{P^{xki^{x))) 

= {2Try~^ifi^dk{P,k - (27r)®fi®^fe * '0(6 

= (27r)^-izr4(6 lcr-'0(6) - i2nr-hn^ik l6^-'40(6 
= {27Ty-HH^dk{ikyr^)m 
= {2ny-yn^ (16*-' + (s - ml icr-") 0(6- 

Consequently, by summing up, 

P{[Q, pyy) = (27r)^-izfi^ 16*-' (n + s - 1) 0(6. 

So, by taking the anti-transform, 

^ [Q, ny = 7fi* (n + s - 1) .F-' ((2^|6)*-'0) 

= ih'^ {n + s — 1) (—A)t^ 0. 

Notice that, as s ^ 1, this formula reduces to the the classical Heisenberg Uncer¬ 
tainty Principle in (7.0.31). 


7.0.1. Fhom the nonlocal Uncertainty Principle to a fractional weighted 
inequality. Now we point out a simple consequence of the Uncertainty Principle 
in formula (7.0.34), which can be seen as a fractional Sobolev inequality in weighted 
spaces. The result (which boils down to known formulas as s ^ 1) is the following: 
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Proposition 7.0.6. For any u e 5(]R"), we have that 
(—A)~ u 

Proof. The proof is a general argument in operator theory. Indeed, suppose 
that there are two operators S and A, acting on a space with a scalar Hermitian 
product. Assume that S is self-adjoint and A is anti-self-adjoint, i.e. 

(u, Su} = {Su, u) and (u, Au) = —{Aw, u), 
for any u in the space. Then, for any A G K, 

IKA+AS'juf = IIAuf-f A^ll^uf+ a(<Au,S'u> + <S'u,Au>) 

= \\Auf + A^IIS-up -f A<(5'A - AS)u,u}. 

Now we apply this identity in the space c taking S := Qk = Xk 

and A := iP^ = A)“ 2 “ (recall (7.0.32) and notice that iP^ is anti-self- 

adjoint, thanks to the integration by parts formula). In this way, and using (7.0.34), 
we obtain that 






77, -I- .S — 1 


7,2 (-Rn I 


V(-A)V-u 


r,27Rn'l 




k=l 


Xj [ll'Pfc^lli2(K") + "^^IIQfe^lli2(B") -f 7A<[Qfc, 


k=l 


V{-A)^u 


L 2 (R>i) 


+ A^ 


\x\u 


-1-7^ A (n-I-s — 1) fi®i((—A) 2 


= V(-A)^u 


H-A" 


\x\u 


L2(R") 
) 

2 

L2(R»i) 


-A(n + s —I)fi® (—A) 


0 


Now, if u ^ 0, we can optimize this identity by choosing 


(n-I-s — 1) fi® (—A) ’4 u 


A := - 

and we obtain that 

0 ^ V(-A)^7i 

which gives the desired result. 


L2(R»i) 


\x\u 


2 

L2(R»i) 


(n-I-s — 1)^ (—A) 4 u 


L2(R»i) 


\x\u 


L2(Rn) 


□ 

















































APPENDIX A 


Alternative proofs of some results 


A.l. Another proof of Theorem 3.4.1 

Here we present a different proof of Theorem 3.4.1, based on the Fourier trans¬ 
forms of homogeneous distributions. This proof is the outcome of a pleasant dis¬ 
cussion with Alexander Nazarov. 


Alternative proof of Theorem 3.4.1. We are going to use the Fourier 
transform of \xY‘ in the sense of distribution, with q e C\Z. Namely (see e.g. 
Lemma 2.23 on page 38 of [97]) 

(A.1.1) = 

with 

(A.1.2) Cq := -2(27r)-«-ir(l -f q) sin y. 

We remark that the original value of the constant Cq in [97] is here multiplied by 
a (27r)“'^“^ term, in order to be consistent with the Fourier normalization that we 
have introduced. We observe that the function K 3 cc i-^ is locally integrable 
only when g > — 1, so it naturally induces a distribution only in this range of the 
parameter q (and, similarly, the function K 3 ^ is locally integrable only 

when q < 0): therefore, to make sense of the formulas above in a wider range of 
parameters q it is necessary to use analytic continuation and a special procedure 
that is called regularization: see e.g. page 36 in [97] (as a matter of fact, we will do a 
procedure of this type explicitly in a particular case in (A.2. 7)). Since K 3 x i-^ |x|® 
is even, we can write (A.1.1) also as 

(A.1.3) T-\m = Cq\x\-^-‘^. 


We observe that, by elementary trigonometry, 

. 7r(s + 1) . 7r(s — 1) , 7r(s — 2) . tts 

sin- = — sin- and sm- = sm —. 

2 2 2 2 

Moreover, 


r(2 4- s) = (1 -t s)r(l 4- s) and r(s) = (s - l)r(s - 1). 
Hence, by (A. 1.2), 

(A.1.4) 


^ • Cs+i Cs -2 = ^ • 4(2^)-2*-ir(2 + s) r(s - 1) sin 
1 4- s 1 4- s 2 

= 4r(l 4 - s) r(s) sin ^ sin y = C's Cs-i- 

Moreover, 

|x|® 4-= 2x+. 

s 4- 1 


sin 


7r(5 — 2) 
2 
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So, taking the Fourier transform and using (A.1.1) with q := s and q := s + 1, we 
obtain that 

2j-«) = :F{\xn + j^:F{dM^+^) 

s + 1 

As a consequence, 

= a + ^Cs+i ici-2+^ 

Hence, recalling (7.0.33), 

= a j-i(ici-i+^) - • f 

s + 1 Ztt 

= C..F-I(|e|-1+^) + ^5,.F-i(|e|-2+^). 

Accordingly, exploiting (A.1.3) with q := —1 + s and q := —2 + s, 
2.F-i(|Cp*J-(x^)) = CsC,-i\x\-^ 

= Cs a_i |a;|-^ + ^ • Cs+i Cs-2 X 
1 + s 

So, recalling (A. 1.4), 

2J-i(|Cp*J'(x^)) =CsCs-i (|t|-^-x|x|-i-'*). 

This and (3.1.7) give that 

(-A)««) = (7,(N-«-cr|xri-^), 

for some Cs- In particular, the quantity above vanishes when x > 0, thus providing 
an alternative proof of Theorem 3.4.1. □ 

Yet another proof of Theorem 3.4.1 can be obtained in a rather short, but 
technically quite advanced way, using the Paley-Wiener theory in the distributional 
setting. The sketch of this proof goes as follows: 

Alternative proof of Theorem 3.4.1. The function h := x^ is homoge¬ 
neous of degree s. Therefore its (distributional) Fourier transform J-h is homoge¬ 
neous of degree — 1 — s (see Lemma 2.21 in [97]). 

Moreover, h is supported in {x > 0}, therefore J^h can be continued to an 
analytic function in C_ := {z 6 C s.t. Sz < 0} (see Theorem 2 in [129]). 

Therefore, g{^) := |^|®J^/i(^) is homogeneous of degree —1 and is the trace of a 
function that is analytic in C_. 

In particular, for any y < 0, we have that 
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where c := That is, g{z) coincides with ^ on a half-line, and then in the 

whole of C_, by analytic continuation. 

That is, in the sense of distributions. 


9{i) 


c 

^-iO’ 


for any ^ e K. 

Now we recall the Sokhotski Formula (see e.g. (3.10) in [17]), according to which 


where 5 is the Dirac’s Delta and the identity holds in the sense of distributions. 
By considering this equation with the two sign choices and then summing up, we 
obtain that 

= — 2z7r5. 


As a consequence 


iO 

9{i) = 




c - iO C + *0 


-I- 2icn5. 


Therefore 

= ier-7(a = + 2^c7^|er^. 

Of course, as a distribution, = 0, since the evaluation of at ^ = 0 vanishes, 
therefore we can write 

m 

Accordingly, £ is homogeneous of degree — 1 + s and it is the trace of a function 
analytic in C+ := {^ G C s.t. '^z > 0}. 

Consequently, is homogeneous of degree —s (see again Lemma 2.21 in [97]), 

and it is supported in {a:: ^ 0} (see again Theorem 2 in [129]). Since (—A)®a;+ 
coincides (up to multiplicative constants) with J-~^£, we have just shown that 
(—A)®a;^ = CoxZ’^, for some Co G K, and so in particular (—A)®a:^ vanishes 
in {a; > 0}. □ 


A.2. Another proof of Lemma 3.1.3 

For completeness, we provide here an alternative proof of Lemma 3.1.3 that 
does not use the theory of the fractional Laplacian. 


Alternative proof of Lemma 3.1.3. We first recall some basic properties 
of the modified Bessel functions (see e.g. [4]). First of all (see formula 9.1.10 on 
page 360 of [4]), we have that 


+ CX) 


57 L 


(- 1 ) 


k^2k 


2 S Aj 22^ k\ r(s + k + 1) 2" r(l + s) 

k=0 

as 1^1 ^ 0. Therefore (see formula 9.6.3 on page 375 of [4]), 

Is{z) := J^(e^z) 


+ 0{\zY+^) 


= e 2 


.2«r(i -t s) 


+o(izr«)) 


2®r(i -f s) 


+ 0{\z\^+^), 
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as \z\ 0. Using this and formula 9.6.2 on page 375 of [4], 


Ks{z) := 


2 sin( 7 rs) 

TT 


(/-(^)-4L)) 


2 sin( 7 rs) V2-^r(l-s) 2«r(l + s) 

Thus, recalling Euler’s reflection formula 

r{i-s)r{s)= ^ 


+ oi\zrn ■ 


sin(7rs)' 


and the relation r(l + s) = sr(s), we obtain 


Ksiz) = 


r(i-s)r(s) 


2-^r(i-s) 2"r(i + s) 


+ o(i^r-^) 


r(s)z-'’ r(i-s)0" 


+ o(|^p-«), 


21-s 2i+®s 

as \z\ 0. We use this and formula (3.100) in [106] (or page 6 in [ 112 ]) and get 

that, for any small a > 0 , 

cos(27rt) 


I 

Jo 


TT 2 

-Y dt = » ^ ^ ^ -p^Xs(27ra) 


0 (t2 + a2)®+5 a®r(s+5) 


(A.2.1) 


’r(s+i) [271 


r(s) r(l-s)7r*a^ 


2 s 


+ 0 (a"-") 




2 a 2 '^r(s+i) 2 sr(s+i) 

Now we recall the generalized hypergeometric functions mFn (see e.g. page 211 
in [ 112 ]): as a matter of fact, we just need that for any b, c, d > 0 , 

i_F 2(^; c, a, 0) r'/M 

r{b) r(c)r((i) 

We also recall the Beta function relation 


(A.2.2) 


B(a,/3) = 


r(a)r(/3) 


r(a + /3) ’ 

see e.g. formula 6.2.2 in [4[. Therefore using formula (3.101) in [106] (here with y := 
0 , u := 0 and /r := s — 5 , or see page 10 in [ 112 ]), 


r+ac 

Jo 


dt 


-2s 


0 ( a 2 + t ^)®+2 


^ { A’ ® 1 1^2 a; 1 ~ 0 


r(i) r(s) 


2 a 2 «r(l + s)' 

Then, we recall that T (|) = tt^, so, making use of (A.2.1), for any small a > 0, 
1 - cos(27rt) r(l - s) 


J 

•JO 


-Y dt = 




io + 2 sr(s+i) 

Therefore, sending a —> 0 by the Dominated Convergence Theorem we obtain 
1 — cos(27rt) 


(A.2.3) f 
Jo 


tl + 2s 


■ dt = lim 


I"'"” 1 — cos(27rt) 7 r 2®+2 r(l — s) 
oJo (t 2 + a 2 )-+i * ^ 2 sr (s + i) 
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Now we recall the integral representation of the Beta function (see e.g. formu¬ 
las 6.2.1 and 6.2.2 in [4]), namely 


(A.2.4) 


r(^)r(i-fs) 
r(t + 5) 


= B 


n — 1 1 

2+''* = 


p + 00 

Jo 


n-3 

r 2 


(1 + r)^ 


dr. 


We also observe that in any dimension N the {N — l)-dimensional measure of the 

K 

unit sphere is ^ x , (see e.g. [90]). Therefore 


j 

JR- 


dY 


Ntt 2 


p + CX) 

Jo 


oAf-1 


■ dp. 


I / -I I T 7-1 0\ ^ + l + 2s -pi / jV I 1 ^ I / -1 O \ ^ + 

(1 -f \Y\^) 2 r (-^ -f ij Jo (1 4 - p^) 2 

In particular, taking N := n — 1 and using the change of variable p^ =■ t, 


J 

Jr"-i 


dp 


(n — 1) TT 2 


-1 (1+ |^|2)=V 


r( 

c+co 


- + 1) Jo 


2r(^ 


I TT 2 p 

- + 1) Jo 


71-3 

r 2 


l) Jo (1 4- p2) 

dr. 


dp 


l) Jo (1 4- r) 

By comparing this with (A.2.4), we conclude that 

dp _ (n-l)7rV r(^^)r(i4-s) 


J 


-1 (l + |ryP)"^ 


2r(^ + i) r(t + s) 


r(f+ s) 


Accordingly, with the change of variable p := |a;i| ^(w 2 ,..., w„), 
1 — cos(27rwi) 


J 

JW 


■ du! 


1 — cos(27ra;i) 


|^|n+2s 
JR \JR' 

JI 

JR \JR'' 

7r^r(i4-s) r 1 - cos(27rwi) 


‘-1 (wf 4- a;| 4- • • • 4- w^)" 

1 — cos(27ra;i) , \ , 
- „ , dll I duji 

|wi|i+2^(l-f |ry|2)^ J 


duJ2 ■ ■ ■ duJn dull 


r{f+ •) 


J 

Jr 


|a;i|i+2^ 

27r^ r (i 4- s) 1 — cos(27rt) 


J 

Jo 


r(f + s) Jo 

Hence, recalling (A.2. 3), 

1 — cos(27ra;i) 


^l + 2s 


duji 

dt. 


\ 

Je' 


|n+2s 


duj = 


27rVr(i-4s) 7r2®+5r(l-s) 


r(f 4 -s) 2sr(s + i) 

7rt+2sr(l-s) 

^r(t + 5) ’ 


as desired. 


□ 
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To complete the picture, we also give a different proof of (A. 2. 3) which does 
not use the theory of special functions, but the Fourier transform in the sense of 
distributions and the unique analytic continuation: 

Alternative proof of (A. 2. 3). We recall (see^ e.g. pages 156-157 in [141]) 
that for any A G C with SftA G (0,1/2) we have 




T(^) kA 


in the sense of distribution, that is 


r |a:|^ ^${x)dx= ^ f |a;| ^(t){x)dx 

Jr J dR 

for every G C® (M). The same result is obtained in [23] , in the proof of Proposition 
2.4 b), pages 13-16. 

As a consequence 


(A.2.5) 


J(J 


\x\^-^(l>{y)e-'^^^^y dy] dx = 




a:| (/(x) dx. 


By changing variable x —x in the first integral, we also see that 
(A.2.6) dy'^ dx = - \x\-^(j){x) dx. 


By summing together (A.2.5) and (A.2.6), we obtain 

f f f cos{2TTxy) dy') dx = ^ f \x\-^(j){x) dx. 

Jr \Jr / 1 l“2“j 

It is convenient to exchange the names of the integration variables in the first 
integral above: hence we write 

f f f |y|^"V(a;) cos{2TTxy) dx') dy = ^ f \x\-^(j){x) dx. 

JR \JR / i [^-) JR 

Now we fix i? > 0, we point out that 

\ \y\^~^Xi-R,R){y)dy = 2 { y^-^dy = 

Jr Jo 

and we obtain that 


2R^ 


(A.2.7) 


( lyl^ ^<l}{x){cos{2TTxy)-Xi-R,R){y))dx) dy 
JR \JR / 


r(i^) Jr 


x\ ^<j){x)dx — 


2R^ 


4>{x) dx. 

Jr 


This formula holds true, in principle, for A G C, with IftA G (0,1/2), but by the 
uniqueness of the analytic continuation in the variable A it holds also for A G (—2, 0). 
Now we observe that the map 

{x,y) ^ |2/|^"V(a;) (cos(27ra;j/) -X(-R,R)(y)) 


^To check (A.2.5) the reader should note that the normalization of the Fourier transform on 
page 138 in [141] is different than the one here in (3.1.1). 
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belongs to L^(]R xK) if A G (—2,0), hence we can use Fubini’s Theorem and exchange 
the order of the repeated integrals in (A.2. 7): in this way, we deduce that 

( \y\^~^<l}{x)icos{2TTxy)-X{-R.R)iy))dy) dx 

JR \JR / 




r(^] 


J 

Jk 


2i?^ 


x\ ^4>{x)dx -— 

A 


4>{x) dx. 

Jr 


Since this is valid for every (j) G C® (M), we conclude that 

(A.2.8) r \y\^-^ {cos{2TTxy) - Xi-R,R){y)) ^y = - 

Jm i ) 

for any A G (—2,0) and x ^ 0. In this setting, we also have that 


2R^ 


r+oo 

I 

JR 


f \y\^ ^ {x(-R,R}{y)-^)dy =-2 

JR 

By summing this with (A.2.8), we obtain 

r (cos(27ra;y) - 1) dy = 

Jr 

Hence, we take A := —2s G (—2,0), and we obtain that 


dy = 


2R^ 


r(i^) 




1 . 


cos(27rxy) — 1 


dy = 


_2s+ „ 


r(-s), 


r(i + s) ' ' ■ 

By changing variable of integration t := xy, we obtain (A. 2. 3). 


□ 
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